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This  work  deals  with  the  synthesis  of  feedback  systems  to  achieve 


strained  part  of  the  system,  denoted  as  the  'plant'.  Part  of  this  work 


deals  with  linear  time-invariant  (lti)  plants  where  the  'cost  of  feed 


back 


if  lti  compensation  is  used,  is  primarily  in  the  bandwidth  of 


the  feedback  loop  being  much  larger  than  that  of  the  system  as  a whole 


making  the  system  very  sensitive  to  sensor  noise.  Here,  the  objective 


is  to  reduce  the  loop  bandwidth  by  means  of  non-lti  compensation.  A 


nonlinear  first-order  reset  element  (FORE)  is  introduced  and  with  it 


a quantitative  synthesis  procedure  which  permits  design  to  specifications 


despite  large  but  bounded  plant  uncertainty.  The  result  is  a very 


significant  reduction  in  loop  bandwidth  and  with  it,  system  sensitivity 


to  sensor  noise.  Stability  criteria  are  included  which  helps  generalize 


Another  method  of  non-lti  synthesis  is  by  means  of  linear  time 


varying  compensation,  applicable  to  a certain  problem  class.  The 


solution  is  not  in  general  available  analytically,  but  is  found  for 


certain  cases  and  exhibits  in  these  reduced  'cost  of  feedback 


final  part  of  this  work,  the  plant  can  be  nonlinear  with  large  uncer 


tainty  giving  a set  of  nonlinear  plants  W . The  concept  is  use  of 


a lti  set  P which  is  precisely  equivalent  to  W providing  the  output 


is,  in  both  cases,  a member  of  an  acceptable  output  set  A 


thesis  procedure  is  presented,  based  on  this  concept,  for  achieving 


time-domain  specifications  of  a specified  nonlinear  character.  All  of 


the  methods  are  illustrated  by  detailed  design  examples 
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1. 1 Generalities 

An  important  problem  area  in  systems  theory  is  to  satisfy 
quantitative  specifications  despite  parameter  uncertainty.  It  is 
assumed  here  that  there  is  a single-input-single-output  constrained 
part  denoted  as  the  'plant'  and  described  by  a linear  differential 
equation  with  fixed  coefficients.  (Fig.  1.1) 


input 


Pf output  yinj+a.y'n  ~^+ — -fa  y = * 

-1  y 1 n n > rr 

Figure!. 1.  Plant  characterization 


tlBl+b1xtB"11. 


There  exists  a finite  set  of  parameters  {k^kj, . . . ,k^}  such 

that  each  coefficient  is  a function  of  those  parameters,  i.e. 

“ (J^dCj^kj,...,^)  for  i = l,n  , b^  = (l^.k^ . . . .k^)  for 

j»l,m  and  K = cfetk^jk^, . . . ,kf)  . Furthermore,  it  is  assumed  that 

each  parameter  k^  , i=l,  l , is  associated  with  a known  range  of 

uncertainty  (k^  ,k^  ] which  can  be  small  or  large.  We  want  to 

min  max 

guarantee  time  domain  specifications  (T.D.S.),  for  example  unit  step 
response  of  the  type  shown  in  Fig.  1.2,  despite  the  uncertainty  of  the 
plant  parameters.  This  means  that  for  any  plant  parameter  combination 
P^  belonging  to  the  set  P of  possible  combinations,  the  system 
output  Cj(t)  is  to  lie  within  some  prescribed  time-domain  bounds 
defining  the  region  (Fig. 1.2)  c ld>  . 


Figure  1.2.  Example  of  time-domain  specifications. 


It  should  be  noted  at  this  point  that  it  is  implicitly  assumed 
that  we  restrict  to  "reasonable"  signals  c^(t)  € 

Mathematically  this  means  that  the  designer  restricts  himself  to 


i dxc . (t) i 


signals  c.(t)  such  that 
3 


< y . (t)  for  i = l,2,...,o  for 

i 


some  given  functions  v ^ (t)  . Although  the  synthesis  technique  uses 

frequency  domain  specifications,  it  was  shown [HI]  that  frequency 
domain  specifications  suffice  for  T.D.S.  of  even  more  general  type. 

It  is  assumed  possible  to  measure  the  command  input  r(t)  and 
the  plant  output  c(t)  and  therefore  the  most  general  structure  is 
that  with  two-degrees-of-freedom,  of  which  one  is  shown  in  Fig.  1.3, 
where  F and  G must  be  chosen  by  the  designer. 

It  should  be  underlined  that  because  of  large  plant  uncertainty 


that  is  assumed  here,  it  is  impossible  in  general  to  use  an  open-loop 
system  and  therefore  feedback  is  definitely  needed.  This  necessitates 
use  of  a sensor,  which  in  turn  introduces  sensor  noise  \ (Fig.  1.3). 


Figure  1.3.  A two-degree  of  freedom  structure 

1.2  Linear  time  invariant  design. 

If  the  plant  is  minimum-phase  there  exists  [Hi]  an  infinite  set 
of  (F^GJ  which  satisfies  the  T.D.S.,  despite  plant-uncertainty 
The  restrictions  on  the  range  of  uncertainty  have  been  discussed  [H7], 
[H8] . Furthermore,  given  a region  a (Fig.  1.4)  of  width  at  most 
where  e > 0 is  arbitrarily  small,  it  is  possible  to  show  [HI]  that 
the  set  of  (F^,G^)g  which  satisfy  those  T.D.S.  is  not  empty. 


system 

output 


Figure  1.4.  Example  of  time  domain  specifications. 


If  a feedback  loop  gives  us  such  benefits  we  should  expect  to  pay 
a certain  price  for  these  benefits. 


8i 


Figure  1.5.  Typical  linear  time  invariant  design, 
both  a)  and  b) . 

Therefore  in  a general  sense,  the  less  sensitive  we  want  the 
output  system  to  be  to  plant-parameters  uncertainty,  the  bigger  will 
be  the  sensor  noise  effect  at  the  plant  input,  and  vice-versa,  and 
this  is  the  real  trade-off  in  any  feedback  system  design.  So  a reas- 
onable definition  of  optimum  is  to  select  (F,G)  € A with  minimum 
® G 2 

o„  = E / \ - — — | 0 dm  where  0 (w)  is  the  power  spectrum  of  the 

P P o nn  nn 

sensor  noise.  However,  this  is  a very  difficult  problem,  as  yet 
unsolved.  Our  definition  for  optimum  (F ^,GJ  t € Ar  is  one  for 
i ki 

which:  Min  k where  as  s -*•  ®>  , G.  (s)  -*•  — for  a-priori  chosen 

* * 1 se 

value  of  the  integer  e . Obviously,  this  definition  is  closely 


related  to  o minimization  but  it  is  definitely  not  the  same  thing 


Our  definition  (besides  giving  a solvable  problem)  is  quite  realistic 


bandwidth  minimization.  This  is 


especially  important  for  higher  order  (e.g.,  bending)  modes  avoidance 


It  is  worth  noting  that  our  definition  is  essentially  the  same  as 


Bode's  in  his  derivation  of  the  so-called  'ideal  Bode  characteristics 


[Bl]  . Using  such  an  optimum  criterion,  Horowitz  [H3]  has  proven  the 


Limitations  in  a L.T.I.  design 


The  nature  of  the  optimum  L is  such  that  Min  k 


(Fig.  1.5)  at  which 


essentially  the  minimization  of  the  frequency  w 


L(jw)  reaches  its  final  asymptotic  slope.  However,  for  minimum- 

3 £n ! T,1  / 

phase  L.T.I.  networks  — — J — L is  related  through  Bode  Integrals  to  /L 


A(w)+j  B(w)  with  A in  nepers  and  B in  radians 


f.n  w/w  then  [Bl] 


in  coth 


In  coth 


,»s  stability  requires  ^LfjcO  > -180°  for  iLfjw^)!  * 1 , then, 

because  of  large  uncertainty,  the  maximum  rate  of  decrease  of  L is 


It  will  be  shown  (section 


T.l),  however,  that  the  disturbance  specifications  are  translatable  into 


a certain  phase  margin  0..  for 


requency 


7 


range  !w  ,o>,J  (Fig.  1.6),  thus 
a b 

increasing  the  cost  of  feedback. 
Indeed,  over  that  frequency  range 


dtn  |L  | . 


do) 


is  then  at  most 

; 

M 


- ( 2 - — ) x 40  db/decade  and 
recalling  (1.1)  and  Fig.  1.5,  the 
level  of  noise  at  the  plant  in- 
put is  then  accordingly  increased. 
This  can  be  visualized  as  follows 
(Fig.  1.7) : 


Figure  1.6.  Typical  L.T.I.  design. 


Figure  1.7.  Trade-off  in  a L.T.I.  design. 


One  possible  way  to  break  this  vicious  circle  is  to  concentrate  on 

din! l| 


devices  for  which  the  Bode  relationship  relating  phase  lag  to 


do) 


V- 


■miIM 
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does  not  hold,  namely: 

NON  - (LINEAR  TIME  INVARIANT)  networks. 

I . 5 Non  - Linear  Time  Invariant  compensation. 

We  will  investigate  here: 

- NON  LINEAR  (N.L.)  networks 

- LINEAR  TIME  VARYING  (L.T.V.)  networks. 

I . 5a  N.L.  networks. 

As  was  seen,  Bode's  gain-phase  relationship  for  LTI  networks  is 
responsible  for  the  large  'cost  of  feedback'  in  the  synthesis  of 
systems  with  large  uncertainty. 

It  is  therefore  logical  to  seek  nonlinear  networks  for  which  that 
relationship  does  not  exist,  and  for  the  purpose  of  comparison,  to 
characterize  them  by  some  mathematical  tool  which  reveals  their  gain- 
phase  relations. 

The  most  commonly  used  tool  is  the  describing  functions  [G2],[M1] 

because  it  uses  the  same  frequency  domain  concepts  as  in  LTI  networks. 

The  notation  |n  ' and  / N will  be  used  to  characterize  respect- 
eq  £__eq 

ively  the  "equivalent"  gain  and  phase  of  the  nonlinear  element  N. 
Obviously,  the  equivalent  gain  and  phase  concept  is  just  a means  for 
us  to  try  to  explain,  visualize  and  invent  nonlinear  devices,  which 
are  hopefully  useful  for  our  purpose.  We  certainly  do  not  intend  to 
rely  on  Describing  Functions  to  give  a quantitative  sythesis  procedure. 
Using  the  above  notation  we  are  then  looking  for  N.L.  devj  • s such 
that: 

d'Ne  I / 

— = -m  db /decade  with  / N > - x 90® 

do:  /ea  20 


(recall  that  for  a L.T.I.  system  a rate  of  decrease  of  -m  db/decade 

is  associated  with  an  average  phase  lag  of  -m  * 90/20°)  . This 

would  then  lead  to  some  open  loop  transfer  function  L (s)  for  which: 

eg 


d|  L| 


du 


= _ 


k db/decade  with 


II  > - 

Leq 


20 


x 9Q® 


which  then  implies 


d L ! d|  L 

-SSL  < 


du 


L.T.I  I L 

— with  |Leq  = [LLTI' 


and  therefore  L (s)  would  be  as  diown  in  dashed  lines  in  Fig.  1.5. 
eq  1 

The  improvement  to  expect  in  terms  of  sensor  noise  rejection  is  then 
obvious  by  inspection  of  Fig.  1.5. 

I . 5b  Linear  time  varying  networks. 

It  is  well  known  that  time-invariant  filters  are  optimum  for 
stationary  processes,  i.e.,  if  all  signals  are  stationary  [PI] 

The  fact  that  we  concentrate  on  specific  command  inputs  starting  at 
time  t=0  , implies  non-stationary  processes.  Therefore  we  are 

induced  to  think  that  L.T.V.  filters  ( F,  G)  can  give  better  perform- 
ances than  L.T.I.  (F,G)  . Physically  it  means  that  the  feedback 

properties  are  tuned  to  follow  the  time-varying  sensitivity  needs, 
rather  than  be  time-invariant.  < 

1.6  Previous  work. 

The  above  topics  constitute  a very  broad  area  of  investigation. 
Long  ago  Holzmann (H10)  mentioned  that  "control  engineers  in  the 
chemical  process  industries  have  long  recognized  the  possibility  of 
achieving  superior  responses  by  means  of  nonlinear  control".  And 
Chesnut  [C2]  quoted  John  Xocre  "that  for  any  linear  control,  he  could 


10 


always  build  a nonlinear  control  which  would  be  better".  Clegg  [Cl] 
and  others  have  noted  that  N.L.  compensation  can  give  better  gain- 
phase  relations  than  LTI  networks  (from  the  describing  function  view- 
{>oint).  But  no  quantitative  synthesis  techniques  exploiting  this  idea 
were  reported.  Recently,  however,  Krishnan  and  Horowitz  [h4]  proposed  a 
quantitative  design  procedure  using  a N.L.  element. 

On  the  other  hand,  it  was  shown  that  "variable  structure  systems" 
have  adaptive  properties  (see  for  example,  [El],  [Ul] , [U2] , [Kl] ) . 

In  these  systems  the  LTI  compensation  changes  when  system  state  crosses 
a switching  surface.  However,  no  quantitative  synthesis  technique  for 
satisfying  assigned  specifications  over  a given  range  of  uncertainty  has 
appeared  for  this  class.  A search  of  the  literature  reveals  that 
besides  [h4]  no  quantitative  synthesis  procedure  using  N.L.  compen- 
sa*-  ion  has  ever  been  given  for  uncertainty  plants  and  which  is 

quantitatively  related  to  the  performance  spe  cifications  and  to  the 
range  of  uncertainty. 

Much  work  has  been  done  on  the  filter  problem  in  bo'h  stationary 
and  non-stationary  processes  (see,  for  example,  [Wl],  [B2],  [SI],  [K2;l. 
Very  little  has  been  done  in  this  area  for  uncertain  plants.  Fleisher 
[FI]  did  consider  the  problem  of  stochastic  plant-parameters  and  gave 
a synthesis  procedure  that  uses  lti  compensation.  However,  because 
of  the  linearization  that  he  made  in  order  to  obtain  a solvable 
problem,  his  synthesis  technique  is  approximately  valid,  at  best,  for 
relatively  small  arameter  uncertainty. 

In  summary  v»  can  say  that  the  field  of  quantitative  synthesis  of 
• icv  svstem.  ainn  N.L.  or  L.T.V.  networks  has  hardly  been  inves- 

* ■'ted  an  of  row. 
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CHAPTER  II . BASIC  PRELIMINARIES 

The  work  presented  in  this  dissertation  rests  on  the  following 
foundations : 

a)  The  synthesis  procedure  for  L.T.I.  systems  to  satisfy  T.D.S. 
despite  large  plant  ignorance,  devised  by  Horowitz  and  Sidi  [H2]  . 

b)  The  derivation  of  optimum  linear  time  varying  filters  for 
non-stationary  processes  due  to  Booton  [B2] . 

c)  The  concept  of  set-Equivalent- Plant-Linear  Time  Invariant 
(E.P. L.T.I.)  used  to  characterize  a large  class  of  non-linear-time- 
varying  plants,  introduced  by  Horowitz  [HI]. 

We  shall  now  briefly  review  those  topics: 

II. 1 Synthesis  procedure  for  L.T.I.  systems. 

The  method  presented  in  [H2]  is  restricted  to  minimum  phase 
plants,  but  has  recently  been  extended  to  non-minimum  phase  plants 
(Horowitz  and  Sidi,  [H9] ) . in  a first  step,  the  T.D.S.  are 
translated  into  frequency  domain  specifications  (F.D.S.).  Although 
no  general  rigorous  translation  between  the  two  domains  is  known, 
there  is,  in  practice,  little  difficulty  in  effecting  one  which  is 
satisfactory  for  practical  engineering  purposes.  Thus,  for  example 
[S2],  the  T.D.S.  of  Fig.  1.2  are  translated  in  the  F.D.S.  of  Fig.  II. 1 
and  therefore  at  each  frequency  w € f!  = ( 0,®]  , the  maximum 

tolerable  variation  A ln[T(iw)j  on  the  system  overall  transfer 

Q 

function  T £ — (see  Fig.  1.3)  is  defined.  The  synthesis  technique 
is  mathematically  rigorous  for  such  w-domain  specifications. 

Secondly,  at  each  u>?  € H , P ^ ( ju>£ ) denotes  a point  in  the 

- 


is  the  sensitivity  of  T(jw)  to  uncertainty  in  P . This  is  some- 
times denoted  as  the  'equality  of  positive  and  negative  feedback  areas'. 
Hence,  if  1 S | < 1 over  some  range  (the  primary  objective  ii  a feedback 
system),  it  must  be  balanced  by  another  in  which  !s|  > 1 . It  is 

easy  to  live  with  this  constraint  in  the  two-degree-of-freedom  system, 
because  at  sufficiently  large  ui  , ! T is  negligibly  small,  so  that 
large  relative  changes  n ' T ' are  inconsequential.  This  means  that 
at  large  u , it  is  essential  that  the  response  tolerances  A ln' T(jui)i 
exceed  the  range  of  the  P(j u>)  template  in  order  that  L(joi)  may  go 
to  zero  as  in-*-®  . 's'  >>1  is  tolerable  in  this  high-frequency  range, 
as  far  as  its  effect  on  T(ju)  is  concerned,  because  the  prefilter  F 
in  Fig.  1.3  attenuates  the  resulting  high  peaking  in  L/(l-*-L) 

However,  the  disturbance  response  T.  « C/D  = (1+L'  1 = S of  (2.3), 

G 

is  then  also  very  large,  which  is  generally  not  tolerable  because  there 
is  no  equivalent  filter  available.  Although  the  parameter  ignorance 
problem  is  assumed  to  dominate,  it  is  necessary  to  consider  the 
disturbance  response,  at  least  to  the  extent  of  adding  the  constraint 

3y>  0,  3 [ Td|  4 i||  = | (1+L)‘X|  £ y , Vu  (2.4) 

The  value  of  y may  be  related  to  the  damping  factor  C of  a second- 
order  response  function  by  the  relation 


, 1 , _ , 1 
Y=  max  | 1 = ■ -■  -v. 

u>  ; l-w  + j 2C<u  1 ! 2c*/’(l-c  ) 1 

The  above  leads  to  the  boundary  of  Fig.  II. 2 with  0 a function 

* m 

of  y , e.g.  6 = 50°  if  y = 2*3  dB  . At  low  frequencies  the 

rn 

parameter  factor  dominates  for  w < in  Fig.  II.2,>  so  these  boundaries 
containno  cart  of  . There  is  an  intermediate  frequency  region  in 


% 


which  part  of  the  boundary  contains  a portion  of  r„  , e.g. 
in  Fig.  II. 2.  At  sufficiently  large  u>  becomes  the 

complete  boundary. 

In  a fourth  step,  using  the  optimum  criteria  of  1.3,  i.e.. 


Min  k,,,  , where  as  s-*"  , G.  (s)  -*■  for  some  given  integer  e 

i e 

s 


the  optimum  L^tjw)  was  proven  to  be  the  one  which  lies  exactly  on  its 


bound  at  each  € ft  . (see  [H2 ] , [H3]  and  [H9]). 


G(jw)  only  guarantees  that  A In!  T(jui)  | is  satisfied.  The  pre 


filter  F is  needed  in  order  for 


i PL 


, . to  lie  in  between  IT... 

1+L  Min 


and  ' T 

Max 


F is  thus  derived,  for  example,  by  matching  ' F ' 
the  specified 


|-£-l 

1+L 


with 


Max 


T (jio)  ! for  all  u>  € Q . 
Max 


Figure  II. 2.  Typical  bounds  on 


Ln(jw)  at  different  frequencies . 


II. 2 Design  example. 


The  plant  is  P(s)  = — with  k € [1,1000] 

s +2c  a;  s+w 
P P p 


Cu  € [-5,3] 
P P 


u /.  2 € [2,10] 

p i-cp 


the  parameter  being  independant.  It  is  required  that  the  system  step 
response  lie  within  given  T.D.  bounds  over  the  entire  range  of  plant 
parameter  values.  This  problem  was  solved  by  Horowitz  and  Sidi  [H2] 
and  [S2]  and  the  resulting  F.D.S.  are  shown  on  Fig.  II. 1.  A peaking 
Y = 2 dB  was  tolerated  which  corresponds  (Appendix  Al)  to  C = .44  and 
to  a tolerable  overshoot  of  21%  in  the  step-disturbance  response. 

Using  the  procedure  described  above, 

. 47.191015(s+.905)  (s+,92)  (s+17.3)  (s+26.1)  (s+200.Ws+2220) 

hi'-s>  - o : 

s (s-t.512)  (s+5.28)  (s+6.1)  (s+37.8)  (s+50)  (s+1000)  (s  +10500s+(15000)  ) 


associated  with  the  nominal  plant  P (s)  = 5 

nom  - 

s(s  +2s+5) 

Lj(s)  is  plotted  in  both  Figs.  III. 13  and  III. 14.  The  prefilter 


is  then  calculated  to  be  F(s)  = 


9 (s+,94) (s+2 .4) 
(s+1.5)2(s2+4.2s+9) 


C or responding 


command  outputs  and  step-disturbance  responses  are  plotted  in  Fig.  III.  18 
and  Fig.  III. 19  respectively  for  different  plant  parameter  values. 

II. 3 Optimum  L.T.V.  filter  for  non-stationary  inputs. 

Given  a set  of  input  signals  that  are  to  be  applied  to 

a system  characterized  by  its  impulse  response  h(t,t)  , i.e. 

(Fig.  II. 3 ) : 


x„(t)  = f h (t,C)  x_(0  d? 
R 0 1 


Furthermore  let  us  assume  that  any  input  signal  Xj(t)  £ {x^}  is  the 


16 


17 


of  the  integral  equation: 

t 

Vl0<t,T)  = / h(t,S)  Y U,T>  dc  for  t > T (2.5) 

0 

Given  the  deterministic  signal  s(t)  , how  can  one  approximate 

the  autocorrelaticn  function  Yss(t^,tj)  for  practical  purposes?  We 

first  note  that  Y^ft^.t^)  a symmetric  function,  continuous  if 

Xs  is  continuous.  Let  us  consider  signals  existing  on  [0,T]  where 

T can  be  as  large  as  desired.  We  can  consider  y (t, ,t_)  as  a 

ss  1 2 

definite  symmetric  kernel  and  by  application  of  Mercer's  theorem  [C3] , 
[C41 


» <P  (t_)<P  (t.) 
r <t  t > = l -r  -1-  r 2 
IX  1 2 r=l  *r 


(2.6)  for  all  t ,t  € [0,T] 


where  X^  and  ( • ) are  respectively  eigenvalues  and  corresponding 
eigenfunctions  of  y , i.e.,  non  trivial  solutions  of: 

1 T 

— «>r(s)  =/  yi]:(s,C)  <P  (C)  dC  (2.7) 

r o r 

For  practical  purposes,  we  can  certainly  approximate  (2.6)  by: 


YII(tl't2) 


N <p  (t.)4>  (t  ) 

i — V— 


r=l 


(2.8)  for  all  tj.tj  € [0,TJ 


where  N is  such  that  the  norm  on  the  space  L (0,T)  of  the  error 
is  less  than  a given  tolerated  error  e , i.e.. 


TT 
//  ! Y 


(t  )ip  (t,)  I » . 2 

„(t  ,t  ) - l j-Z i dtldt  - l I—  | < e (2.9) 

00  1 r 1 2 r=N+l  >X 

r 

II. 4 Concept  of  set-Equivalent-plant-linear  time  invariant. 

It  was  suggested  long  ago  to  replace  non-linear  or  linear  time 
varying  networks  by  an  approximate  equivalent  linear  time  invariant 
ones, for  convenience  in  analysis  and  synthesis.  The  Describing  function 


(G2 1 is  a well  known  example  of  such  a substitution.  However,  it  was 
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definitely  not  used  as  a tool  for  quantitative  synthesis  in  uncertain 
systems.  The  concept  of  Equivalent  riant  Linear  Time  Invariant 


(E.P.L.T.I.)  set  which  can  be  properly  used  for  quantitative  synthesis 

was  recently  introduced  by  Horowitz  [HI,  H6] . This  concept  is  now 

briefly  reviewed.  Given  fig.  II. 4)  R={r}  and  D={d)  - two  finite 

sets  of  deterministic  command  and  disturbance  inputs. 

Let  I = {ia(t) ,a=l, . . .N}  be  the  given  set  of  system  inputs  for  which 

a 

the  design  is  to  be  executed,  where  i may  consist  of  a command,  a 


disturbance,  or  a combination  of  any  two  of  these.  Given  the  set 
Ca[t]  =fca(t)}  of  acceptable  system  output  signals  associated  with 
ia(t)  , it  is  possible  to  find  the  corresponding  set  za[t]  = (za(t)} 
of  acceptable  plant  output  signals  and  thus  the  set  7a[s]  = {#za(t)>}. 
Let  the  plant  be  characterized  by  a set  W = {w}  (because  of  plant 
uncertainty)  of  nonlinear  time  invariant  (the  latter  for  the 
sake  of  simplicity  - the  extension  to  nonlinear  time  varying  is 
straightforward  ) continuous  mappings  w : X(t)  -*z(t)  , with  unique 
continuous  inverse  w 1 (31) . Then  for  any  plant  w^  € W there 
exists  a plant  input  signal  X*  (t)  which  produces  the  plant  output 
signal  z^(t)  € Za[t]  and  we  can  define: 


P.  (s)  = 

IV 


_ zv(s> 

"«£&Si(t>l 


as  the  a.  - equivalent  plant  linear  time  invariant  (E.P.L.T.I.) 

transfer  function,  in  the  sense  that  the  input  xa.(t)  into  the  linear 

vi 

n i 

Piv  (s)-  gives  the  same  output  as  the  input  x^(t)  intq  the  non-linear 


plant  w. 


Thus  the  method  is  restricted  to  nonlinear  w for  which 


the  set  X is  Laplace  transformable.  Restriction  is  made  in  the  mean- 


time to  design  procedure  for  w such  that  P (s)  is  minimum  phase 


(p  (t)}  = R 


{»'))} 


0={<J(t)} 


Figure  II. 4.  Feedback  structure  with  nonlinear  uncertain 
plant  set  W 

(B2 ) for  all  a,  i,  v 

As  i ranges  over  all  w^  € w and  v ranges  over  all 
z*(t)  € Z0[t]  , P.^(s)  describes  the  set: 

Pa[s]  = {P.“(s)  , z“(t)  € Za[t]  , w.  e W} 

which  is  defined  as  the  a -E.P.L.T.I.  transfer  function  set  associated 
with  the  input  signal  ia(t)  over  the  output  set  Ca[t]  and  the  non- 
linear plant  set  W 

Once  this  is  done,  the  quantitative  design  procedure  used  for 
L.T.I.  systems  [H2]  can  be  applied  to  the  E.P.L.T.I.  problem  where  the 
specifications  on  the  acceptable  output  are  those  assigned  to  the  non- 
linear problem.  Pa[s]  becomes  the  uncertain  L.T.I.  plant  set  in  this 
equivalent  L.T.I.  problem.  Schauder's  fixed  point  theorem  is  used  to 
prove  that  the  solution  to  the  latter  problem  is  valid  for  the  nonlinear 
plant  set.  The  synthesis  procedure  of  Paragraph  II. 1 can  be  easily 
extended  to  nonlinear  time  varving  plants. 


I 


f " 1 

* 
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TN  SYSTEMS  WITH  LARGE  PLANT  IGNORANCE . 

III.l  INTRODUCTION . 

The  Clegg  Integrator  (C.I.)  is  a linear  integrator  (L.I.)  whose 

output  y is  reset  to  zero  whenever  the  input  x crosses  zero  ( [Cl] , 

see  Fig.  III.l).  If  it  is  properly  inserted  in  a feedback  loop  system, 

it  can  reduce  considerably  the  system  step  response  overshoot. 

The  C.I.  response  y(t>  to  a sinusoidal  excitation  x(t)  is 

shown  on  Fig.  III. 2 and  its  describing  function  N for  sinusoidal 

inputs  is  plotted  on  Fig.  III. 3.  It  is  seen  that  the  rate  of  decrease 

of  | N j and  of  the  L.I.  are  the  same.  However,  j N = -38°  compared 

to  -90°  for  the  L.I.  Therefore  the  phase  margin  9 associated 

M 

with  the  "non-linear"  open  loop  transfer  function  L^ts)  containing  a 

C.I.  is  bigger  (by  52°)  than  the  one  associated  with  L (s)  contain- 

LTI 

ing  a L.I.  in  place  of  the  C.I.  As  we  know  phase  margin  is  closely 
related  to  the  overshoot  in  the  system  step  response.  This  explains 
qualitatively  the  differences  mentioned  earlier  noted  by  Clegg.  However, 
quantitative  synthesis  can  not  rely  on  Describing  Function  Theory  (as 
already  stated)  and  in  any  realistic  design  one  has  to  overcome  a big 
instability  problem  if  one  wants  to  use  a C.I.  (see  Chapter  IV).  To  a 
large  extent  the  'First  order  reset  element'  (FORE)  (which  is  the 
linear  element  1/s+b  whose  output  y is  reset  to  zero  whenever  the 
input  x crosses  zero  (Fig.  III. 4))  overcomes  the  instability 
problem  mentioned  before,  and  this  will  be  explained  in  detail  in 
Chapter  IV. 

It  is  worth  noting  (Figs. III. 3,  III. 2, 


III. 5),  that  FORE  appear'- 


Figure  III.l.  C.I. 


Figure  III.4.  FORE. 


(COM  6* 


w) 


I 

0 
, i 

-2 


Normalized 
reset  value 
(%  of  the  Cl 
reset  value) 


.01 


.98 


.1 


.86 


] 


.26 


.10 


.05 


.01 


Figure  I I 1.2.  C.I.  response  to  Figure  III. 5.  FORE  response  to  a sine 
a sine  wave.  wave. 


Figure  III. 3.  Normalized  describing  functions  of  C.I.  & FORE. 

as  a linear  element  to  slow-varying  signal  and  as  a non-linear  element 
to  fast-varying  signals, because  the  value  of  the  reset  (Fig.  III. 5)  is 
definitely  a function  of  the  input  signal  frequency.  This  gives  a flex- 
ibility, not  available  in  the  C.I. 

Furthermore  the  frontier  between  the  two  characters  can  bo 
adjusted  to  a particular  design. 

III. 2 Analysis:  Characterization  of  'FORE1. 


Some  tools  are  first  derived  that  are  used  to  solve  the  synthesis 


problem. 


_ «... , , 


-*-r 


III. 2a.  Equivalent  linear  representation. 


Let  9 = {t^  , i=l,2...  } be  the  set  of  zero-crossings  of  the 
input  signal  x(t)  to  a FORE  (Fig.  III. 4),  i.e.  t.  € 9 is  such  chat: 
x(t.)  = yf(t.+)  = 0 

Then  for  t„ < t < t„ , , , where  t„  , t„ , , € 9 , we  have : 

N N+l  N N+l 


yf(t)  = ; X(ae-b(K)  dC 


If  u(t-tK) 


f 1 for  tit, 
0 for  t < t, 


K denotes  the  Heaviside  function 


r _b(t .c)  i 

then:  yf(t!  = [u(t-tQ) -utt-t^)  ] f / X(C)e  dsj  + 

t2 

!u(t-t1)-u(t-t2) ] f X(?)e'b(t‘C)dC  + ..  + 

fcK+l  ..  , 

(u(t-tK)  -u(t-tK+1)J  / x(C)e  ^ ^;dc  + ..  + 

t bK 

u ( t— t ) f x(()e'b(W!d c = u(t-t  ) f i(Oe'b(t'Ud?  - 

tn 

t N 0 


[ I / X(Oe' 
K=1 u t„  , 


(3.0a) 


Let  x * & f X(Oe  b^K  ^d£  (3.0)  and  notice  that 

' Vi 

y^(t)  = / x(Oe  ^d£  is  the  output  signal  of  the  linear  element 

b0 

1/s+b  for  the  same  input  x(t) 


n rt 

Then  y,(t)  * y (t)  - ] ■ f e T 5(x-t  )dr  X * u(t-t  ) (3.1) 

' K=1  Lt  K * 

0 

where  S (t)  is  the  unit  impulse  function. 

At  this  stage  several  comments  are  necessary. 


i 
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a) . Not  I*  that  x * represents  the  reset  value  that  occurs  at  time  t=t  . 

K K 

Thus  the  nonlinear  aspect  of  FORE  is  represented  (3.1)  by  a train  of 

* 

impulses  of  value  -XR  at  tR  € 9 , added  to  the  input,  i.e.  FORE  can 

be  replaced  (Fig.  III. 6)  by  1/s+b  and  adding  to  the  input  the  train  of 
N 

impulses  - Y x*  6(t-t„)  for  t < t < t„._  . B)  . We  also  learn  from 

^ K K N N+l 

(3.1)  that  if  x(t)  produces  yf(t)  , then  k-xtt)  (where  k is  some 
real  constant)  gives  the  output  k*y^(t)  and  therefore  the  linear  chara- 
cter is  preserved  in  the  multiplication  by  a constant.  This  is  quite  an 
important  feature  because  it  is  then  sufficient  to  characterize  FORE  for 
unit  signal,  command  or  disturbance,  in  order  to  know  its  behaviour  for 


the  non  unit  one.  Y) . From  (3.1)  one  notes  that  FORE  is  nonlinear 

from  the  additive  point  of  view.  If  X^(t)  , y^(t)  , and  x^(t)  , 

2 

yf(t)  are  respectively  paired,  then  in  general,  X3<t)  = x^(t)  + X2(t) 

3 1 2 

gives  rise  to  an  output  signal  yf(t)  ^yf(t)  + yf(t)  . The  equality 
3 12 

yf (t)  = yf(t)  +yf(t)  occurs  if  and  only  if  the  sets  of  reset  instants 
6^  and  02  associated  with  X^(t)  and  X^ft)  are  equal.  Thus  FORE 


is  a linear  element  over  the  class  X of  input-signals  X(t)  that  have 
the  same  zero-crossings.  However,  even  when  FORE  acts  on  X,  it  has  not 
the  commutativity  property  with  linear  elements,  i.e.  FORE  * LTI f LTI  * FORE . 
All  these,  then,  legitimize  as  a notation,  the  use  of  (^g)  * to  chara- 
cterize FORE.  From  now  on,  any  FORE  will  be  followed  (even  if  not 
specifically  mentioned)  by  the  element  's+b)  (Fig.  III. 7)  and  we  define 


/-La* 

y — x m 

i 

Vs+bJ 

s+b 

K-lV(t-V  fcK€ 

Figure  III .6.  Equivalent  representation  of  FORE. 


acre?* 
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the  equivalent  open  lrv-n  transfer  function 


L (s)  ^G(s)  P (s)  . Note  that 

eq 


L (s)  w'-if'n  no  reset  occurs.  In  a similar  manner  the  eauivalent 

oq  ETI  ^ 


overall  transfer  function  is  defined  as  T (s)  =F(s)  L !s)/(l+L  (s)) 

eq  eq  eq 


By  using  the  equivalent  representation  of  FORE  shown  on  Fig.  III. 6 and 


by  application  of  the  superposition  theorem.  Fig.  III. 7 becomes  equiva 


lent  to  Fig.  III. 8 from  the  output  signals  point  of  view.  Thus  the  non1 


inear  output  C (t)  is  the  combination  of  the  linear  system  output 


For  example,  in  Fig.  III. 7,  let  R(s)  * 1/s  , D = N=0  , F(s)  = 1 


with  to 


first  reset  instant  t is  such  that  C„(t,)  =1  (because  *(t,)  = 0) 


the  nonlinear  system  output  C (t)  = C (t)  ■‘■C  (tl  = 1 


1/1+L 


Figure  I I I. 8.  Equivalent  system  block  diagram 


ir  m 

T«’F  w • L*PG 
-t 

M 0 
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In  summary: 

f CP(t) 

for 

t ( 

W"  - 

1 1 

for 

t > 

This  result  is  very  interesting  because  the  steady  state  is  reached 
in  a finite  time  t , which  is  impossible  to  acheive  with  L.T.I. 
systems  and  because  the  result  (0%  overshoot  in  the  system  step  response) 
is  independant  of  £ and  . However,  it  is  of  little  practical 

use  because  of  the  pure  2nd  order  system  considered  here. 

Physically  the  above  is  easily  explained.  After  the  first  re- 
set occured,  the  steady  state  (*=0  , y « 0 , c * 1)  is  reached  which 
eliminates  the  dynamics  for  t > t^  , despite  possible  ignorance  in 
u 

N 

III. 2b.  Non-linear  step  response  overshoot. 

Using  some  of  the  results  of  the  previous  section,  it  is  poss- 
ible to  calculate  the  nonlinear  overshoot  that  occurs  in  the  step 
disturbance  response  of  the  system  shown  in  Fig.  III. 7.  For  practical 
purposes  because  of  the  universal  character  of  the  'optimum'  L.T.I. 
open  loop  transfer  function  L(jo>)  in  problems  with  large  parameter 

uncertainty,  it  is  possible  to  approximate  the  resulting  high  order 

V e'tds 

L/(l+L)  by  a delayed  second  order  function  T»pp(s)  “ — r . 

s%2Co)ns-hon<: 


Figure  III. 9.  Linear  and  nonlinear  system  step  responses. 


26 


(3.1b) 
For  the 


This  approximation  is  justified  in  Appendix  Al.  In  the  second-order 

-cv/i-s2 

system  the  overshoot  is  related  to  £ by  = e 

and  to  the  magnitude  peaking  by  k-Tftjw)  ! = — — 

'1+L  1 max  2i;/jR2 

remaining  of  this  paragraph,  t^  is  assumed  to  be  zero.  Therefore  the 
following  results  will  only  be  applied  to  those  systems  for  which  the 
time  delay  can  be  neglected  (see  Appendix  1) . 


Consider  the  system  of  Fig.  III. 7,  where  r = n = 0 and  d(t)  is 


a step.  Then  in  Fig.  III. 8 : C.  (t)  = — sin  (tu^t+cos  C)  with 


/l^2.  w„ 


0 - ~N  _x 

IT-  cos  K 

The  first  reset  instant  occurs  at  t,  = and  therefore 

1 


the 


first  reset  value  is:  ( = -C^) 

|*  = - /c„( C)e'b(tl"°di;  - + i e~btl  + ^-r-i  / ebC  C 


\ ■ ■ ; ''A""'’  ~ - ■ b " - ’ b 

0 -Cu„t 

V N 

integrating  by  parts.  As  C5  (t)  = — ■ — — sina)Qt 


C^OdC  after 


we  have 


Y = + 


-bti  u 
e 1 N 


1-2?^  + (!f-)2 

b b 


✓1=5- 

^ e^N*!  + — e'btl  ^ (3.2) 


The  nonlinear  system  step  disturbance  response  is  therefore: 

for  t 6 t. 


CNL(t) 


CJl(t) 

C4(t) -X1*Cjl(t-t1)  for  tx^t<t2 

with  t^  , tj  £ 0 . 


(3.3) 


There  exists  at  least  one  extremum  at  tQ  E [t^t^l  at  which, 

Sl(V  = W - y Ct(t0-tl)  = 0 (3.3a) 

Replacing  each  term  by  its  respective  value  gives: 

!,inVo  • V <^2  c°sV VV  - CsinWY,e  » *j 


>inw„  (t„-t  ) ) 1 ] 


Using  sinw. t,  * /1-C  and  cosu>  t,  =■  we  have: 
Cl  o 1 


1 
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- Cut 
we  NO 

— r (l-X*eWN  Iid  I sinw  t = 0 

■2  1 N 0 0 


(3.4a) 


'1-? 
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As  t-  t [t.,t  1 , there  is  only  one  extremum  at  t.  = — (3.4b)  which 

U 1 ^ 0 Wq 

corresponds  to  the  instant  at  which  (t)  has  its  first  minimum. 

The  value  of  the  extremum  is: 

°V  * 'CNL(t0)|=  IVV  - ^'VVl'1  = 0VUN  - A (3‘5) 


with  A 


1-2CM  +M 2 
^ 1 1 


(M^e-0*  - M1(1-2CM1)  e'U/Ml) 


R = e 


-Cccs 


(3.7) 


M = /2cZ+i/4C4  + 1 

1 b (3.8) 


(3.6) 


u _ ff-cos  C 


(3.9) 


and  u denoting  the  crossover  frequency  of  L (s)  , i.e.,  u 

c eq  c 

such  that  ! L (ju  ) ; = 1 
eq  J c 1 

Using  (3.5)  and  (3.6),  charts  are  presented  in  Fig.  III. 10, 

giving  the  net  overshoot  or  undershoot  in  the  step  disturbance 

u 

Q 

response,  for  the  non-linear  design,  with  M « — as  a parameter. 

D 

It  is  reasonable  to  hypothesize  that  for  values  of  t very 

d 

small  relative  to  t^  and  tQ  , Fig.  III. 10  can  be  used  with  good 
accuracy.  Experimental  results  support  this  and  will  be  given  later. 

III.2.C  When  is  it  possible  to  consider  FORE  as  a linear  element? 
Assume  FORE  is  located  as  Siown  in  Fig.  III. 7.  and  let  d(t)  be 
a step-disturbance  and  r = 0.  Let  C^tt)  (Fig.  III.8J  be  as  shown  in 
in  Fig.  III. 11,  then  for  t < t^  we  have: 

tl 

x(t)  = (t)  = -Ce(t).  Also  xi*  - - / Ct(t)e-b(t"tl)dt  . (3.10) 

If  b = 0 (C.I.  case),  x^»  is  the  area  under  from  t«0  to  t^ 

in  Fig.  III.  11.  With  b^O  , X^«  is  the  area  under  curve 

_ _ (J_.  -b(ti-t)  . If  bt,  >1  bv  a factor  of  3 or  more,  B is  very 
B = C,  (t)  e 1 1 1 1 
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satisfied  it  implies  that  'FORE'  can  be  considered  as  behaving  like 
1/s+b  from  a practical  point  of  view. 

III. 3.  Synthesis. 

III. 3a.  Philosophy  of  the  design  procedure. 

We  use  the  two  degree  of  freedom  structure  of  Fig.  III. 7.  The 
open  loop  transfer  function  L (s)  has  to  satisfy  bounds  [H2]  at 

LTI 

each  frequency  oj  € (Fig.  III. 12).  In  general  in  the  low- 

frequency  range  (w  < w ) L (jui)  lies  on  bounds  due  to  the  specifi- 

D LTI 

cations  on  the  command  input  responses  and  in  w > io  L (s)  is 

o LTI 

determined  by  the  maximum  tolerable  overshoot  in  the  step  disturbance 


response.  This  merely  means  that  for  u>2  > u>  > u>b  , the  bound 

imposes  a certain  phase  margin  9 for  L ( jto) 

M LTI 
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Recall  that  the  smaller  the  phase  margin,  the  bigger  is  the 

tendency  to  instability  characterized  by  a larger  overshoot.  The  N.L. 

character  of  the  FORE  in  this  higher  frequency  range  (tu  ' to  } , will 

b 

permit  a smaller  phase  margin  6'  on  L (s)  for  acheivinc  the  sane 
specified  maximum  tolerable  overshoot  in  the  step  disturbance  response. 

This  allows  a faster  decrease  of  ' I,  ( j<*>) I (Fig.  III. 12)  and  is  accom- 
panied by  the  reduction  of  the  sensor  noise  effect  at  the  plant  input. 

We  can  adjust  the  value  "b"  of  FORE  such  that  in  the  lower  frequency 

ranae  (u < w.  ) FORE  acts  like  a linear  element,  so  that  the  bounds 
b 

on  L (s)  due  to  the  command  input  are  unchanged.  Consequently  the 
prefilter  F remains  L.T.I.  which  is  certainly  an  interesting  feature 
and  one  of  the  important  advantages  of  FORE  over  the  C.I.  Therefore  a new 
set  of  bounds  on  L (s)  has  only  to  be  derived  in  the  higher  frequency 
range  (w > w.  ) . This  is  done  by  using  Chart  III. 10  as  will  be 

D 

explained  in  Part  I.  For  high  plant  gains,  an  undershoot  problem  arises 
which  is  treated  later  in  Part  II. 

III.3.b.  Nonlinear  Design  Details  - Part  I. 

Without  loss  of  generality,  the  design  details  are  better  presented 

by  means  of  an  example.  We  will  take  the  example  of  section  II. 2 

treated  in  [H2!  . The  bounds  T (u)  due  to  command  response  specifications 

and  the  bound  , due  to  the  disturbance  response  specifications, 

are  plotted  in  Fig.  III. 13.  To  simplify  the  design  with  little  loss  in 

exploitation  of  the  benefits  of  nonlinearization,  the  parameter  b of 

FORE  is  chosen  approximately  at  the  first  w value  (denoted  by  w ) 

b 

at  which  E?(jw)  lies  on  r.  in  Fig.  III. 13,  i.e.  at  which  the 
disturbance  response  dominates,  in  the  present  case  at  60  rps  (Fig.ITI.13.) . 
Tills  definitely  assures  that  the  (nonlinear  response  to  r ) = (linear 
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« 


Figure  III. 13.  Linear  and  nonlinear  designs  on  Nichols  chart  (at  K=l) 
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i I 


L 


response  to  r ) 


From  Fig.  II. 1.  the  maximum  bandwidth  of 


T(jto)  » 7 rps.  From  the  relation  Bw  = u [1-2^  2+ (2-4j;  2+4r  4) 2 ) ^ , <x 

n n 


is  certainly  less  thail4,  so  (3.11)  is  easily  satisfied.  The  above 
choice  for  b considerably  simplifies  the  design  details,  because 
only  the  disturbance  response  need  be  considered,  meaning  that  only 


r„  is  changed  into  r „ ; the  other  r (w . ) boundaries  need  not  be 

l ML  1 


recalculated,  in  general,  in  problems  with  large  parameter  uncertainty, 
this  method  of  choosing  b will  have  (3.11)  easily  satisfied  which  then 
assures  that  the  (nonlinear  response  to  r ) = (linear  response  to  r ). 
It  is  for  such  problems  that  the  added  complexity  and  resulting 
bandwidth  saving  of  a nonlinear  design  are  justified. 


How  is  r . obtained?  Of  the  total  ensemble  of  plants,  consider 
n SL 


that  subset  whose  loop  transmission  L=PG  has  crossover  frequency 

Since  there  may  be  an  infinite  number  of  plants  in  this  sub- 

From 


to  =60 
c 


set,  pick  the  one  with  the  largest  value  of  | L/(l+L) j 


Fig.  III. 10  at  M=1  , 20%  overshoot  gives  8 * 180-155  - 25°  . We 


position  the  plant  template  at  to=60  (shown  shaded  in  Fig.  III. 13), 
such  that  one  extreme  point  (X)  at  which  it  cuts  the  zero  db  line 
is  at  9m  = 25°  . Another  extreme  position  is  at  Y at  which  6 =-138°, 
for  which  the  overshoot  is  about  3%,  from  Fig.  III. 10.  Because  of  the 
shape  of  the  template,  it  is  seen  in  Fig.  III. 13  that  the  corresponding 
extreme  position  of  L (nominal)  is  on  the  vertical  line  9 =-144°  , 


i.e.,  in  the  nonlinear  desion  at  <o=60,  L , ,,  may  lie  on  or  to 

(nominal)  1 


the  right  of  r (60) 
n?. 


, the  latter  being  the  vertical  line  9 =-144°. 


The  procedure  is  illustrated  again,  at  to  = 100.  Consider  that  plant 


value  whose  loop  transmission  has  u = 100  and  gives  the  largest 


value  of  L/(1*L)  ' . 

shoot  oives  9 


'100  ’ From  Fiq'  at  M= 100/60,  20%  over- 

180-165.5  = 14.5°  . We  position  the  u = 100  plant 


template  such  that  the  extreme  point  at  which  it  cuts  the  zero  db  line 

is  at  0 =14.5  . The  corresponding  positions  of  L (nominal)  is  on 
m 

the  vertical  line  9=-160°  , which  becomes  r (100) 

nx. 

Some  of  the  experimental  results  are  given  in  Figs.  III. 14, 

III. 15  at  this  point  to  support  the  above.  Fig.  III. 14  presents 
Bode  plots  of  three  nonlinear  designs  whose  differences,  significant 
only  for  u>>1000  , will  be  explained  later.  Consider  the  response 

for  k = 10  in  Fig.  III. 15.  In  Figs.  III. 13,  III. 14,  the  corresponding 
a>c  is  obtained  by  raising  the  nonlinear  L ^ (ju)  by  log  10=20  db 
or  alternatively  by  lowering  the  zero  db  line  by  -20db  , which  cuts 
very  close  to  60  rps  and  -140°  . Hence  M=1  and  from 
Fig.  III. 10  at  140°  on  the  M=1  curve,  the  net  overshoot  predicted 
is  5%,  almost  precisely  that  obtained.  Similarly,  consider  k=25  in 
Figs.  III. 13,  III. 14,  III. 15.  The  corresponding  wc  is  obtained  by 
examining  the  -20  log  25=-28  db  line  which  cuts  L^fju)  close  to 
(o=1000  at  -160°  . In  Fig.  III.  10,  on  the  curve  M=100/60 

(extrapolated  from  M=l,2)  predicts  10%  overshoot  compared  to  9.1% 
obtained.  At  k = 100  , on  the  corresponding  -40  db  line:  = 200  , 

0=-172°  and  in  Fig.  III. 10  the  corresponding  extrapolated  M=200/60 
curve  predicts  about  12%  overshoot,  compared  to  11.4%  obtained  in 
Fig.  III. 15. 

III.3.C.  Nonlinear  design  - Part  II. 

Suppose  the  above  approach  based  on  a second-order  model,  is 
pursued  over  the  entire  u)  range,  and  the  resulting  nonlinear  design 
implemented,  as  was  in  fact  originally  done  in  this  research.  It  is 
found  that  the  results  are  as  predicted  for  command  inputs  r(t)  for 
all  plant  conditions,  and  for  disturbance  inputs  corresponding  to  the 
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smaller  plant  gain  factor  values,  from  k»l  up  to  approximately 
k = 100  . As  the  plant  gain  factor  is  further  incresed,  there  is  a 

transition  from  the  second-order  type  response,  in  the  manner  indicated 
in  Fig.  III. 15. 

■ There  are  two  important  phenomena  to  be  noted.  One  is  the  local 

maximum  marked  A in  Fig.  III. 15,  which  gradually  increases  as  k is 
increased,  until  k>100  approximately,  it  becomes  the  major  maximum, 
as  it  clearly  is  for  the  case  k»1000  . The  second  phenomenon  is  the 

local  minimum  immediately  following  A , which  eventually  (at  k=100) 
leads  to  an  undershoot,  growing  by  the  time  k=1000  , to  31%  under- 

shoot; both  phenomena  not  at  all  as  predicted  by  the  second-order  model. 

1 

j III.3.d.  Explanation  of  divergence  from  second-order  model ■ 

i 

The  second-order  model  is  satisfactory  for  all  plant  parameter 

i 

combinations  P.^  such  that  the  effective  important  part  of  L.  has 
relatively  few  poles  and  zeros,  even  though  L has  exactly  the  same 
number  of  poles  and  zeros  over  all  P parameter  combinations. 

Consider  with  say  7 zeros  and  12  poles,  such  that  most  of  the  pole 

and  zero  corner  frequencies  occur  at  'L^(jw)1  very  small.  Thus  in 

%.  Fig.  III. 14,  the  loop  transmission  at  k=25  is  obtained  by  lettina 

V 

1 I 

the  present  -28  db  line  be  the  zero  db  line.  ^-25^“^  can  be 
reasonably  approximated  by  a pure  second  order  system  for  almost  a decade 
* beyond  this  point.  The  value  of  L25<jl000)  is  approximately 

-70+28  = -44  db  and  so  is  ' L25^^1+I,25^  jlOOO  ' since  1+L.,5  ( jlOOO)  =1  . 
Therefore,  the  large  number  of  poles  and  zeros  introduced  in  the 
w=(1000, 3000)  range  are  comparitively  far-off,  and  therefore  the  time 


delay  t^  can  be  neglected  in  our  approximation. 

However,  at  larger  k values  these  corner  frequencies  become 


/l  (LINEAR) 


I L ! (LINEAR 
Jv  DESIGN) 


Transition  f rom  second  order  type  of  responses 
as  k increases. 
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less  far-off  and  have  a decided  influence  on  the  response,  because 
they  occur  in  a frequency  range  where  L/(l+L)  is  no  longer  very 
small.  As  k increases,  the  corresponding  L/(l+L)  is  a high-order 
system,  and  t^  can  no  longer  be  neglected.  The  response  of  a time- 
delayed  second  order  system  is  compared  in  Fig.  III. 16a  with  the  actual 

L 

linear  step  response  of  — — - . Note  how  well  the  two  agree, 

1+L1000 

despite  the  fact  that  is  a very  high-order  system.  The 

nonlinear  step  response  is  obtained  from  the  linear  step  response,  by 
mans  of  Eq.(3.1).  The  first  reset  is  at  t^  , at  which  the  impulse 
response  due  to  an  input  -x^*5 (t-t^)  must  be  added  to  find  the  effect 
of  FORE.  Because  of  the  time-delay  t^  , the  effect  of  the  impulse 
is  n.  i felt  until  t > t^ + t^  , thus  explaining  the  local  maximum, 
which  grows  with  the  time  delay,  which  in  turn  grows  with  k 
Decrease  of  the  time  delay  by  postponing  the  fast  change  in  L magnitude 
and  phase,  now  at  (1,000-3,000)  to  a higher  frequency,  would  be  very 
wasteful  in  terms  of  L bandwidth.  Actually,  it  is  found  that  a small 
increase  of  phase-margin  by  a few  degrees,  suffices  to  satisfy  the 
overshoot  limitation. 


The  second  important  phenomenon  which  does  not  agree  with  the 
second-order  model,  is  the  large  undershoot  in  the  step  response  of 
L/(l+L)  at  large  k , as  in  Fig.  III. 15  at  k=1000  . This  is  assoc- 

iated with  the  approximately  horizontal  character  of  'L^(jw)  (see 
Figs.  III. 13,  111.14)  in  the  range  w~ (1500, 3000)  . The  frequency 

response  of  !L^/(1+L^)!  at  k=1000  , can  be  obtained  by  shifting  the 

zero  db  line  to  -60  db  in  Figs.  III. 13,  III. 14  and  is  shown  in  Fig. 
III. 14.  Besides  the  principal  peak  at  , there  is  a minor  peak 

(or  almost  horizontal  segment)  at  . It  has  been  noted  [S?1  that 


I- 


3": 

the  inverse  transform  of  L^/(l+L^) 
atively  predicted  from  a sketch  of 

reflecting  the  latter  about  the  vertical  axis.  The  shape  of  ! 

at  predicts  an  increase  in  the  slope  of  the  inverse  transform 

followed  by  a decrease,  as  compared  with  the  time  response  if  M2 
was  not  present.  This  prediction  is  borne  out  in  Fig.  I I I. 16b,  which 
gives  the  derivatives  (impulse  responses)  of  the  two  curves  in 


can  be  approximately  and  qualit- 
:L  ( j«)/(l+L  (j«>)  I vs  « , by 


Figure  III. 16.  The  two  significant  phenomena  of  Part  2 of  design. 
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Fig.  III. 16a.  has  a minor  effect  on  the  step  responses  in  Fig. 

III. 16a,  but  a major  effect  on  their  derivatives  in  Fig.  III. 16. b.  The 
larger  peaking  in  the  impulse  response  (due  to  , gives  the  under- 

shoot observed  in  Fig.  III. 15  at  k=1000.  If  one  uses  these  impulse 
response  curves  together  with  Eq.(3.1),  it  is  found  that  the  nonlinear 
step  response  obtained  is  precisely  as  predicted  by  Eq.(3.1). 

The  above  explanation  sugae  hat  the  minor  peak  in  Fig. 
ITI.14  should  be  eliminated,  by  using  a characteristic  such  as  or 

in  csgs.  III.  13,  III.  14.  The  larger  undershoot  in  the  step  response 
was  thereby,  in  fact,  considerably  decreased  as  can  be  seen  on  Fig. 

III. 17.  We  have  gone  into  this  matter  in  great  detail,  because  of  the 
universal  nature  of  these  results:  Use  of  nonlinear  compensation  is 
justified  and  highly  worthwhile  when  large  high-frequency  plant  gain 
factor  uncertainty  leads  to  a large  frequency  range  over  which  L(j“)! 
(in  . linear  desicn.)  must  decrease  relatively  slowly.  Nonlinear  design 
per, .-its  a significantly  larger  rate  of  decrease.  The  loop  character- 
istics have  a universal  nature  in  the  high-frequency  region,  so  it  is 
important  to  thoroughly  understand  their  properties  in  the  nonlinear 
design . 

III.3.e.  Detailed  Design  - Part  II. 

While  the  preceding  explain  the  results  obtained,  they  do  not  by 
themselves  provide  the  data  needed,  (equivalent  to  Fig.  III. 10  for  the 
first  part  of  the  design),  for  quantitative  design  to  specifications. 

For  this  purpose  a model  is  used  to  represent  as  accurately  as  possible 
the  system  step  disturbance  response  and  to  provide  the  mass  of  data 
that  is  needed  in  order  to  permit  the  desianer  to  systematically  find 
satisfactory  design  parameters  for  his  specific  problem.  Dependinc 


X’W’%  v- 


’VI" 

r- 
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upon  the  nature  (odd  or  even)  of  the  excess  of  poles  over  zero  e 


use  can  then  be  made  of : 

(s+z) 


u„4n 

1 


(3.12)  Lft(s) 


s2(s2+2C)li1)fs+aj)l2)2n 


(s+z) 


if  eT  is  odd 
Li 


s2(s+u)l)  (s2+2C^io^s+Uj2)  2n 


u>„4n+l 

„ l if  e.  is  even. 

U L 


9 being  such  that  9 = w where  w is  the  crossover 

cmax  cmax 

frequency  of  the  L (s)  which  corresponds  to  k = k 

eq  r max 

In  this  'part  II'  frequency  range,  it  suffices  to  obtain  a 

satisfactory  design  at  k = , corresponding  to  the  zero  db  line 

shifted  downward  by  60  dbs  (Fig.  III. 13,  III. 14).  This  is  because  the 

two  phenomena  previously  noted  have  their  maximum  effects  at  kmax 

It  is  upto  the  designer  to  make  sure  that  his  actual  L (s)  (at 

eq 

k = k ) approximates  L^ts)  quite  closely  for  a decade  or  so  on 
max  A 


each  side  of  the  crossover  frequency  u> 


If  this  is  not  so,  the 


model  L^fs)  should  be  changed  accordingly.  Without  loss  of  gener- 


ality, following  [H2]  we  select  e =5  , i.e.  n = 2 

L 


From  Fig. 


2 4 

III.  13,  9 = «»  32.10  . The  design  parameters  are  presented  in 

Tables  III.l,  III. 2.  (1)  The  CP  parameter  is  related  to  the  peaking  in 

L,  i.e.  as  i changes  from  1 to  4,  the  nature  of  L changes  in  the  manner  in 
Fig.  III. 14  for  L, 


at  k = k 


j = l to  3.  (2)  G is  the  gain  margin  of  L. 

m l 

(3)  (P  is  the  phase  margin  at  the  point  at  which 


Arg  begins  to  rapidly  decrease.  There  are  three  data  values  given 
for  each  combination  of  the  above  parameters.  These  are  (a)  OV,  being 
the  first  peak  overshoot  value  as  well  as  the  second  and  third,  when 
the  latter  two  are  also  significant  (b)  UN,  the  first  peak  undershoot 
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Table  III.l. 


Cl 

c2 

c3 

c4 

G = 

m 

10 

dB,4  =2.5® 

ov 

40,77,90 

30,44 

29,32 

27 

UN 

62,66,70 

42 

29 

20 

GF 

0.39 

0.68 

1.04 

1.6 

G * 
m 

12 

dB,  0 - 2 .5® 

OV 

23,48,50 

25,33 

23,24 

24 

UN 

42,35 

28 

22 

8 

GF 

0.65 

1.07 

1.55 

2.2 

G = 
m 

14 

dB,  4 = 2 .5° 

OV 

26,42,38 

26,27 

24 

18 

UN 

26 

18 

12.5 

3.5 

GF 

1.06 

1.80 

3.17 

9.20 

G * 
m 

10 

dB , 4 = 5 ° 

OV 

28,64,79 

35,57,58 

31,38 

30 

UN 

61,61,62 

52 

42 

22 

GF 

0.68 

0.72 

1.02 

2.00 

G * 
m 

12 

dB,  4 = 5 ® 

OV 

18,43,46 

22,37 

27,28 

21 

UN 

40 

38 

28 

13 

GF 

1.07 

1.16 

1.60 

4.80 

G = 
m 

14 

dB,  4 = 5° 

OV 

26,40 

21,28 

21 

20 

UN 

28 

28 

21 

4 

GF 

1.67 

2.00 

2.92 

12.3 

G * 
m 

10 

dB,  4 = 10  ° 

OV 

15,39,44 

27,36 

26 

UN 

43,34 

50 

20 

GF 

2.12 

1.63 

4.61 

G = 
m 

12 

dB,  4 = 10° 

OV 

15,33,33 

22,22 

19 

UN 

40 

34 

14 

GF 

2.38 

2.72 

10.4 

G * 
n 

14 

dB,  4 • 10® 

OV 

20,29 

17 

13 

UN 

37 

30 

10 

GF 

5.01 

4.74 

28.8 

as  well  as  the  second  and  third  when  significant  (c)  GF,  the  normal- 


zed  high  frequency  gain  factor  of  LeQ(j  to)  . (Multiplication  by  9 
gives  the  real  gain.)  The  .-mailer  GF  , the  more  economical  is  the 
design. 

It  is  necesssary  to  relate  the  design  parameters  of  Table  III.l 
to  design  numbers,  e.g.  poles  and  zeros  of  which  correspond  to  any 

entry  in  Table  III.l.  This  information  is  available  from  Table  III. 2 
which  relates  the  entries  and  parameters  of  Table  III.l  to  those  of  the 
loop  transmission  of  (3.12). 

For  our  design  example  let  us  assume  that  in  addition  to  20%  over- 
shoot, (design  ) , a maxim-urn  of  30%  undershoot  is  tolerable,  in  the 
system  step-disturbance  response.  Using  Table  III.l,  a phase  margin 
tp  = 2.5°  is  then  first  considered.  A satisfactory  design  can  be  obtain- 
ed with  the  combination  C.  and  G = 14  db  , however  there  is  an 

4 m 

overdesign  with  respect  to  undershoot  specifications.  At  D = 5*  , it 

is  seen  that  the  combination  C,  with  G = 14  db  slightly  violates 

3 m 

the  overshoot  specification,  but  the  undershoot  is  here  too  only 

21%  <30%  max.  The  last  design  is  certainly  more  economic  because  its 

GF  “3.10^  , compared  to  9.1016  at  ip  = 2.5°  . If  ip  = 10°  , it  is 

seen  that  combination  C,  for  G = 12  db  , is  almost  satisfactory. 

3 m 

It  corresponds  (Table  III. 2)  to  the  parameters  = .29  , z = 1500  , 

2 7 

to,  = 1.13  10  . By  cut  and  try  it  is  found  that  the  design  for 

2 7 

which  = .29  , z = 1750  , u>^  = 1.15  10  satisfies  exactly  the 

specifications,  as  shown  on  Fig.  III. 17,  (curve  1). 

If  only  a very  small  undershoot  (Q<0%)  (design  ) is  tolerable,  it 
can  be  found,  following  the  same  procedure,  that  the  best  design  lies  around 
crr-jbir.ation  , with  to  = 2.5°  and  G^  = 14  db  , corresponding  to 
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2 7 

Cj  = .50  , = 2.4  10  and  z » 2000  . By  cut  and  try  the  design 

L2  is  then  obtained  with  C.  * .53  , » 2.25  107  , z » 1750 

for  which  the  step  disturbance  response  satisfy  20%  overshoot  and  0% 
undershoot  specifications,  as  can  be  seen  in  F;g.  III. 17.  If  we  want 
an  r,j  design  for  which  the  step  disturbance  response  is  similar  to 
the  one  shown  in  Fig.  III. 17,  the  same  procedure  gives  “ .65  , 

■J  = 2.5  107  , z = 1750  . 

The  three  different  designs  that  were  obtained  (Fig.  III. 13, 

III. 14)  are: 


. k(s+.905) (s+.92) (s+17.) (s+17. 3) (s+26.1) (s+200) (s+1750) 

L (S)  « u . 

s(s+512)  (s+5.28)  (s+6.09)  (s+23)  (s+37.8)  (s+50)  (s+130)  (s +2t  w s+u  ) 


k * 2.42  1013  , uj  - 11.5  106  , 2?^  - 1950  when  i«l 

k « 9.2  1033  , Wj  * 22.5  106  , 2C^(i)1  » 5000  when  i«2 

k * 11.4  1033  , a.'2  * 25.  106  , » 6500  when  i = 3 

The  prefilter  for  all  three  designs  is,  not  unexpectedly,  the  same 

as  the  one  used  in  the  linear  time  invariant  design  (see  II. 2 for 
details) . 

III. 4.  Design  results. 

The  responses  to  command  inputs  for  different  plant  parameter 
combinations  are  shown  in  Fig.  III. 18  for  the  nonlinear  design.  Not 
unexpectedly,  these  hardly  differ  from  those  obtained  in  the  L.T.I. 


design  (see  [H2,  S2] ) . Recalling  II. 2,  let  F(s)  = 


a=  3 rps  as 


a crude  approximation.  Using  a delayed  second  order  approximation 


(Appendix  Al)  for  the  closed  loop,  we  get  (Fig.  III. 7): 

2 

a)  -stj 

. , n e a a 

* ~ 2 2 s(s*ij  when  the  command  is  a unit  step. 

s +2 Cu  s+u 
n n 


«2,-7 


2 l-7 

V° 


zlO-3 

2^-7 

V° 
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Table  III. 2. 
C2  Cj 

= io  dB,  $ =; 


12  dB,  <J>  =2.5° 


= 14  dB,  4>  i 


G 

m 

= 10  dB, 

* =5° 

zlO-3 

3 

2.5 

2.0 

1.5 

2 TO-7 

v° 

0.12 

0.8 

0.165 

0.75 

0.23 

0.8 

0.37 

0.97 

G 

m 

= 12  dB, 

4 =5° 

zlO-3 

3 

2.5 

2.0 

1.5 

2"  1-7 

V° 

0.14 

1.0 

0.185 

0.95 

0.26 

1.0 

0.462 

1.5 

G 

m 

= 14  dB, 

* =5° 

zlO-3 

3.0 

2.5 

2.0 

1.5 

2>-7 

V° 

0.155 

1.25 

0.21 

1.25 

0.3 

1.35 

0.59 

2.4 

G 

m 

= 10  dB, 

4 « 10° 

zlO-3 

2.0 

1.5 

1.0 

2'’* -7 
U?.10 

0.16 

1.15 

0.25 

0.875 

0.53 

1.2 

G 

m 

= 12  dB, 

4 =10° 

zlO-3 

2.0 

1.5 

1.0 

2in-7 

U).  10 

Z 

0.185 

1.22 

0.29 

1.13 

0.65 

1.8 

G 

m 

= 14  dB, 

4 =10° 

zlO-3 

2.0 

1.5 

1.0 

C0 

2 l-7 
u)‘lO 

0.195 

0.77 

0.325 

1.49 

0.835 

3.0 
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for  t,  < t < t„ 


cT(t)  = -X*  = - 


2a 


, — 2 2 
u Sl-C  (l -c  ) 
n 


'l'sin(u)  /T-C  (t-t, )) 
n X 


2a 


"n(t'tl)  2 

sin(u  /l -£  t) 
n 


Therefore,  the  nonlinear  system  response  is,  for  t^  < t < 


CNL<t)ftSl-e 


■Qt  t Ci  / — 2 

+ - - sin(u)  /1-C  t)  (3.14)  because  C is  very  small. 


n 

Noting  now  that  X * » X*  for  K = 2,3,...,  it  is  seen  that  the  nonlinear 

K 1 

-ci  t 

system  response  is  C (t)  «s  1-e  for  t not  too  small,  i.e.,  when 
NL 

a/o)  can  be  neglected  in  comparison  of  1-e  at  . For  instance  at 
n 

k = 1000  , to  »600  rps  and  V to  » 1/200  . Therefore  when  t>  .04 

n n 

seconds  (compared  to  a settling  time  of  approximately  1 second) , 

1-e  at  > .1  = 20  a/to  and  the  sinusoidal  component  can  be  neglected. 

n 

However  for  small  t ( < .04  seconds)  this  is  not  true  and  the  nonlinear 
and  linear  system  responses  differ  one  from  another.  Fortunately,  this 
cannot  be  noticed  because  all  the  quantities  involved  are  then  very 
small. 


For  intermediate  values  of  the  gain  factor  (SfkflOO)  we 

have  30  rps  < < 100  rps  and  . 5 < b/u^  < 1 . 66  giving  (Fig.  III. 5) 

at  most  a reset  value  of  30%  of  the  reset  value  obtained  with  a C.I. 

in  place  of  FORE,  which  can  once  more  be  neglected  because  here  too 

>>  a . 
n 

For  small  gain  factors  (k  < 5)  the  linear  system  is  similar  to 

the  nonlinear  one  because  w <30  ms  and  therefore  >2  giving 

n n 

(Fig. III. 5)  a reset  value  which  is  very  small  ( < 5%)  and  can 
therefore  be  nealected. 
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Assume  now  that  a command  input  is  applied  whose  bandwidth  is 

smaller  than  the  bandwidth  of  the  prefilter.  At  small  gain  factors 

oiving  b/u  > 2 the  reset  value  ( < 5%  of  the  C.I.  reset  (Fig. III. 5)) 
n 

can  certainly  be  neglected.  At  large  gain  factors  giving  bandwidth  of 

loop  » bandwidth  of  prefilter,  we  have  by  anology  with  (3.12a)  and 

(3.13)  C ( t) w command  input  + A sin(ui  /f-C^t  + 9)  where  the  residue 
£ n 

i F( jw  ) R(jcu  ) | 

A is  =21  — jm-j |.  Therefore  A is  very  small  as  before, 

J n 

and  the  arguments  made previously  are  easily  extended  here. 

At  intermediate  gain  values,  it  might  happen  for  some  problems 

that  the  corresponding  range  w is  not  big  enough  compared  to  the 

n 

bandwidth  of  the  prefilter,  preventing  us  therefore  from  making  general 
statements.  However,  for  these  problems  (as  in  the  present  case)  where 
<o  » bandwidth  of  the  prefilter,  even  in  the  intermediate  range  of  plant 
gain  factor  values,  the  previous  assertions  can  be  extended. 

We  can  therefore  conclude  that  the  superposition  theorem  holds 
effectively,  at  all  gain  factors  for  command  inputs,  provided  that  the 
bandwidth  of  the  command  signal  f bandwidth  of  the  prefilter.  It  is 
probably  unlikely,  in  most  cases,  that  the  system  will  be  subjected  to 
command  inputs  whose  bandwidth  exceeds  3 times  the  bandwidth  of  the 
prefilter,  (here  3 w «j  10  rps)  because  the  amplitude  of  the  command 
signals  at  the  output  of  the  prefilter  would  then  be 'very  small  ( < 1/10 
of  the  input  command  amplitude  in  the  present  case  because  (see  II.  2) 

! F ( j 10) ! « .1  ) . 

The  step  disturbance  responses  are  shown  on  Fig.  III. 19. a & b at 
different  plant  parameter  values  for  the  L.T.I.  design  and  in  Fig.  III. 
20.  a-c  for  all  N.L.  designs. 


lb 
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In  all  cases,  the  overshoot  is  S to  the  maximum  tolerable  20%. 

Comparing  Figs.  III. 19. a and  III. 20. a,  one  notes  that  besides  the  f 

case,  there  exists  only  slight  differences.  This  is  expected  because  for 

small  gain  factor  values,  giving  b/w  > 2 , FORE  acts  like  a linear 

n 

element  (Fig.  III. 5).  It  is  then  easily  understood  why  there  exists 
only  minor  differences  between  the  and  designs  of  Fig.  III. 

13  and  14  for  those  plant  conditions.  Therefore  we  can  say,  that  super- 
position holds  for  small  gain  factor  values  for  disturbance  inputs  whose 
bandwidth  is  smaller  than  the  bandwidth  of  the  loop.  When  the  latter  is 
untrue,  C (tl^dlt)  because  the  feedback  is  ineffective.  Therefore 
r(t)«)-d(t)  and  the  first  reset  is  determined  by  the  first  zero  cross- 
ing of  d(t)  . Let  d(t)  “Asinw.t  with  u>  >>  oj  by  assumption. 

„ d d n 

n/u 

° 2An 

Therefore:  x*  » - / d(£)cU  n « - (where  n is  the  reset  value 

1 ui. 

0 d 

percentage  given  in  Fig.  III. 5 as  a function  of  b/w  ),  and 

n 

CO  _ , 0y 

C _ (t)  *»d(t)  + 2An  — e Un  sin(u  /l-?  t + 0 ),  or  eouivalently , 

\L  u n 

d 


C.,r(t)  «d(t)«CA(t)  because  o> . >>  u 
NL  £ d n 


Therefore  superposition  also 
holds  at  small  gain  factors  when  the  bandwidth  of  disturbance  >>  band- 
width of  the  loop.  The  differences  that  exist  for  high  gain  factors 
between  Fiq.  III. 13. b,  Fig.  III. 20. b and  c have  already  been  explained, 
and  were  therefore  expected,  so  the  designer  need  only  consider  the  non- 
linear response  at  high  gain  factors  for  disturbances  other  than  steps. 
For  instance.  Fig.  III. 21  presents  responses  at  large  k to  a set  of 

disturbances  d(t)  = e -lu,dtcoswdt  , «.  =*  30,300,1200,6000  rps. 

s 


Using  the  approximation  D(s)  « 


we  have : 


2 2 
s +U£j 

Q(s)  £ t~“  D»-  w 2s  e"Std/«s2+(o2)  (s2*2Coj  s+u  2))  with  t.  < .15  for 
1+L  n a n n 

high  gain  factors  as  mentioned  above.  Therefore:  (t  = TIME  - t ,) 


1.00 


(a) 


TIME 


Figure  I I I. 20.  Response  to  step  disturbance  (nonlinear  design 
-.j ( t ) = +A  e r'UJntsin(u  /l-£2t  + 0 ) -Bsin(w,t  + g ) 


with  A = w 


we  have 


t + cos  C)  which  is  the  L.T.I.  system 


step  disturbance  response  when  using  a second  order  approximation 


30  rps  (compared  with  w = 6000  rps 


disturbance  (compare  Figs.  III. 21. d and  III. 20. b,  case  f ) . Therefore 


the  reader  notes  that  when  the  bandwidth  of  the  loop  » bandwidth  of  d*st 


inputs  the  superposition  theorem  holds  even  for  large  plant  gains 


course,  the  disturbance  attenuation  capability  decreases  as  w 


6000  rps 
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then  the  feedback  is  ineffective  (A« — (— ) 2 , B (s  (— )2)  and 

/IV  ud  wd 

C ( t)  ta D ( t)  = coa)  t . Therefore  i.  (t)  = -cosa)  t and  the  first 
D t X.  a x,  q 

reset  is  determined  by  the  first  zero  crossing  of  D(t) . Therefore, 

X * » 2/  u>  and  I*  C_(t-t  )w  — e (’“n^t  tl^sin(a)  /l^2(t-t.))  . This  is 
id  1 l 1 co  n 1 

d 

obviously  a very  small  quantity  which  can  be  neglected  in  general. 

Fig.  III. 22  presents  the  system  output  for  simultaneous  command 
(unit  step)  and  the  disturbance  and  white  sensor  noise  inputs  shown. 

It  is  seen  to  be  satisfactory.  Indeed  one  notes  that  superposition  holds 
here,  in  a general  sense,  due  to  the  fact  that  the  bandwidth  of  the 
disturbance  signal  <<  bandwidth  of  the  loop. 

In  summary  we  can  say  that  the  superposition  theorem  holds  in  a 
general  sense  for  command  inputs  and  disturbance  inputs  provided  the 
bandwidth  of  the  signals  is  smaller  than  the  bandwidth  of  the  loop. 

Furthermore  it  should  be  underlined  (as  stated  many  times)  that 
the  reset  instants  are  primarily  determined  by  the  loop  capabilities 
when  bandwidth  loop  » bandwidth  signal,  while  there  are  determined 
by  the  signal  (command  or  disturbance)  when  the  converse  is  true. 

Fig.  III. 23. a, b present  the  responses  due  to  a truncated  ramp 
disturbance,  and  to  a unit  ramp,  respectively.  One  notes  that 
(Fig.  III. 23. b)  for  some  plant  conditions  a limit  cycle  occurs  in 
presence  of  a ramp  disturbance.  Much  attention  will  be  devoted  in 
chapter  IV  to  the  stability  of  feedback  systems  which  include  reset 
elements,  and  it  will  be  shown  how  such  instability  can  be  eliminated. 
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H 
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from  those  of  the  FORE  input  X(t)  , but  even  if  the  statistics  of 
T ( t)  are  known,  calculating  those  of  the  FORE  output  y^(t)  would 

I j 

require  computation  of  the  expected  values  of  a stochastic  integral, 
one  of  whose  limits  is  a random  process  whose  statistics,  in  turn, 

i I 

depend  on  the  integrated  process,  obviously  a tremendous  task.  Using 

the  fact  that  for  a zero-mean  Gaussian-Markov  process  X(t)  , the 

mean  and  variance  of  T tend  to  zero  as  the  process  becomes  more 

uncorrelated,  following  references  [PI] , T is  approximated  by  its 

mean  value  T giving: 
m 


rfc  -rirs  -Mt-O  \ Wr,  -b(t-C) 

y,(t)  = S X(C)  e dC  “ / xf 5 ) e d? 

f t-T(t)  t-T 

m 

The  cross-correlation  of  x(t)  and  y(t)  is 

t+T  -b(t+T-C) 


<p  (t)  = E{y  (t+T)x(t) } - E { / X(?)e 

yx  f - 


X(tldC> 


t+r-T 


m 


T 

"I  , , b^p-Tm'  , • Taking  Fourier  Transform  gives  [Gl] 

= / <J>  (p+T-T  )e  dP 

0 XX  m 

-T  (b+ jw) 

« ^(w)  = r-r-  (1-e  )=  («)M(«,t.) 

yx  XX  b+]W  XX  m 

The  normalized  random-input  describing  function  is  therefore 


N (<u,T  ) 


N(oj,T  ) 
in 


n ' ‘m  [1/ (b+jcu) ] 


1-e 


-T  (b+jw) 
m 


(3.15) 


so  that  IN  1 >1  over  some  w and  <1  over  others. 

n 

If  the  input  X to  FORE  is  a stationary  white  noise  process  of 
2 

power  spectrum  “ °,j  > then  the  root  mean  square  value  of  the 


noise  output  y^  is  given  by: 

+a> 

(y,>  * 157  / lN(",T  >! 2 o 

f rms  L2it  ' m 


\ ♦**.-**.  ■ 


r + cd  7 MT  , 1L 

v II  ....  -bTm.  . -bTm  2 m,  1 I ‘ 

(yfW  ; [(1+e  ’ - 4 e COS  ~ d“  I °N 

L — op  zh  +uj  J w 


'Vras  = [ ^ <l+e"bTm)2  - ^ f ^"1 


-bT„ 


_ ojT 

2 m 


b +oj 


cos  — 2 du 


°N  (3 


Let  us  consider  the  signal  f(t)  of  Fig.  III. 24,  we  have 


, 1 , (s+b)  T /2  -(s+b)  T /2  , 2 ...  ..  , m 

F(iai)  = — — (e  m - e m ) = — — sinh(s+b)  — 

S+D  S+D  2 


+m  2 i *- 

Using  Paseval's  relation:  f f (t)dt  = — / F(jw)  F(-jw)du 

ZTT 


2*  ' 


T T 

yy  sinh(b-jw)  ysinh  (b+ju)  y do) 


00  'ji  +b 

<o 

4 


— / 

2ir  2 .2 


-=>  a)  +b  '■ 
bT 


, bT  - ooT 

m 2 m 

cosn  — r — - cos  — — 

2 2 


du 


= £.  cosh 


m 


, +~  . 2 _ wT 

1 , 4a  2m, 

I y — g cos  2 d(,l) 


2it 


(3. 


-oo  (jo  +a 


•H® 


bT 


bTm  “bTm 

2 0 m _ g m 

But  / f (t)dt  = ^ =sinh  b 

—00 


m 


(3. 


Using  (3.16a)  and  (3.16b)  in  (3.16)  we  get: 

T 1 -bT  2 2 bTm  ebTm  1 c 

= ■—  (1+e  bT"V  +(sinhbT  - 2cos)i  -r~-  ) — r—  N 


.16) 


16a) 

16b) 


(VJ 


and  therefore: 


ty,) 
f rms 


’ (1+e  2 2e”bTm  _ 12 

= + —5—  (smh  bT^  - cosh  bTm  - 1)  ! ffN 


or  <yf>  = 0 jr-  (l-e'bTm) 
I IT'S  N 2d 


(3.17) 


For  the  linear  element  1/s+b  , it  can  be  shown  that: 


/T 

^lin^rms  °N  2b 


(3.18) 


It  can  then  be  seen  from  Fig.  III. 25  that  in  this  case  the  noise 

level  at  the  output  of  FORE  is  always  smaller  than  the  noise  level  at 

the  output  of  the  linear  element. 

The  theoretical  results  are  in  accordance  with  those  obtained 

experimentally  by  controlling  the  interval  between  resets.  It  is  noted 

that  when  T^b > 1 by  a factor  of  3 or  more,  i.e.  when  the  noise  signal 

is  relatively  slow  varying,  then  there  is  practically  no  difference 

between  FORE  and  1/s+b  . However,  the  difference  becomes  tremendous 

as  T b <<  1 . Here  too  we  can  say  that  FORE  can  discriminate  between 

m 

fast-varying  noise  signals  and  slow-varying  ones. 


O 1 

f rms  f 


1 

N 2b 


. P,n  — I 

P'  2b  r 

f /r  I 

•G5k  jb  r 


*ylin^  rms 


4 /T 

-4aN  2b 


0.  .15 


T75  IT  1.25  1.5  17s 


Figure  TTI.75.  R.M.S.  values  of  the  noise  at  the  output  of  FORE 
5 1/s+b  in  the  Dresence  of  white  noise  input. 


The  effect  of  white  sensor  noise,  in  the  closed  loop  system,  at 


the  plant  input  Z is  given  by 


$ (oi) 

22 

$ (U>) 


| L/P  ;2 
1+L 


and  the  m.s.  noise 


nn 

2 °°  du) 

at  z is  o = / <?  (u>)  — , where  the  $ are  the  power  spectra. 

2 22  2TT 

If  L = N ( uj,T  ) L (s)  of  (3.15)  is  used,  then  at  some  T , the  above 
n m eq  m 

becomes  infinite,  as  shown  in  Fig.  III. 26.  However,  experimental 

2 

determination  of  a (T  ) gave  not  unexpectedly  finite  results  for 
z m 

all  T , as  shown  in  Fig.  III. 26.  In  this  run  T was  controlled 
m m 

by  having  FORE  reset  determined  externally  at  period  T , rather  than 

m 

by  the  2ero  crossings  of  FORE. 

When  the  actual  nonlinear  system  was  simulated,  the  experimental 


Figure  III. 26.  Effect  of  sensor  noise  at  plant  input  versus 
Tm  at  k = 103  . 
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result  was  a (normalized)  value  of  0.37,  corresponding  to 
z 

= (2.5)  10  ^ sec.  in  Fig.  III. 26.  The  linear  design  of  [ H2 ] for 

2 4 

the  same  oroblem  resulted  in  a (normalized)  = (1.66)  10  , both 

z 

theoretically  and  experimentally.  Thus,  the  nonlinear  design  which 

achieves  the  same  system  output  tolerances  to  command  and  disturbance 

inputs  as  the  linear  design,  does  so  with  rras  sensor  noise  effect 

/ 4 

smaller  by  a factor  of  0.37/(1.66)10  = .0047.  The  improvement  is 

very  spectacular  in  this  specific  numerical  example  because  of  (1)  the 
large  uncertainty  in  the  plant  high-frequency  gain  factor  k , which 
gave  the  nonlinear  design  a large  f regency  range  in  which  to  exploit 
its  advantage  of  larger  phase  lag,  and  (2)  the  great  difference 
between  ! L ! and  'p!  in  the  high-frequency  range,  which  leads  to 
serious  noise  effect  at  the  plant  input.  The  advantage  of  nonlinear 
over  linear  design  will  be  less  to  the  extent  that  these  two  factors 
are  lessened. 


v ***** ' 
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CHAPTER  IV. 


IV. 1.  Introduction. 


General  stability  criteria  are  known  for  only  a small  class  of  non- 
linear and  time  varying  feedback  systems.  The  circle  criterion,  for 
example,  does  not  apply  to  the  C.I.,  FORE  or  the  more  general  reset 
elements.  It  was  necessary  therefore,  to  develop  our  own  criteria  for 
these  elements. 

In  chapter  III  we  implicitely  restricted  ourselves  to  stable  open 

loop  transfer  function  L (s)  (Fig.  IV. 1. a)  and  will  therefore  derive 

eq 

sufficient  conditions  for  Bounded  Input  Bounded  Output  stability  of 
the  nonlinear  feedback  system  of  Fig.  IV. 1. a.  We  will  then  investigate 
possible  limit-cycles  in  the  nonlinear  feedback  system  of  Fig  IV. ’. a. 
containing  the  element  (1/s+b)*  whose  output  y is  reset  to  zero 
whenever  the  input  X equals  the  real  value  a ; a = 0 , b ^ 0 corres- 
ponding to  FORE  and  a =0  , b = 0 to  C.I.  Only  stable  limit-cycles  with 
at  most  two  reset  instants  per  cycle  will  be  considered.  More  general 
limit  cycles  can  be  considered,  however  todate,  limit-cycles  with  more 
than  2 resets/cycle  were  never  encountered  experimentally. 

Our  primary  objective  will  be  to  ensure  stability  when  FORE  is 
used,  because  of  its  usefulness,  demonstrated  in  chapter  III.  In  the 
following,  intensive  use  of  the  2-transform  is  made  with  the  usual 


notations 


'(s)  ,•*(),  2 ( • ) 


Note  that  XK*  introduced  in 


chapter  III  is  a number.  In  this  chapter,  9 will  denote  the  set  of 
reset  instants,  i.e.,  9 * {t  3 x(t  ) = a}  . 

K K 


srws  -s- 


Figure  IV.l.b.  Equivalent 

representation . 


TV. 2.  Bounded  Input  Bounded  Output  (B.I.B.O.)  stability 

sufficient  conditions  for  a class  of  nonlinear  feedback 
systems . 

Lemma  ) : Consider  the  nonlinear  feedback  system  of  Fig.  IV. 1. a, 
where  L^ts)  denotes  the  open  loop  transfer  function  of  the  L.T.I. 
system  obtained  in  absence  of  resets.  If  the  L.T.I.  feedback  system 
is  asymptotically  stable  and  if  the  sequence 
fcK  -b(t  -O 

{jx  * ! s.  | / XtOe  di;  | , t , ,t  € 9}  of  the  nonlinear  system 

K-l 

is  bounded,  then  the  nonlinear  feedback  system  considered  when  resets 
occur  is  stable  in  the  B.I.B.O.  sense. 

Proof : It  was  shown  (equivalence  of  Fig.  III. 7 and  III. 8)  that 

Fig.  IV. 1. a.  is  equivalent  to  Fig.  IV.l.b. 

00 

Therefore  CCt)  = C^t)  +CfJL(t)  = C (t)  - £ j^*C,(t-t  ) (4.0)  , 

K — 1 

where  C^(t)  denotes  the  L.T.I.  system  impulse  response.  By  assumption 
the  LTT  system  is  asymptotically  stable.  So:  V t , 3 , a € IR  such 

that  ! ( t ) | f M e 'lt  (Al)  and  (B.I.B.O.  of  the  LTI  design 

which  follows  from  asymptotical  stability)  VR  such  that  | R ( t ) ! t 
3 M such  that  :Cj ft)  < M (A2)  (see  Fig.  IV.l.b). 

Let  us  first  assume  that  the  set  9 Qf  reset  instants  is  finite, 
Then,  if  y § sup  and  usi- g (Al)  we  have:  for  t > t„ 


Thus  the 


lemma  i 


Assume  now  that  9 is  infinite  but  countable.  By  assumption 


with  t 


9 would  be  finite,  while  if 


would  not  be 


We  have  then 


£ 9 such  that  t € [t 


Therefore  VK,  (using 


Using  then  (4.0)  and  (A. 2) 


Stability  is  therefore  proven  when  9 is  countable 


Assume  now  that  9 is  not  countable.  It  means  that  there  exists 


at  least  one  instant  t 


By  definition 


on  which  x(t)  E 0 


1 
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t + r. 

X*  = lin  ! X(!;)e"b(tP-0dC  = 0 
P1  c=0  t 


The  contribution  to  C„, (t) 

NL 


(Fig.  IV.l.b)  of  the  interval  ]t  ,t  + n]  is  therefore  null  and  the 

P P 

contribution  to  C (t)  of  Q is  therefore  the  same  as  a single 
NL  *d 

reset  occuring  at  time  t = t^  . A new  equivalent  set  9 is  then 

obtained  by  removing  all  the  possible  intervals  from  9 and  by 

replacina  them  by  the  time  instant  t . Obviously  9 is  countable 

P 

if  and  only  if  {Q  } is  countable.  Given  two  intervals  n and  C) 
p "s  t 

we  have:  p^  n Q = 0 for  all  s , t and  U c TR  • Using  the 

fact  that  between  two  real  numbers  there  always  exists  a rational 

number,  each  Q is  associated  with  one  rational  number.  The  rational 

P 

numbers  being  countable,  so  is  {Q  } . The  lemma  being  true  for  the 

P 

countable  set  9 we  have  proven  the  lemma  for  an  uncountable  set  9 

Lemma  1 rests  on  the  assumption  of  a bounded  sequence  { ! X * } . 

x * and  t are  defined  by  the  set  of  equations: 
mm 

m-1 


x(t ) — x (t ) + y x * c.  (t  -t.)  = o 

m 2.IH  “ V A m V 


K=1 

t , t 

= / xp  (c)e-b(tm-,;)dc  + l f e"  { -m-55  Cfi  (?-tK)d!; 

, Nn-l  K~  1 ^-1 


6'  m V 


with  X (t)  ^R(t)-C.  (t) 
1 l 


In  absence  of  a rigorous  proof  we  conjec- 
ture that  is  bounded  if  the  LTI  feedback  system  is  asymp- 

totically stable.  It  should  be  emphasized  that  once  the  above  is 
rigorously  established,  fruitful  results  on  stability  will  follow  such 
as:  a)  generalization  of  lemma  1 to  any  reset  element  g*  (defined 
in  chapter  V) , because  g*  will  be  seen  to  be  a finite  combination  of 
element  (1/s+b)*  . b)  asymptotic  stability  of  the  nonlinear  feed- 

back system. 


V* 


IT 


i. 


■ * i 
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IV. 3.  Investigation  of  limit  cycles  with  one  reset/cvcle. 
Sections  IV. 3, 4 are  exclusively  devoted  to  such  limit  cycles, 
even  when  not  explicitely  so  stated.  Recalling  (3.0a)  the  system  of 


Fig,  IV. 2. a is  described  by: 

y (t)  = / e'b(t'ax(OdC  - I [ / 


-b(t-C) 


K=1  Lt 


i(Od;  u(t-tK)  {4-1) 


K-l 


with  t € 8 = {t  :X(t  ) = a} 

1\  lx  K 


Equivalently: 

t 
/ 
t 


y(t)  = / e‘Mt‘UI(Od(  - [ 

k=l 


/e‘b(t'0x(Odcl  u(t-t  ) 


- I 

K=1 

t 

y (t)  = / e 
b0 
t 


u(t_V 


(4.2) 


-b(t-C) 


X(5)d?  - l 
K-l 


f e-b(t"°X(i;)dcl  Fu(t-t  ) -u(t-t  1 
j | K K+l)j 

0 

t (4.3) 

K 


y (t) 


/ e_b(t_!:)x(^)di;  - l e"b(t-tK)  ; e-b(tk-C)x(c)dCru(t-t  )-u^t-t  ,)1 


K— 1 


K'  K+^ 

(4.4) 


(4.5a)  , i.e.,  y? (tR)  is  the 


Let  y*(tK)  £ f e"b(tK-J:,X(C)di; 
b0 

output  at  t = t of  (1/s+b)  due  to  the  input  X(t)  (Fig.  IV. 2. b), 


Therefore: 


y (t)  =y,(t)  - l e_b(t‘tK)y|l(tk)  futt-t^)  - u(t-tK+1)  ] (4.5b) 


K— 1 


K+l 


J 


Note  that  y (t  ) [u(t-t  ) - u ( t-t  ) ] can  be  considered  as  the  output 

Z lx  K K+l 

of  a sampler  followed  by  a zero  order  hold,  with  sampling  times  corresponding 
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to  the  reset-instants  t E 9 « when  the  input  iSy^(t)  , (Fig.  IV. 2. b). 
The  output  of  the  hold  is  then  multiplied  by  e"b(t-tK*  which  is  the 
impulse  response  of  1/s+b  . Therefore  by  taking  Laplace  transform, 

we  have : 

BJs)  A [u(t-tv)  - u(t-tK+1)]  e-b(t_tK) j = (B0<S)  * i^b  ) 


1-e 


“SOu 


where  3Q(s) 

zero  order  hold  and  a A t -t 

K K K- 1 


, i.e.,  BQ(s)  is  the  Laplace  transform  of  a 


Therefore:  B (s)  = 


c+j° 


1-e 


-(S-C)Ov 


and 


B (s)  = 
e 


2 it  j f 

J c-}“ 

1-e-(s+b)CTK 

s+b 


s-L 


L+b 


dC 


(4.6) 


(4.7) 


Relation  (4.5b)  can  therefore  be  considered  as  the  sum  of  two  signals 
as  in  Fig.  IV. 2. b,  which  is  equivalent  to  Fig.  IV. 2. a from  the  input- 
output  point  of  view. 

If  a limit  cycle  with  one  reset  per  cycle  occurs,  then  as 
t “ , a -*■  T , a constant  period.  Therefore  in  the  steady  state,  the 
nonlinear  feedback  system  of  Fig.  IV. 3. a becomes  equivalent  to  the 


Figure  IV. 3. b.  Fquivalent  representation  when  a limit  cycle 
with  one  r^set/cycle  occurs. 


i 'l 

t 


L 
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linear  samplec.  data  feedback  system  shown  in  Fig.  IV. 3. b.  The  real 
number  X^*  was  defined  in  chapter  III  as:  x£  k ^ e b*tK  ^x(i;)df; 


K-l 


Under  the  assumption  of  a limit  cycle  of  period  T , 


KT 


X * = / e KT  r,'x(C) di;  (4.8).  Recalling  the  equivalent 

(K-l ) T 


representation  of  Fig.  III. 6,  xR+  represents  the  strength  of  an 
iipulse  whose  effect  is  the  same  as  the  reset.  Therefore,  if  a limit 
cycle  sustains,  all  the  reset  values  are  equal  and: 

. . = X*  = u (4.9),  where  u is  independant  of  K 


jr*  = r* 

K K-l 


What  relation  then  exists  between  XK*  and  e (kT)  ? We  have: 


e* (z)  i Z[  l e (kT)  6 (t-KT) ] (4.10).  However,  from  (4.5a) 


K=1 


e (kT)  = y^ (kT)  by  definition.  Replacing  (4.5a)  in  (4.10)  gives: 


KT 


r K. 

:*(z>  = Z J [ / 
tC=1  0 

-.[!< : 

k-l  0 


e”b(KT‘C)x(q)dc!  fi(t-KT)  j 


.-blk-1,T  • ” * ... 


yT 


(K-l) T 


(4.11).  Using  (4.9)  in  (4.11)  gives: 


’ (z)  = zf  [ y (l+e”bT+. . .+e'b(K-1)T)  6 (t-KT) j (4.12a) 


■ K=1 


-b(m-l)T 


l [ ue  •'*  6 (t-KT)  i (4.12b)  and:  e*(z) 


uz 


K=1  m=l 


(z-1)  (z-e_bT) 


(4.13) 


Let  us  now  consider  the  system  of  Fig.  IV. 3. b.  We  have: 


(1+L  (s))  E(S)  = R(s)G(s)  — —7”  + B (s)  L (s)  e*  (s) 

eq  s+b  e eq 


or  equivalently: 

R(s) G ( s) 


L (s) 


€ (s) 


J!3_ 


(1+L  (s) ) (s+b)  + Be(s)  1+L  (s) 


t*(s) 


(4.14) 


eq 


eq 


STOfc'V  IT 


1 
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Taking  the  z-transform  gives: 


E*(s) 


As  B (s) 
e 


_ (. 


RG 


1-e 


V(l+L  )(s+b) 
eq 

-(s+b)T 


s+b 


, -1  -bT. 

= l-(l-z  e ) 


and  (4.14a)  becomes: 


)*  + (Be  e‘(s)  (4-14a)  • 

eq 

_ f^eq  V _ -1  -bT  Z'  Lec  y 

’ 1 \ 1+L  / 1 ) V (s+b) ( 1+L  )J 


( 1 1 V M -1  "bT,  ( 1 \ * 

Vs+b  ' (i+L  )(s+b>y  ~ 1 e ' Vd+L  ) (s+b) y 
' eq  ' eq 


i'.C 


e*(s)  = 


yw  Va^eq)(s+b)  ; 


(d+L^Hs+b)  ) * 


(4.15) 


e*(s>  4 


z-e 


-bT 


A*  (s) 


with: 


RG 


A*  (s)  A VU+L^Hs+b) 
1 


(1+L  ) (s+b) 

eq 


r 

r 


(4.16) 


(4.17) 


As  X( s)  = (s+b)  e(s)  , (4.14)  gives  then: 

L (s) (s+b) 


r(s)  = f(s)G(?|-  + B (s)  — 

1+L  (s)  e 1+Leq(s) 

eq 


Taking  the  z-transform  gives: 


J RG 


**(z)  -fcf-)*  fz)  + A*(z,lirT-  i (z) 

eq  eq 

Comparing  (4.16)  and  (4.13)  we  can  conclude  that  A*(z)  represents 


£*{s) 


(4.18) 


the  train  of  impulses  due  to  the  resets,  and  therefore  if  a limit 

cycle  exists,  A*(z)  = f(z)  + ~rr  where  f(z)  has  all  its  poles  strictly 

inside  the  unit  circle,  and  equivalently  at  the  steady  state  A*(z)  = . 

Conversely  if,  as  t -» ® , the  reset  values  tend  toward  a constant 

y * 0 with  periodicity  T,  i.c.,  A*(z)  -»  as  z - 1 , and  this  occurs 

when  the  input  x tends  toward  a with  the  same  periodicity,  i.e., 
a 

2*  (z)  ■+  — - as  z -*•  1 , then  we  can  conclude  to  a sustained  limit 

z-l 


67 


cycle  of  period  T . We  have  thus  proven: 

Theorem  2 : A nonlinear  feedback  system  which  contains  an 

whose  output  y is  reset  to  zero  whenever  the  input 
X crosses  o , for  a given  a , sustains  a limit  cycle  of  period 
T with  one  reset  per  cycle  if  and  only  if  there  exists  some  finite 
nonzero  T , u € such  that  in  the  neighborhood  of  t = 00  , we 
satisfy: 


element 


(-L- 1 

Vs+b  > 


I 1:  m (z-1)  A*  (z) 

z=l 

IX  lim  (z-1)  r*  (z) 

z=l 


(4.19a) 


(4.19b) 


Qualitative  statements  about  possible  limit  cycles  can  be  made  by 

inspection  of  (4.17).  From  (4.19a)  a limit  cycle  exists  if 

A*(z)  -*■  — ^-r-  as  z -*■  1 . However,  the  number  of  Doles  at  z = l of 

z-l 

A*(z)  is  only  related  to  the  number  of  poles  at  s=0  of  A*(s)  , 

( RG 

or  in  other  words,  to  the  excess  e of  poles  at  s = 0 of  - — ■ . , — - . 

\ (s+b)  (1+L) 

/ i \. 

over  ^ (s+b) (1+L)  ) ' e = 0 limit  cycle  can  sustain,  while 

if  e=l  a limit  cycle  is  predicted.  4.19a  and  b give  then  the  two 
unknowns  y , T which  characterize  the  limit  cycle.  Furthermore,  as 
can  be  seen  from  (4.17)  the  value  of  e is  only  a function  of  the 
number  of  poles  at  s = 0 of  R and  number  of  poles  at  s = 0 of  L , 
under  the  assumption  that  G has  no  such  pole.  It  is  therefore  very 
striking  how  for  such  a nonlinear  feedback  system  the  stability  problem 
is  related  to  the  "type"  of  input  R applied  and  "type" of  L , where 
by  "type"  we  mean  the  number  of  integrations.  (Type  "0"  R corresponds 
to  a signal  R without  integration,  type  "1"  to  one  integration. 


)• 


etc.  . . ) 
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I V . 4 . Examples  of  applications  of  Theorem  2 . 

It  will  be  assumed  that  L (s)  has  one  integration. 

eq  / \ 

f X \* 

Therefore,  in  both  cases  (b  - 0 , b / 0)  , \i"jl7+bj~(T+L — TJ  (wl"'ich  i; 

eq  ' 

the  denominator  of  A*(s))  has  no  pole  at  z = l 


IV. 4. a. 


oe  ”0"  inputs . 


In  this  case,  e = 0 and  theorem  2 predicts  that  no  limit  cycle 


can  sustain. 


IV. 4. b.  Type  "1"  inputs . 


If  b ^ 


° '( 


(s+b) (1+L  ) 

eq 


has  no  pole  at  s = 0,  therefore  e = 0 


and  theorem  2 predicts  no  limit  cycle.  However,  if  b-0  , then 


e = l and  a limit  cycle  can  be  sustained. 


In  :he  latter,  let  G = 1 , L (s) 

eq 


A 

s(s+2o) 


R = — 
s 


We  have,  from  (4.17): 


r ! s+2c  ] 


2 L*  7 


s +20S+A 
s+2c  T 


z r-=-s-+-<g-  l 

[s  +2cts+A  J 

Suppose  that  o^  - A<  0 and  let  = /A-o^  , then: 

2 o,  ,,  -oT  _ -oT  A,sln“oT, . -2cT  -aT  A -aTSlnu0T 

-[z(l-e  cosWoT-e  (o-  )+e  ‘e  cosu^T-Ma-^e  — zr~ 


A*  <s)  = 'i ^ ZU b : 

-aT  o -aT 

(z-l)  (z  - e cosujpT  + — e sina)QT) 

and  / z(z-e  oTcosu  T+-^-  e oTsinio  T z — e °TsinwnT  A*(s) 

( 0 wo  0 ; “o  0 


x*(S)  = 


2 . -OT  _ -20T 

z - 2z  e cosWgT  + e 


2 _ -oT  _ , -2oT 
z -2z  e cos(oQT  + e 


According  to  theorem  2,  a limit  cycle  exists  if  and  only  if,  3 T>  0 , 


such  that:  [condition  II] 


lim  (z-l)  x*  (s)  = — 


-oT  . _ 

e smUgT 


Z-l  “n  i_e  aTcosoo/vT  + — e °Tsinaj„T 

0 «0  0 

In  the  case  of  a C.I.  (a=  0)  , it  is  easily  found  that: 


or  T 


v * w 


-2CT  -2CT 

1-3“  cosoj^T  + e 


and  u=  lim  (z-1)  A*( z) 


, , -OT  m o . 

z=l  1-e  cosw„T + — sinio„T 

0 „0  0 


1+e 


-OT 


If  o - A > 0 ,no  reset  occurs,  and  the  system  is  stable  if 
o > 0 , and  unstable  if  c <0  . The  reader  should  note  that  if  a > 0,  a 
reset  occurs  at  t =*°°  and  therefore  a 'limit  cycle’  then  occurs  of  period 
T = «>  . These  limit  cycles  were  not  predicted  by  using  the  Dual  input 
describing  functions  [Gl] . As  mentioned  before,  if  b ^ 0 , no  limit  cycle 
can  sustain.  Thus  the  regions  of  stability  with  respect  to  a step 
input  can  be  compared  in  the  a - A plane  on  Fig.  IV.4.a,b,c  for  the 
L.T.I.,  the  C.I.  and  the  FORE  case  respectively.  This  proves  then  that 
FORE  is  superior  to  the  C.I.  from  the  stability  point  of  view  (at 
least  for  type  "1"  inputs)  as  already  stated  many  times  in  chapter 
III.  Furthermore  for  the  C.I.  (b = 0)  , the  criterion  suggests  to 

introduce  a zero  at  s = 0 in  G(s)  in  order  to  stabilize  the  nonlinear" 


UNSTABLE  H 


STABLE 

REGION 


A 

a 
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> 

A 
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i\  I A 
,?/  0 
A 
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feedback  system  of  Fig.  IV. 3. a.  Another  alternative  is  to  introduce 

two  poles  at  s = 0 in  the  ooen  loop  transfer  function  L (s) 

eq 

Those  results  can  be  physically  easily  understood.  Indeed, consider 
a L.T.I.  feedback  system  containing  a linear  integrator  (L.I.)  in  its 
loop,  whose  input  is  x and  output  y.  When  a type  "0"  input  R 
is  applied  to  such  a system,  as  t -*•”  , R-+0  , X-»  0 and  y ■*  0 
If  the  L.I.  is  replaced  by  a C.I.,  one  notes  that  the  steady  state 
(x  = 0 , y = 0)  is  compatible  and  thus  no  limit  cycle  can  sustain.  If 
now  R is  type  "1"  , then  as  t , R -*■  k , X-*-0  and  y-*-k  for 

r y 

the  L.T.I.  system.  Unfortunately  the  steady  state  (X=0  , y * k^  ) is 
impossible  to  sustain  in  the  nonlinear  mode  with  a C.I.  in  place  of  the 
L.I.  and  therefore  a limit  cycle  occurs.  (Note  that  this  mode  is  comp- 
atible with  FORE,  thus  explaining  the  difference  frcm  the  stability 
point  of  view.)  To  place  some  derivative  in  G (fig.  IV. 3. a)  as  is 
suggested  by  the  criterion  implies  to  place  another  integration  in  the 
loop  (because  of  the  factor  1/ G ) and  therefore  the  state  (X=  0 , 
y = 0)  for  the  C.I.  becomes  a compatible  one  and  a limit  cycle  is 
therefore  avoided. 

IV. 4. c.  Type  "2"  inputs. 

In  this  case,  e = 2 when  b = 0 and  e = l when  b ^0  . This 
suggests  that  a limit  cycle  only  sustains  when  bWO  . if  b = 0,  amplitude 
of  oscillations  would  increase  with  time.  For  example,  if  G =1  , 


L (s) 

on 


A 

s (s+2a) 


, R(s)  = — , then  for  b^O,  we  have: 
s 


r x x i 

' s 2 

s +2cs+A  J 


7 7 

[ s +2os+A 


we  have: 
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If  o2  -A<0  and  using  = A-o2  , 


AMs)  =■ 

(z- 

and : 

X(z)  = 

7 [32Z 
L s 

so,  as 

z 1 

X(z)  -*■ 

¥ «*• 

..  -oT  o -oT  . -oT,  -oT  _ a . 

z(l-e  cosco„T e sincj-T)  + e (e  -cosid-T  + — sinunT) 

0 cd  0 0 0Dq  0 


a -oT 
e s 


2 o 4o 

T 5 + “ « 1 

A A 


-oT 


s +2os+A 


« -<7T  ■ 

Az  e sinojgT 


+ «(z) 


2 - -OT  _ , -20T 

z -2ze  coscdqT  + e 


m „ 2 ,2 

-OT  2o  -A 

^ — r sinu.T) 

2o(a)q o 


-OT 


-OT  . 


(z-1)  (1-e  cosuIqT  - e sinai  QT) 


2 -1 

and  with  o = , , A = (i>„  and  4 = cos  C we  know  that  a limit 

N N 


cycle  exists  if: 


3 T such  that  lim  (z-l)X(z) 
z=l 


/l-e'^V  sin<u0T  + 2t 


-f1- 

“jlL 


sin2^ 


. _ sin(iDoT  x 4) 

l-e’^  "Vin*-'" 


Let  a = 0 (FORE)  then  we  have  to  satisfy: 

sin2  $ = e ^UNTsin(cD0T  + 2 <M  and  3 T^O  , if  -1  < C < .19  , 

for  which  case  a limit  cycle  exists.  Region  of  stability  with  respect 


Figure  IV. 5.  Stability  regions  of  FORE  with  ramp  inputs. 


to  a ramp  input  are  then  shown  on  Fig.  IV. 5 for  the  FORE.  In  orde: 


to  overcome  the  unstability  problem  noted  here,  the  criterion  suggests 


like  G 


seems  then  to  be  reasonable  in  order  to  avoid  possible 


limit  cycles  in  the  nonlinear  ramp  input  response  and  more  generally 


inputs 


s +2cs+A 


t s +2as+A  J 


We  should  underline  here  that  theorem  2 only  predicts  limit  cycle 


with  one  reset/cycle.  This  type  of  cycle  usually  only  occurs  with 


2 are  not  satisfied,  as  it  does  usually  for  high  order  systems,  one 


should  investigate  possible  limit  cycles  with  2 resets/cycle  which  is 


the  most  commonly  found  experimentally 


If  a limit  cycle  of  period  T with  two  resets/cycle  sustain,  it 


means  that  the  set  of  reset-instants  is  composed  of  two  susbsets 


deduced  from  9. 


Fig.  IV. 2. b can  be  considered  as  the  sum  of  two  sequences 


each  of  them  of  period  T 


£ y^Cfn+DT)  6(t-(n+X)T)  *y 


I y. (nT) z nT  + l Y. ( (n+X)T) z 
n=0  n=0 


Figure  IV. 6.  Equivalent  representation  of  FORE  when  a limit  cycle 
with  two  resets  per  cycle  occurs. 


Therefore  Fig.  IV. 2. b is  equivalent  to  Fig.  IV. 6,  with 


1_e-(s+b)0K 

B (s)  = (recall  4.7)  where  o„  = t„  - t„  , . When  such 

e S+D  K K K-l 


a limit  cycle  sustains,  o = t - t . has  two  distinct  limits 

K K K- 1 


a*  = (n+A)T  - nT  = AT  (4.22)  and  a2  = (n+l)T  - (n+A) T = (l-A)T  (4.23) 
K K 


Therefore  (recall  4.7  and  the  definition  of  y*(t)),  B^fs)  = B^ts)  = 


l-e'(s+b)aK 


s+b 


(4.24)  when  acting  on  the  sequence  yj  1 ( t ) while 
2 


1-e‘(,+b)°K 

B (s)  = 3_  (s)  = (4.25)  when  acting  on  y?  _(t).  There- 

e 2 s+b  x.  ,2 


fore  the  nonlinear  feedback  system  of  Fig.  IV. 3. a becomes  then  equiv- 
alent to  the  linear  sampled  data  feedback  system  shown  on  Fig.  IV. 7 


from  the  input  output  point  of  view. 


•A  TS 

Let  then  Z [K(s)  e ] = K(z  , m)  with  m + A =1  denote  the  so-called 

modified  Z-transform  (see  [Ll] , [Jl] , [S3]).  (4.26)  becomes: 

_XT<? 

x(s)  = (s+b)P  + e*(s)B1(s)A(s)  (s+b)  + e*(s)B2(s)A(s)  (s+b)e  (4.27) 


and: 

X(d  -ATS 

£,  <s)  = = p + E?(s)B,  (s)A(s)  + e*(s)B_(s)A(s)e  (4.28a) 

| 1 s+b  11  22 


XTq  XT^  XT6? 

e (s)  = e (s)e  = Pe  + e*(s)B,  (s)A(s)e  + e*  (s)  B (s)  A(s)  (4.28b) 
2 1 11  2 2 

Therefore  after  taking  the  Z-transform: 

— XT9 

E*  = P*  + E*  (B^A)  * + (B2Ae  )*  E* 
e + = (PeXTS)*+  E*  (B]AeXTS)*  + (B2A)* 


\ 


or  equivalently: 

P*(1-AB2)»  + (AB2e“XTS)*  (PeXTS)* 

E*  = 

1 A 

XT9  XT0? 

(Pe  )*  (l-AB  )*  + (AB  e )*  P* 
e. 

XT^  — XT6? 

with  A = (1-AB2>*  (1-AB^*  - (ABje  )*(AB2e  )* 

Noting  that : 


(4.29) 

(4.30) 

(4.31) 


(s+b)X  T 


A*(Z) 


A* ( z,m) 


(s+b) (l-X)T 


A*  (z ) 


A*(z,X) 


(s+b) XT 


(l-X)ST 


AB*(z,X) 


A*(z.X) 


A*  (z) 


s+b 

■ ( s+b)mT 


AB* (z,m) 


A*(z,m)  -z  e 


A*  (z) 


P(z) (1-A  (z)  + e 


( 1 - A (z)+e  ^bTA  (z,m) ) ( 1-A  (z) +e 


zP(z,X) (1-A  (z)+e 


A (z  jn))  + P (z)  (zA  (z,X) -A  (z)e 


Taking  the  .’-transform  of  (4.27)  and  using 


XT  (s+b) 


( B A ( s+b) ) 


•mT(s+b) 


(BA (s+b) ) + = Z 


JT(z>  = ( (s+b)  P)  + + (A(z)-e 


A(z))c  (z) 


and  by  taking  the  z-transform  (I(s)e  ls)*  (4.34b) 


x(z,X)  = ((s+b)P)»(z,X)+(A(z,X)-e~  A(z))e,  (z)  + (A(z)-e 


If  a limit  cycle  with  two  resets  per  cycle  sustains,  then  the  equence 


of  reset  values  {jr 


/ r(c)e  k dc ) defined  in  chapter  III 


is  composed  of  two  subsequences 


(n+X)T 


r(Oe-b(,n+X)T-9d, 


an'’ : 


*4  * 


X 2 ( n+1) 
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(n+1) 

(n+X)T 


z(;)e-b((n+l)T-Odc  , 


with  and  independant  of  K 

Therefore,  by  anoloqy  with  (4.10)  to  (4.13),  at  the  steady  state: 


e*(z)  = 


V 


!z-l) (z-e  b‘) 


(4.35)  E*(z)  = 


V 


(z-1) (z-e 


-bT) 


(4.36) 


Conversely,  if  it  exists  p,  , u2  , T such  that,  for  the  system  of 

Fig.  IV. 7,  (4.35)  and  (4.36)  are  satisfied,  it  implies  that  as  t ■* » , 

the  reset  values  X*  have  two  distinct  limits  U , U_  and  therefore  a 
K ± ^ 

limit  cycle  of  period  T with  two  resets  per  cycle  sustains. 

Therefore  we  have  proven: 

Theorem  3 : A nonlinear  feedback  system  which  contains  the  r.on- 

1 * 

linear  element  ( — — ) , whose  output  y is  reset  to  zero  whenever 

S+D  Ct 

the  input  X crosses  a , for  a given  a , sustains  a limit  cycle  of 
period  T with  two  resets  per  cycle,  if  and  only  if,  there  exists 

and  X £]0,1[  such  that  in  the 

neighbourhood  of  t = 00  we  have  : 


finite  nonzero  T , , V2  ^ 


lim  (z-1) (z-e  bT)e^(z)  = 
z=l 

lim  (z-1) (z-e  bT) e 2 (z)  = v2 
z=l 


(4.37a) 

(4.37b) 


1 im  (z-1)  X(z)  = a 
z=l 

lim  (*-l)r(t,  X)  « 


(4.38a) 

(4.38b) 


ir»  more  complicated  than  those  of 
• here  too,  as  shown  below,  to 

• to  whe-  a limit  cycle  occurs, 
r * »*  »od  4.3B)  if  one  is 


1 

l 

! 

i 

i 

i 

I 

i 


j 

i 
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A*  (z) 

A 2 

= , -bT  (4.40) 

(z-e  ) 

Two  cases  are  now  investigated. 

a)  . b = 0 and,  as  it  is  implicitely  assumed  that  L (s)  has  only 

eq 

one  pole  at  s = 0 , Z(s)  has  no  pole  at  s = 0 , and  therefore  all 

Z*  have  no  pole  at  z = l . As  lim  (e  m^TI*(z,X)  - Z*(z))  ? 0 (4.41) 

z=l 

lim  ( zZ*  (z,m)  - e_mbTZ*(  z) ) t 0 (4.42)  , lim  (e_AbTZ*  (z,m)  -Z*  (z)  ) *0  (4.43) 
z=l  z«l 

lim  (zZ*(z,A)  - e‘XbTZ*(z)  ) *0  (4.44)  and  lim  (z_1Z*  (z)  2 - Z * (z  ,m)  Z*  (z , A ) ) 
z=l  z=l 

/ 0 (4.45)  , we  conclude  that  the  number  of  pdes  of  AJ  , AJ  at  z = 1 

are  exactly  equal  to  the  number  of  poles  at  z = l of  P(z)  (which  has 

the  same  poles  as  P(z,x>) 

b)  . b/0  and  then  £(s)  has  one  pole  at  s = 0 , whose  residue  is 

1/b  . The  inequalities  (4.41  to  4.45)  are  here  too  satisfied.  Further- 
more, all  terms  Z*  in  the  numeration  of  both  e*  and  e*  have  one 
pole  at  z = 1 which  is  cancelled  out  by  the  one  pole  at  z = 1 of 

z_1Z* (z) 2 - Z*(z,m)  Z* (z, A)  . 

v (2) 

Indeed,  one  notes  that  the  latter  has  no  such  term  as  — , with 

(z-1)  2 

-1  2 1 z-1 

Y ( 1 ) ^0  because  y(z)  - z ■ —r  - — - — j and  therefore  y(z)  has 

b b b 

-1  2 

one  zero  at  z-1  , implying  then  that  z Z*(z)  - Z*(z,A)  Z*(z,m) 
has  only  one  pole  at  z=l  . Therefore  we  can  conlcude  here,  too, 

♦•hat  the  number  of  poles  at  z = l of  AJ  , A*  is  exactly  equal  to  the 
number  of  poles  at  z = 1 of  P(z) 

Therefore  in  all  cases  (b  = 0 , b ^ 0)  condition  I of  Theorem  3 
is  satisfied  if  P(s)  has  one  pole  at  s = 0 , imolying  then  that  a 


I 


limit  cycle  with  two  resets  per  cycle  may  sustain.  The  reader  has 
certainly  noted  that  P(s)  is  exactly  the  same  quantity  than  the 
numerator  of  A*(s)  in  (4.17)  [Theorem  2].  Therefore  the  qualitative 
discussion  made  in  section  IV. 4 for  Theorem  2 is  easily  extended  here 
and  it  is  seen  that  the  occurence  of  a limit  cycle  is  entirely  deter- 
mined by  the  number  of  integrations  (type)  of  the  system- input,  while 
the  limit  cycle  itself  is  mainly  characterized  by  the  loop  transfer 

function  L (s)  . If  G has  no  zeros  at  s = 0 , we  can  conclude  then, 

eq 

that  a nonlinear  feedback  system  containing  a C.I.  sustains  in  general, 
a limit  cycle  with  "type  1"  inputs  but  such  limit  cycle  does  not  occur 
when  the  C.I.  is  replaced  by  FORE.  However  a nonlinear  feedback  system 
containing  FORE  sustains  a limit  cycle,  in  general,  with  "type  2" 
inputs . 

It  remains  to  show  that  condition  II  of  Theorem  3 is  a consequence 
of  condition  I.  This  is  now  established. 

Replacing  (4.39),  (4.40)  in  (4.34a)  we  get: 


Jrtz)  = ( (s+b)  P)  * + A ( z) 


r“inbT 

P(z)(e  S»(z,X)  - E«(z) ) + P(z,X) (zE»(z,m) 


_e  mbTr * ( 2 ) ) -z  ^e  mbT(zP(z,l) (e~XbTE«(z,m)  -E*(z))  + zP (z) (£♦ (z, 1 ) 

-1  -XbT  * i A 

-z  e I ( z ) ) ) 


+A(z,m) 


f zP ( z , X ) [e'XbTZ*(z,m)  - I«(z) 1 + P(z) [z£*(z,X)  - e"XbT£«(z) ] 


-e"*bT(P(z)  [e-mbTj:*(z,l)  - l*(z)  ] + P(z,X)  [z£*  (z,m)  - e-mbTE*  (z)  ] ) 


or  X (z)  = ( (s+b)P) * + 


-1  -bT 

A(z)  (~P(z)E*(z)z  x + P(z,l)E«(z,m))(z-e  oi) 

- 1 -HT 

(z  E*  (z)  - £*(z,m)E*(z,X))(z-e  D1) 


+ A(z,m) 


(P(z)E*(z,X)  - P(z,X)E*(z))(z-e~  ) 

-1  2 -bT. 

(z  Z*(z)  - E* (z,m) I* (z,X) )(z-e  ) 


(4.46) 


Equivalently  from  (4.34b)  we  have: 
X (z,X)  = ( (s+b)P)*(z,X)  + 


-1  A (z)  (P  (z)  E (z ,1 ) - P(z,X)E(z))  + A(z,X)  (P(z,X)E(z,m)  -P(z)E(z)) 


z~lT(z)2  - E(s,m)E(z,X) 


Defining  n = lim  (sP(s))  we  have  then: 
s=0 


(4.47) 


C(T,X)  = lim  x(z)  (z-1)  =n 


A*  (1)  (E*  (1  ,m)  - E«(l))  -t-A*(l,m)  (E*  (!»!)- E*  (ID 


E* (1)  - E* (1 ,m) E* (1,X ) 


(4.48a) 


and  limxtz, X)  (z-1)  =n 
z=l 


A*  (1)  (E*  (ljX)  - E*  (1) ) + A*  (1 , X)  (E*  (1  ,m)  - E»(l)) 


E* (1)  - E* (1 ,m) E* (1 ,X) 


(4.48b) 


=C(T,1-X)  (by  definition  of  ( ) 


recalling  that  condition  X is  satisfied  if  P(s)  has  one  pole  at  s=0. 


we  conclude  (in  that  case  n^O  ) that  condition  II  of  Theorem  3 is 

a consequence  of  condition  I. 

IV. 6.  Example  of  application  of  Theorem  3. 

2 1 
Let  L (s)  = q and  G = 1 . Consider  R(s)  * — . 

eq  s(s+l)2 

2 

Then  A (s)  = - and  thus  it  can  be  seen  that  the  L.T. I. system 

(s+2) (s  +1) 

obtained  is  unstable  (oscillatory) . The  corresponding  system  response  is 

plotted  on  Fig.  IV. 9. 

Suppose  that  b = 0 , a = 0 (C.I.)  then: 


P(s)  = 


s(s+2)  (s  +1) 


so  n = lim  sP  (s)  = — 
s*0 


A (s)  = — = — + E (s)  with  E (s) 

1 s(s+2) (s  +1) 


.2  .4 (2s+l) 

s+2  " 2 


and  thus: 


£ (z,X) 


2z(z-cosT)+ZsinT 

2 

z -2zcosT  + 1 


2(z  cosXT-cosmT) +zsinlT+sinmT 

2 

z - 2zcosT  + 1 


S' 
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Besides : 


\z-e 


A(z,m) 


z(z-cosT)  - 2zsinT 


z - 2zcosT 


zsinT  \ 
+ 1 ) 


(zcosmT  - cosXT)  - 
2 

z - 2zcosT 


2 (sinmT  + sinXT)  \ 
+ 1 ' 


C(T,X)  0f  (4.48a)  and  (4.48b)  is  thei  plotted  on  Fig.  IV.8.a  cor 
different  values  of  T and  X . (The  reader  should  remember  that 
m * 1-X  so  relations  (4.48a)  at  X = .2  say  corresponds  also  to 
relation  (4.48b)  at  hi  * .8  ) . 

Condition  II  of  Theorem  3 implies  that  we  should  satisfy: 

3 T,X  such  that  C(T,X)  ™ C(T»1-X)  = 0 (A. 49) 

From  Fig.  IV .8  .a  (4.49)  is  satisfied  for:  3<T<5  and 
.1<X«.2  , and  therefore  this  region  is  magnified  on  Fig.  IV. 8. b. 


.04  4.05 


Figure  iv.8.  Plot  of  versus  T,  for  the  characterization  of  a 

limit  cycle  with  two  resets  per  cycle. 
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Figure  IV. 9.  Comparison  of  the  three  system  step  responses. 

Thus  Theorem  3 predicts  a limit  cycle  of  period  T such  that 

4.02  <T<  4.03  with  \ such  that  . 176<X<  .177  . This  is  then 

confirmed  by  an  analog  simulation  of  the  nonlinear  system  includirg  a 

C.X.  as  shown  on  rig.  IV.  9. 

If  now  b^O  ,a«0  (FORE)  it  was  already  mentioned  that  no 
limit  cycle  can  sustain  because  P(s)  does  no  longer  have  a pole  at 
s = 0 implying  then  that  A*  and  A*  do  not  have  any  pole  at  ~ » 1 . 
This  is  conf irmed  by  an  anolog  simulation  of  the  nonlinear  system 
including  FORE  as  shown  on  Fig.  IV.  9,  with  b ■ 1 . 

One  should  be  aware  that  the  Dual  Input  Describing  runctions  fc2 ’ 
give,  here  too,  wrong  results,  except  for  the  period  of  the  limit  cycle 
sustaining  in  the  N’L  system  including  the  C.I.  Indeed,  it  can  be 
shown  that  when  b=0  , the  D.I.D.F.  predicts  a limit  cycle  of 

frequency  w w 1.5  rd/s  (T«4s)  with  ratio  .7  < A/B  < .75  (instead  of 
5)  where  the  assumed  input  is  X=  B + Asinui^t 

For  b/0  , the  D.I.D.F.  predicts  either  a limit  cycle  or 

instability  (compare  with  above)  and  for  instance  with  b = 1 , a 

limit  cycle  is  predicted,  characterized  by  1.  <uj0<1.1  and  a ratio 
.8  < A/B  < .9  . 


r*i 

i 
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V.l.  Introduction. 

The  C.I.  has  been  generalized  into  TORE.  It  i*  therefore  natural 
to  think  about  possible  generalizations  of  FORE.  One  way  is  by 
resetting  the  output  y whenever  the  input  Z • a*  0 . However, 

nothing  is  gained  by  doing  this  when  both  positive  and  negative  '-orevr 

| I 

inputs  are  to  be  applied  to  the  system.  Resetting  for  z * 1 * is 

conceivable  but  the  performances  for  two  input  signals  Z(t)  and 
k-Z  (t)  , with  k real  * 0 , would  no  longer  be  proportional  with 

ratio  k and  this  is  certainly  a big  weakness  if  the  objective  1st 
guarantee  T.D.S.  of  a linear  type. 


Another  way  to  generalize,  is  to  extend  FORE  to  LLRE  A** ' • lag 

1 


s*  a 


(third-order  ...)  ■ ( 


SORE  (second-order  ...)  * < 

1 


TOPE 


-) 


(s*b)  ( s*al 
etc...,  and  in  a very  general 


is*  a)  ( s*b)  (s*c) 

manner  to  g(s)*  , where  g(s)  is  any  rational  transfer  function 

whose  output  y is  reset  whenever  the  input  z is  zero.  (Fig.  V.l) 
V.2.  First  equivalent  representation  of  g*  . 


Let  9 ■ (tR  : x( t^)  • 0}  be  the  set  of  reset  instants.  Then  ■ 


” tK 

-:)d(; 

- 1 ;z(C)  qlt-cldc 

u ( t— 

K=1  K-l 

1 V 

g 

■ 

g* 

_JReset  when  x = 0 


Fiaure  V.l.  General  reset  element  (Equivalent  notation) 


84 


Suppose  we  restrict  G to  elements  with  m real  poles,  so 
m t 

?(t)  - (t)  ♦ l A e'  1 

i-1 


I . lUI*  • k 

y(t)  - /y(C)q (t-OdC  - £ *i  I /Z( C ) e'^1 ' dC  / e"ai  ' ' T,«(t-tK)d 

0 i-1  K-l  tK-1  0 

orif  1J,  J / e"1i,tK'r’z(i:)<5C  with  tR  € 0 

Vi 

y (t)  - /Z( C)*(t-C)4C-J  [ *,  Ij  / e'ai(t'T,S(T-tK><lT 

0 K-l i-1  ’ 0 

-kj  to  the  equivalent  representation  ehown  on  Fig.  V.2. 

jf  q*  is  imbedded  in  the  feedback  system  of  Fig.  V.J.a,  it 
li,  be  replaced  by  its  equivalent  representation  of  Fig.  V.2.  There- 
fore, by  application  of  the  superposition  theorem  to  all  the  signals, 
he  output  c(t)  can  be  computed  by  using  the  equivalent  represen-ation 


IXK.l*(t-V 


J1XK.ia,t-tK> 

K-l 


l x*  t-t  > 

K-l  K'm  K 


Figure  V.2.  Equivalent  representation  of  g* 


***■ 
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of  Fig.  V.3.b.  Using  such  a representation  it  is  possible  to  derive 
charts  similar  to  Fig. III. 10  for  any  specific  g . If  one  is  then 


from  the  command 


m such  that 


input  point  of  view,  the  design  philosophy  used  with  FORE  can  then  be 
applied  here  too.  This  was  used  for  instance  with  LLRE,  but  no 


improvement  over  FORE  could  be  noticed 


such  that 


then 


we  have 


Using 


and  therefore,  by  anology  with  the  previous  chapter,  ( section  IV.  3)  the 


system  of  Fig.  V.l  is  equivalent  to  that  of  Fig.  V.4 

A 1_e"(S+ai,CK 

a 3 t -t.  , and  B.  (s)  = — • a 


Figure  V.3.b.  Equivalent  representation  of  two 
degree  of  freedom  structure. 


By  anology  with  Theorem  2 we  can  state 


If  in  the  neighborhood  of 


such  that 


V.5,  some  finite  nonzero  numbers  T 


then  a limit  cycle  with  one  reset/cycle  exists 


Application:  'LL  FORE 


value  of  T than  in  the  existence  of  a limit  cycle 


This  will  now  be  illustrated  by  considering  the  simple  example 


L(s)  is  assumed  to  have  only  one  pole  at  the  origin,  and  L is  such 
L (s) 

that  L(s)  ^ . So  L/l+L  has  no  pole  at  s = 0 while  1/1+L 


has  one  zero  at  s = 0 


then 


at  the  origin.  As  there  is  no  term  in  the  denominator  of  both  e*  and 


containing  the  factor  1/s  , there  is  no  such  factor  as  1/z-l 


in  the  denominators  of  both  e*  and  e* 


According  to  the  criterion 
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there  exists  a limit  cycle  if  e*  and  £*  contain  the  factor  1/z-l  , 

Therefore,  this  can  only  occur  here  if  the  numerators  of  both  t*  and 

e*  contain  such  a factor,  i.e.  for  "type  2"  inputs  R , and  especially 

with  ramp  inputs.  Indeed,  R/l+L  has  then  a pole  at  the  origin  and  the 

A. 

residue  at  z = l of  the  number  e*  is  — - — rrr-  (1*1,2)  which  is 

i a . L (0) 

l 1 

finite  non  zero. 

The  stability  criterion  then  implicitely  suggests  the  insertion 


"type  2"  inputs.  Indeed  we  have  then: 


of  a 

zero 

that 

V; 

"type 

2" 

B191 

L 

g 

1+L 

= (i 

-1 
-Z  i 

B2g2 

L 

g 

1+L 

Vl 

L 

g 

1+L 

B g. 
l’2 

L 

g 

1+L 

1 ' 2 1 ' 2 


-1  -otT  A1  S+a2 
(1"Z  6 > — 


Q(s)  ) , with  Q (s)  £ TTfr 


\ (■  <*2Q(0) 

a1+a2  L ~ 


i = , = 1+L(S) 

( s ) 1 (after  fractional  expansion). 


Al+A2  L 1 S 

fa  2I2J 

Al+A2  L 2 S 


-1  -oXT,  __1 


V*2  !'‘l2Ti  * V”  ] 


The  denominator  of  both  e*  and  e*  is  then: 


, _ A,  a,Q(0) 

. ,,  -1  -anT,  1 2 

i-U-z  e 1 ) - — 

A1  A2  1-z 


n -1  -a,T.  A2  al®(°) 

a’z  e ’ rzr  ~^r 

1 2 1-2 


>(0)  (l-z~1e~C‘^T)  (l-z~1e~a2T) 


(a1+a2> (1-z"1) 


(El(2)A2ai 


+ j:2(z)Aia2  “ E1^A2°1 


- S2(2)A1 


(l-z'1e"r'lT)  ( l-z-1e"  °^T)  Aj  A2  Q(0) 


(A1*A21  <1_Z  *'2 


(aia2  ” * P(2) 


■ re  I • ( I ) / "» 

w ; * b a residue 


So  the  denominator  contains  now  the  factor  1/1-z 


-X 


am  (A^e-«xT  + a > .-«2T 


a1+a2 


12 


•lA2e 


) ^ 0 because  by 


assumption  A^a^  + u^A2  = 0 


It  is  easily  seen  that  the  numerators  of  both  e*  and  e*  have 


-1  2 

such  term  as  l/(l-z  ) 


Therefore,  the  numerators  of  both  e* 


ai.d  e*  have  at  most,  a term  in  1/1-z 


-1 


and  therefore  if 


Ai'ij  + ai^i  = 0 b°th  e*  and  e*  do  not  contain  such  a factor  as 


1/1-z 


-1 


preventing  then  the  system  from  limit  cycling  in  presence  of 


"type  2"  inputs,  and  in  particular  with  ramp-inputs.  Relation  (5.7)  is 


for  instance  satisfied  by  choosing  A 


1 , . 2 

and  A_  = 

a, -a.  2 a -a 


21  "2  ”1 

that  is  referred  to 


lending  to  the  element  G(s)*  = (- — -\ 

V S+ai  S+a2  / 

as  a Lead  - Lag  - and  - First-Order  - Reset  - Element  ( LLFORE ) 


/ s i \* 

V.5.  Synthesis  procedure  with  LLFORE  s+b7 

Using  the  results  of  section  V.2  (Fig.  V.3.a  and  b)  the  nonlinear 


feedback  system  of  Fig.  V.6.a  is  equivalent  from  the  output  signals 
point  of  view  to  the  linear  system  of  Fig.  V.6.b  with: 


X*  A 

K,1  - 


/ - *•. 

bK-l 


FORE 


X*K  2 = f e'a(t"C)x(OdS 
fcK-l 
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If  ;<  is  so  chosen  as  to  satisfy  a <<  Bandwidth  of  the  open  loop 
transfer  function  for  all  possible  plant  parameter,  then  obviously 
a <<  b , therefore  we  can  consider  that: 


b s+a  , , a s+b  _ , , 

- — 1 and  - — <«  0 and  therefore  LLFORE 

b-a  s b-a  s 

behaves  practically  like  FORE.  This  is  in  accordance  with  experimental 
results:  the  design  of  chapter  III  was  used  and  on  Fig.  V.7  it  is 

shown,  for  the  maximum  plant  gain  factor,  the  non  linear  step  dist- 
urbance response  for  different  values  of  a . As  expected,  for 
a 5 1 rps.  (compared  to  a bandwidth  of  600  rps.)  there  is  hardly 
any  difference  from  the  result  obtained  with  FORE  (see  Fig.  III. 20. b 
(f  1)),  while  as  a increases  up  to  10,000  rps.,  the  nonlinear  system 
response  tends  to  be  like  the  one  obtained  with  a purely  linear  system. 

Therefore  the  design  procedure  derived  with  FORE  needs  only  to  be 
completed  here  by  a suitable  choice  of  a . Since  (Fig.  III. 13)  the 
minimum  bandwidth  of  L^ts)  is  roughly  6 rps.,  one  can  choose  for 
example  a = .01  rps. 


Signal 

Amplitude 


a 

□ 

1 

o 

10 

A 

100 

+ 

1000 

X 

10000 

Figure  V.7.  System  step  disturbance  response  of  the  nonlinear  feedback 

system  at  k * 1000  for  different  values  of  the  parameter  a 


3HK5E 
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As  expected,  one  notes  from  Fig.  V.8.a,b,c  that  no  limit  cycle 
is  sustained  in  presence  of  a ramp  input  for  all  values  of  a , but  it 

is  seen  that  the  transient  gets  bigger  and  bigger  as  a.  decreases, 
which  is  reminiscent  of  the  limit  cycle  that  existed  with  FORE. 

V.6.  Serial  and  multiplicative  combination  of  G's. 

It  is  worthwhile  to  mention  that  g*  can  itself  be  generalized 
by  considering  a system  of  the  type  shown  on  Fig.  V.9.a  or  of  the  type 
shown  on  Fig.  V.9.b,  or  any  system  which  combines  both  types. 

V.6.a.  Serial  Combination  (Fig.  V.9.a) . 


If  x.  = {tR  , y^(tR  ) = 0}  denotes  the  set  of  reset  instants 
1 ^ 


associated  with  the  input  y^  , it  is  then  obvious  that  xQe:  x^c 


..c  x . . If  we  assume  that  each  g.  has  m.  real  poles,  then 

p-1  1 1 


for  i = 1 , 2 , . . .p 
t 


N 


Rl , 11  ^ 

'i««  ’ l 

fcK.  -1  ‘ 

and  mi 

gi  (t)  = h0i  <5  (t)  + l a., 
j-1  13 


«(C-tK  )dC 
i 


The  equivalent  representation  of  Fig.  V.9.a  is  plotted  on  Fig.  v.10. 
V.6.b.  Multiplicative  Combination  (Fig.  9.b) . 

T suc^  that  y0(tK)  = 0}  denotes  the  set  of  reset  instants. 

It  should  be  first  noted  that  if  in  a)  , T = x = . . . = x = x , class 

12  p 


a)  degenerates  then  in  the  class  b)  studied  here. 

Therefore  the  equivalent  representation  of  Fig.V.9.b  is  given  by 


Fig.  V.10  where  the  substitution  t =t  should  be  made  for  all 

. K 
1 


i = 1,2,...  p . 


• rii  tiiinfom i> 
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Signal 

Amplitude 


Signal 

Amplitude 


Signal 

Amplitude 


Figure  V.8. 


0,X  0 .20  O.«0  0.00  0.00  I.X  1.20  ,*c 

TIME 


..fcC  .00  2.X 


(c) 


Nonlinear  system  response  to  ramp  disturbance  for 

different  values  of  the  parameter  a . 

a)  a = 100  rps  , b)  a = 10  rps  , c)  a = 1 rps 
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A systematic  engineering  design  procedure  for  drastically  reduc- 
ing the  'cost  of  feedback'  of  linear  feedback  systems  has  been  pres- 
ented and  illustrated  in  chapter  III.  It  is  clear  that  the  nonlinear 
compensation  (FORE)  introduced  is  justified  in  problems  of  large  para- 
meter uncertainty,  and  that  the  design  procedure  is  a general  one  for 
this  problem  class.  The  nonlinear  design  permits  the  attainment  of  the 
same  command  and  disturbance  performance  tolerances  as  a linear  design, 
but  with  significantly  smaller  loop  transmission  bandwidth.  An  import- 
ant feature  of  the  design  procedure  is  that  it  permits  design  to 
quantitative  specifications,  a property  generally  lacking  in  present 
nonlinear  feedback  synthesis  techniques,  for  systems  with  significant 
parameter  uncertainty. 

It  is  important  to  note  the  inherent  assumption  in  that  chapter, 
that  the  primary  design  problem  is  that  of  satisfying  the  response 
tolerances  to  command  inputs.  However,  realistically,  consideration  is 
also  given  to  disturbance  inputs  in  the  form  of  steps  at  the  output, 
with  assigned  restrictions  on  the  damping  of  the  resulting  output. 

Since  the  disturbance  response  is  nonlinear  (in  contrast  to  the  command 
response  which  is  essentially  linear) , one  cannot  in  general  guarantee 
acceptable  response  to  all  possible  disturbance  inputs. 


Therefore  emphasis  was  then  placed  (chapter  IV)  on  the  stability 
problem  for  this  class  of  nonlinear  systems.  Suff  conditions  were 

derived  for  B.I.B.O.  stability  and  it  was  conjectured  that  the  nonlinear 


Furthermore,  necessary  and  sufficient  conditions  were  derived  for 
the  existence  of  a limit  cycle. 

The  results  of  chapters  III  & IV  were  then  used  in  chapter  V to 
derive  the  new  element  LLFORE  which  can  be  considered  a generalization  of 
FORE.  Indeed,  the  latter  improves  the  performances  of  FORE  with  respect 
to  ramp  disturbances  (and  more  generally  with  respect  to  "type  2"  dist- 
urbance inputs) , without  harming  any  of  the  quantitative  benefits  obtained 
with  FORE.  The  design  procedure  derived  in  chapter  III  is  therefore 
easily  extended  to  this  nonlinear  element. 

The  philosophy  and  motivation  prompting  the  nonlinear  compensation, 
is  based  on  linear  frequency  concepts,  coupled  with  linear  feedback 
design  techniques  for  guaranteeing  performance  tolerances  despite  large 
parameter  uncertainty.  'This  philosophy  has  thus  proven  itself  as  at 
least  one  approach  worthy  of  pursuing.  One  might  search  for  other  non- 
linear elements  with  even  greater  phase  advantages  than  FORE,  over  linear 
elements  with  the  same  magnitude  characteristic.  As  an  example  the  non- 
linear element  such  that:  y =*  { * /x(c)d£  if  x i < 0 is  a "°°  Phase 
lag  integrator"  from  the  describing  function  point  of  view.  With  FORE,  the 
equivalent  linear  phase  lag  useable  was  shown  to  be  almost  180°.  Non- 
linear elements  can  undoubtedly  be  found  which  in  a describing  function 
charaterization,  would  permit  even  greater  phase  lag. However,  as  with  the  C° 
phase  lag  integrator,  the  more  difficult  challenge  is  to  find  a charac- 
terization of  the  nonlinear  element  useable  for  the  useful  system  control 
signals,  and  so  permitting  design  to  quantitative  specifications. 

It  will  certainly  be  worthwhile  in  the  future  to  extend  this 
research  to  multivariable  systems,  once  a rigorous  synthesis  procedure 
is  available  for  L.T.I.  multivariable  systems  with  large  plant  ignorance. 


CHAPTER  VI ■ 


LINEAR  TIME  VARYING  COMPENSATION  OF  FEEDBACK 
SYSTEMS  WITH 


VI . 1 Introduction 

It  was  noted  (section  II. 3)  that  the  typical  set  of  command  inputs 

to  a system  can  often  be  imbedded  in  the  set  of  non-stationary  processes. 

Therefore,  it  is  expected  that  L.T.V.  compensation  in  a feedback  system 

will  result  in  better  performances  than  L.T.I.  compensation,  with  respect 

to  sensor  noise  effects.  It  is  implied  that  the  instant  of  input  (R) 

application  is  known.  Given  (F,G)  (Fig.  VI. 1),  we  get  different 

output  responses  c^(t)  for  the  same  command  input,  due  to  different 

plant  parameter  combinations  P^  € IP  , as  shown  in  Fig.  VI. 2.  At  each 

instant  of  time  tQ  , the  maximum  spread  is  then  characterized  by 

Ac(tQ)  (Fig.  VI. 2).  The  L.T.V.  networks  F,G  can  be  associated  [S5] 

with  L.T.V.  operators  f,g  , i.e.  (Fig.  VI. 1): 
t t 

v(t)  = / f(t,C)  r (C)d£  and  x(t)  = J g(t,C)  e (C)dC  can  be  written  in 

0 0 

the  symbolic  notations: 

v = f*r  and  x = g*e 

Therefore,  the  effect  of  the  noise  n(t)  at  the  plant  input  x (assum- 
ing R = 0 ) is  given  by:  (Fig.  VI.  1) 


Figure  VI. 1.  System  Block  Diagram 


system 

output 


Figure  VI ■ 2.  System  Output  Response  for  Different  Plant  Parameters 


x = -g*  n -g*p*x 
n n 

x^  * - (1+g*  p)  1*g*n  if  (l+g*p)  1 exists. 

x = 0*n  (6.1. a) 

n 

t 

whore  0 is  a L.T.V.  operator,  i.e.  x (t)  = / 9 (t,£)  n (G)d5 

n 0 

The  mean  square  value  of  the  noise  at  the  plant  input  defined  as 
2 2 

c . (t)  = <x  (t)  > , where  the  bracket  sign  stands  for  the  ensemble- 

P . I . n 

average,  is  therefore: 


y f u 

°P.I.  (t)  = < / 6(t,C1)  n (r.^d!^  / 9(t,C2)  n (C^dl^  > 

t t 

= / / 9(t,C  ) 0 (t,C,)  < n(C  ) n (;  ) > dc.d?  . 

00  1 


By  definition,  the  autocorrelation  function  of  the  noise,  is 


= < r-Uj)  n > ' so 


t 

r r 


1 J eft,;,)  9 (Z  ,r  )dc.di; 


(t) 


( 6 . 1 . b ) 
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Using  arguments  analogous  to  those  used  for  L.T.I.  systems,  it  can  be 

2 

shown  that  the  spread  Ac( t)  is  a conflicting  factor  with  cp  (t) 
namely,  the  bigger  Ac  the  smaller  o2  , and  vice  versa.  The 
challenge  then  is  to  solve: 

Vt  Min  { Ac2  (t)  + W2(t)a2  (t)>  (6.2) 

(f,g)  L.T.V. 

with  W(t)  some  given  weighting  function. 

Because  of  the  complexity  of  P in  general,  solving  (6.2)  over 
the  range  of  plant-uncertainty  is  a difficult  task  for  which  no 
techniques  exist  at  present.  Let  c^(t)  be  the  system  response  to 
input  r when  P = P ^ . It  is  tempting  in  order  to  have  a solvable 

problem,  to  pretend  that  the  maximum  range  Ac(t)  in  (6.2)  consists 
of  the  differences  (c^  - c 2>  and  this  for  all  t , in  the  system 
response  due  to  2 'extreme'  plant  conditions  P^(s)  and  P2<s)  ' to 

the  same  input  r , i . e . , 

vt  , Ac(t;2  = (Cl(t)  - c2(t))2  (6.3)  . 

In  practice  the  extreme  points  of  the  spread  Ac(t)  in  Fig.  VI. 2 do 
not  necessarily  correspond  to  the  outputs  (Vt)  for  any  2 plant 
conditions  chosen  as  'extreme'.  The  weighting  function  W(t)  in  (6.2) 
is  then  an  extra  degree  of  freedom,  to  help  contend  with  the  fact  that 
the  upper  and  lower  bounds  (Fig.  VI. 2)  of  the  system  responses  for 
P.  € , do  not,  in  general,  correspond  for  all  t , to  the  two 

'extreme'  plants  P^s)  , P^Ls)  , i.e.  (6.3)  usually  does  not  hold  in 
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a realistic  problem. 
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Using  (6.4.d)  and  (6.9)  in  (6.8), 
Ac  = d2  - p2*g*Ac  - P2-g"n 

Let 

„ A 

£ = p • g 

2 v2  ’ 

then  (6. 9. a)  is  rewritten  as: 
<1+*2)-Ac  = d2  - Syn 


(6. 9. a) 


(6 . 9 -b) 


or 

Ac  = ( l-*-2. 2 ) ~ 1 * d2  - (l+^2>  -1*  ^2"n  (6.9.c) 

if  ( 1+2- 2 ) has  an  inverse. 

Let  d_  = 0 , then  Ac  = -(l+£  ) "n  which  together  with  (6.4.d) 

z n 2 2 

gives : 

Ax  = g*  C(l+e._) _1* «•_  - 1 ] -n 
n 2 Z 

Ax  = g*  (l+l,)-1'  C - (1+Jtj]*n 
n z z 2 

Ax  = - g*  (l+^_)  1,n 
n z 

which  can  be  rewritten,  (recalling  (6.9.b)),  as: 

Axn  = - p2_1- 2.2-  (l+i2)_1-n  = p2_1*  (l+)l2)'1^2*n 
because  and  (1+5..,)  ^ commute. 

Note  (6.9.b)  that  - p2  1,^2*(l+i2)  1 = -g* (l+p2"g)  1 = - (l+g’p2>  1*g  , 
so  from  (6.1. a),  6 = ” P2  *’*2' * (6.9.d)  when  P = P2 

Hence,  for  the  plant  P = P2  , the  minimization  problem  of  (6.2)  becomes 
(under  the  above  assumption  of  a binary  plant)  equivalent  to: 


Vt  , Min  {Ac(t)2  + w2(t )o2  (t) } (6.10) 

fp  • X • 

-g 

where  Ac(t)  is  the  component  due  only  to  d2  , and  6 of  (6.9.d)  is 

2 

used  in  (6.1.b)  for  op  ^ 


2 


i vat  ion  of  the  Optimum  Filter 


part  of  Fig.  VI.  3 when  an  impulse  is  applied  at  W at  t = T 


effect  of  d is  given  by  letting  n = 0 in  (6.9.c)  , so  in  (6.10) 


ion 


Ac(t)  =<  ^d2(t)  - / h2(t,5)d2(C)dO 


where  the  bracket  indicates  the  ensemble-average  over  the  set  of  command 


From  (2.6) 


inputs  r. (t) 


where  ^ (•)  is  an  N-dimensional  vector 


(6. 12. a)  is  then 


with  k € [k  . ,k  1 a cure  gain  factor  and 
nun  max 


If  the  plant  P = k 


W(t)  = 1 , then  from  (6.9.d,  6. 11. a) ,9 


using  (6.12.c),  (6.1.b)  and  (6.12.d),  (6.10)  becomes: 


Vt  , Min  ^ Y DD  ( t , t ) - 2 / h2(t,OYDD(C,t)dC  + / / h2  (t.C^h.,  (t,^) 


I.  YDD(C1';2)  + 2 YNN(;i'C2>.  d?ldC2  ) 


(6.13) 


By  using  ydd  instead  of  yjd  in  (2.5),  we  conclude  that  the  optimum 
solution  h2  of  (6.13)  is  the  solution  of  the  integral  equation: 


V t f t 


, / h (t,of 


YDD(C'T)  + 72  YNN  (C'TVd?  = YDD(t'T) 
*2 


(6.14) 


We  restrict  ourselves,  without  loss  of  generality  to  white  sensor  noise 
2 2 

of  strength  o , with  resulting  Y._,  (C»t)  = cj  6(?-t)  • Therefore 

N NN  N 

(6.14)  becomes: 


V t5T  ,Uh2(t,T)  t I h2(t,?)J2(C)  l?(T)Td(  = K2(t)l2(T)' 
O „ (6.15) 


a2  a2 

with  U = . (Note  that  p = -y  W(t)  for  W(t)  ?1  .) 

k2  k2 

(6.15)  has  a solution  if: 

h2(t,x)  = <&(t) -i)2(T)Tu(t  - t)  where  is  a N-dimensional 


vector,  and  the  heaviside  unit  function  u(t-r)  = 


duced  to  give  a causal  h2  . (6.15)  becomes: 


Jo  t < r . . 

i is  u 

*■1  t >,  T 


t?T  , u fi-(t)-D,(T)T+  £(t)  J D2(0TD2(C)di;- 3),(t)T  = 


= l>2(t)J)2(T) 


(6.16) 


(6.16)  implies  that: 


G<t)  = 3>,<t)  f ull  + / D,(0TD2(0dc  ] -1 


(6.17) 
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rl  ' '1  if 


»(t)  = ull  + /^>2(C)\,(C)dC  _1 

L 0 J 

hj  (t  ,T  ) = JD2(t)  IT  1 (t)  • 32  (T)Tu(t  - T) 


(6.19) 


(6.18) 


Note  that  it  is  an  N » N matrix. 

From  (6. 11. a)  one  gets  (Appendix  A 3.1)  the  open  loop  impulse  response 


associated  with  P = ?2  , 

*2(t,T)  = 3>2<t)  -rr-1  (x)  • D2  (T)Tu(t  - t)  (6.20) 

Tho  open  loop  impulse  response  fc(t-r)  for  any  P = k 6 is 

(Appendix  A3. 2) 

i(t,T)  = X l (t,r)  (6.21)  with  X = (6. 21. a) 

*2 

and  the  corresponding  closed  loop  impulse  response  (Appendix  A 3.2): 

h(t,T)  = X J)2(t)  A1(t)Al'1(T)TT-1(T)»  ®2(T)Tu(t -T)  (6.22) 


where  A is  a N»  N matrix  solution  of: 
A^  A^-1  + X IT-1  TT  = 0 


(6.23) 


It  is  shown  (Appendix  A3. 3)  that  the  system  response  to  a command 


input  r , when  P = k € is: 


= , + ( imi'i x ^ 

1 ' V ^ V TT(t)  / ) 


(6.24) 


where  c^(t)  is  the  system  response  to  the  same  input  r at  P * k^ 


VI. 2. c Example. 


Let  (r(t);  be  a set  of  step  command  inputs  of  amplitude  a , 

2 

with  average  mean  square  value  < a > = 1 ; P = k€[k  ,k  I is 

min  max 

a real  gain  factor  and  c (t)  is  the  system  response  to  r(t)  =u(t) 


when  P = k^  . From  (6.9),  at  P = k2  , and  r»au(t),  d2  = 


(k2-k,) 


From  (6.12.b),  the  autocorrelation  function  associated  with 
**](  )c  — 

d2(t)  iS  YDD(t'T)  = < “ Cl(t)  -1^  “ °1(T)  > 


k -k  k2-kl  A 

= — Cl(t)  -k-Cl(T)  = D2(t>'D2(T> 

Therefore  N = 1 , and  from  (6.19), 


(6.25) 


Mt)  = 4 + / D K)2dC 
0 


(6.26) 


From  (6.22),  the  closed  loop  impulse  response  at  P = k is: 

TT  (T  ) ^ ~ 1 

h(t,T)  = A D (t)  ■ - D_(T)u(t-T) 

2 *(t)+A  2 


(6.27) 


with  A = kAj  (from  6. 21. a)  , and  the  system  response  to  step  command 
of  amplitude  a is  given  by  (6.24),  i.e. 


c (t)  = 8 Cj(t)  ^ 1 + 


1 * \ 


u * (t) 


(6.28) 


noting  in  (6.26)  that  t (0)  = u . The  effect  of  white  sensor  noise 
at  the  plant  input  is  (see  Appendix  A3. 4) : 


°P.I.(t)  = 


\-l  2 V *(t 


)«  ‘-r 


2 

°2(t)  T 

k2 


En  (tt  (t)/yi) 
it  (t) 


(6.29) 


, k 1 

lf  X = iT  = 2 
2 * 


The  reader  has  probably  noted  that,  as  yet,  nothing  has  been  said 

as  how  to  select  P =k,  , P = k_  , given  [k  . ,k  ] . This  is 

XX  2 2 min  mdx 

one  of  the  major  issues  if  one  is  interested  in  a synthesis  procedure. 

Let  us  assume  in  the  meantime,  that  k^  and  k 2 have  somehow  been 

selected.  Let  the  relative  change  m(k  , t^)  at  a given  time  t , in 

the  system  response  at  P = k € [k  . ,k  ) , be  defined  as: 

’ r min  max 


m (k , tQ) 


Using  (6.28) 


a c(to)-aci(to) 
aci(to> 


m(k,t0> 


( 


” (tQ)  ) 


k/k. 


(6.30 


As  an  illustration,  m(k,tQ)  is  plotted  versus  k , for  tQ = .2,1.,.. 
in  Fig.  VI. 4. a for  the  case  k^  = 100  > k2  = 1 , and  in  Fig.  VI. 4. b for 

the  case  k^  = 1 < k2  = 100  . The  spread  in  the  overall  system  response 

to  a command  input  r , at  time  tQ  , is  obtained  from  these  curves 
as  follows: 


and 

For 


SPREAD  = c_  (t_)  * [ (m(k,tn>)  - (m(k,t.)) 

l u u max  u 

MIN  are  taken  over  all  possible  values  of 

instance  in  Fig.  VI. 4. a,  when  k . = .1  and 

min 


where  MAX 

min  J 


k e 

k 

max 


[k  . ,k  ] . 
min  max 

= 10., 


[m(k,t  ) ) = m(k  ,t  ) , 

0 max  max  0 


[m(k,t  5]  . = i"<k  • .t  ) , 

o min  mm  o 


and  therefore  it  is  easily  seen  that: 


SPREAD  = c^  (2 . ) * (.99-. 24)  = .72  * 0^2.)  at  tQ  = 2 seconds 

= 0^(9.)  * (.76  - 0)  = .76  * 0^(9.)  at  tQ  * 9 seconds,  etc 

Consider  the  case  k < k„  , shown  in  Fig.  VI. 4. b.  If  k . < k, 

1 2 mm  l 

then  very  poor  results  may  be  expected  for  k < k < k^  , since  the 

settling  time  is  seen  to  be  very  large.  For  the  same  reason, 

k < k . < k„  is  not  a better  choice.  However,  if  k,  < k,  < k . < k 

1 min  2 12  mm  max 

then  we  get  a reasonable  spread  in  the  system  command  responses,  for 

large  changes  in  k 6 tk  , ,k  ] . If  k,  = 1 < k * 100  are  selected 

min  max  2 1 

(T'ig.  VI. 4. a)  then,  on  the  one  hand  k . Ik,  should  be  satisfied,  if 

min  2 
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Figure  VI. 4 Plot  of  m(k,tg)  for  k^ > k^  la)  and  k^  > k^  (b) . 

one  is  to  avoid  a long  "tail"  in  the  system  response  to  command  inputs, 

and  on  the  other  hand,  k,  Sk  should  be  used  to  avoid  tremendous 

1 max 

overshoot  at  small  tQ  . Note  that  all  those  results  are  independent 
of  the  nominal  choice  of  c^  (tg)  , and  of  the  actual  values  of 


k . ,k  and  are  therefore  very  general, 
min  max  1 ^ 


To  summarize,  we  have  found  that  from  the  system  command  response 


viewpoint.  The  curves  (B)  in  dashed  lines  were  ontained  by  takinc 


[k  or  k ] (instead  of  IP  = (k,  or  k ]) 


the  binary  plant 


Obviously,  B passes  through  (k  = k ,0.)  for  all  t 


when  k = k_  in  the 


latter  case,  and  the  noise  effect  is  zero  when  no  uncertainty  (k  = k ) 


VI . 5 . a , that  when  k > k 
max 


tremendous  noise  level  at  the  plant  input  for  all  k..  < k < k 


: . < k 
min  max 


From 


Fig.  VI. 5. b,  it  is  obvious  that  (6.32)  is  also  highly  satisfactory 


The  choice  (6.32)  seems  a better  one  than  (6.31),  especially  if 


k.  and  k are  taken  much  smaller  than  the  prescribed  k 


max 


However,  such  a choice  is  very  poor  as  far  as  stability  is  concerned 


when  higher  order  systems  (obtained  when  dealing  with  more  realistic 


Figure  VI . 5 Mean  square  value  of  the  noise  at  the  plant  input  for  diff 
erent  gain  factor  values  of  k and  at  different  time  instant  tn  when 
i , > k_  (a)  and  k„>  v (b) . 
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[ 


plants)  are  considered.  Indeed,  it  will  be  shown  (section  VI. 4. c)  that 

for  high  order  systems,  k A,  should  be  made  smaller  than  some  real 

max  Z 

number,  preventing  (6.32)  from  being  used.  This  result  is  not  obvious 

here,  because  a first  order  system  is  stable  for  all  values  of  gain 

factor  k . Therefore,  one  has  to  use  the  other  possibility  (6.31) 

which  is  more  compatible  with  stability  problems.  Noting  (Fig.  VI. 4. a), 

that  k . /k.  should  be  as  large  as  possible  to  guarantee  satisfactory 
min  z 

spread  (i.e.,  T.D.S.),  we  will  therefore  use  for  the  remaining  of  this 
chapter : 


k,  = k 
1 max 


k.  f k . 

2 min 


(6.35) 


VI. 2. d Results  and  Some  Preliminary  Conclusions. 

We  can  now  complete  the  above  example  (beginning  of  section  VI. 2. c). 

In  addition  to  the  assumptions  already  made,  we  take: 

C,  (t)  = 1-e  t , k . = 1 , k = 100  and  w(t)  = 1 

1 min  max 

Following  (6.35),  k^ = 100  , k^ = 1 are  selected.  The  design  has  thus  been 
completed  [Appendix  A3. 6. a for  details], and  the  overall  system  responses  to 
unit  step  input  are  shown  on  Fig.  VI. 6. a for  different  values  of  P •»  k . 

The  mean  square  value  of  the  noise  at  the  plant  input  is  shown  in 
Fig.  VI. 6. c where  it  is  compared  with  the  results  obtained  in  a L.T.I. 
design  (Appendix  A3. 8. b) which  acheives  roughly  the  same  time  domain 
specifications  (shown  on  Fig.  VI. 6. b)  as  the  L.T.V.  system.  These 
results  confirm  that  much  may  be  gained  in  reduction  of  sensor  noise 
effects  at  the  plant  input,  by  using  L.T.V.  compensation. 


K - lOO/V"  / 


time 


System  Response  to  a unit  step  command  (a)  LTV  design, 

(b)  LTI  design.  (c)  Comparison  of  sensor  noise  rejection 
at  the  plant  input  between  the  LTI  and  LTV  designs. 


Ill 


(t)  - c,  (t)  , but  this  is  true  only  for  the  simple  first  order 

max  min 

system,  and  not  in  general  for  more  realistic  systems,  even  in  the  case 

where  there  is  no  uncertainty  in  the  dynamics  of  the  plant. 

S)  Note  that  (6.18)  and  (6.22)  do  not  have  h(t,t)  = 0 ,Vk  , 

Vt  , implying  that  our  optimum  design  is  a first  order  one,  which  is 

of  course  not  realistic.  Therefore  the  above  should  be  considered  only 

as  an  academic  excercise,  giving  us,  however,  insight  as  to  how  k^  , k^ 

should  be  chosen,  given  the  uncertainty  range  (k  . , k I 

’ nun  max 

y)  Note  that  h (t,i)  is  a stable  design  for  k € ]0  , ®[  , 


which  will  no  longer  be  true  (as  already  noted)  in  high  order  systems. 


VI. 3 Some  Approaches  to  the  More  Realistic  Problem. 


VI. 3. a Generalities . 

It  is  noted  that  when  the  plant  is  a pure  gain  factor  ( as  in  the 

previous  section) , the  minimization  of  the  sensor  noise  (of  strength 
2 

) effect  at  the  plant  input,  is  equivalent  to  the  minimization  of 

2 2 

sensor  noise  (of  strength  u = o /k  ) effect  at  the  plant  output 

N 

(which  is  also  the  system  output).  Therefore,  the  derivation  in  VI. 2 

can  be  considered  as  a direct  extension  of  Wiener  filter  theory  to 

linear  nonstationary  problems.  As  is  well  known,  the  derivation  of 

the  optimum  filter  in  such  cases,  is  always  done  on  the  closed  loop 

impulse  response  h(t,T)  , which  does  not  explicitly  feature  the 

complexity  of  the  plant  and  therefore  a first  order  optimum  solution 

is  usually  obtained  (with  white  sensor  noise) , which  is  unsatisfactory 

if  the  excess  of  coles  over  zeros  of  the  plant  e >.  2 . (Indeed,  as 

P 

h is  first  order,  so  is  the  open  loop  i , implying  in  the  latter  case 
that  q possesses  at  least  one  derivative  — which  is  impossible  to 
implement  precisely) . 
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In  L.T.I.  systems,  P(s)  and  G(s)  commute,  so: 

+OD 

T (u)  9 — / !T7r'2  $ (w)dcu  * ~- 

P.I.  2-n  J 1 1+L  nn  2n 


1+L 


-foo  $ 

= -f—  / 1 H ! 2 $ (w)doj  with  \ = — — 
2 it  1 nn  P 


.<•  4>  2 

f l-£L'2,-nJll  d<D 


(6.36) 


This  means  that  the  mean  square  value  of  the  noise  effect  at  the  plant- 
input  j (id)  is  equal  to  the  mean  square  value,  at  the  system 


output,  due  to  a distorted  sensor  noise  power  spectrum  0 / V 

nn 


It 


is  therefore  seen  that  H then  emerges  with  the  proper  excess  of  poles 
over  zeros,  provided  the  sensor  noise  is  distorted  by  P 1 . This 

reveals  that  the  problem  posed  in  (6.2)  is  a very  realistic  one,  if  it 
can  be  solved  in  general. 

VI. 3. b The  Problem  Posed  by  the  Solution  of  (6.2). 

It  is  assumed  in  the  remainder  of  this  section,  that  the  input 

r is  deterministic,  implying  that  D^tt)  is  also  deterministic 

(N  = l)  , and  P (s)  = k/s  with  k 6 [k  . ,k  ] . Under  such  assump- 

mm  max 

tions : 


, ° t r 3h, 

(6.1)  becomes  a = ——  J -r£-  (t,;) 

p-x-  k22  0 1 3t 

and  (6.2)  becomes  (using  (6.12)  and  u = cf/k2 

N 2 


dC 


(6.37) 
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Using  a variational  argument,  i.e.,  letting  h2  = h’  + f , where 
h°(t,T)  represents  the  optimal  solution  of  (6.38)  and  f(t,T)  is  an 
arbitrary  small  impulse  response  belonging  to  the  admissible  class  of 
solutions,  i.e.,  satisfying  h(t,t)  =0  , it  can  be  shown  (the  deriva- 

tion omitted  here,  as  the  results  are  not  used,  but  they  are  similar  to 
those  of  Appendix  A3. 5),  that  the  optimum  h2(t,x)  is  the  solution 
of : 


Vf  , Vt  , t>c2 


h2  (t'C)YDD<?l 


C2)dC  - y 


DD 


(t,c2) 


dn2  3f  \ 

f(t,C2)  + It  it  (t,C2)J  dC2  = 0 (6.39) 

Mathematically,  (6.40)  is  equivalent  to: 
t 

/ h2(t't,YDD(C'T)dC  = YDD(t'T>  ' Vt  ' f°r  Tft  <6-39-a^ 

{° 

*■  3h° 

(t , t)  =0  Vt  i x (6.  39. b) 

which  is  unrealizable.  This  can  be  understood  as  follows.  For  t=T 
with  T fixed,  the  optimum  solution  exists  and  since  (6.39.b)  should 
be  satisfied  at  t = T , the  optimum  solution  is  dependent  on  T . Therefore 
at  t = T’  > T , the  optimum  solution  is  a funtion  of  T’  , and  unless 
(6. 39. a)  happens  to  be  satisfied  (which  can  only  occur  by  chance), 
this  optimum  solution  for  t=T'  is  no  longer  optimum  for  t=T  , 
and  therefore  there  does  not  exist  a solution  to  problem  (6.38)  and  in 
general  to  (6.2).  The  reader  may  note  that  when  3h/3t  is  replaced 
by  h in  (6.38),  i.e.,  when  the  noise  effect  is  considered  at  the 
system-output  (rather  than  at  the  plant  input),  then  (6.2)  has  a unique 
solution  (derived  in  VI. 2).  This  may  be  surprising  at  first  glance, 


but  note  that  in  the  latter  case,  (6. 39. a)  does  not  need  to  be  satis- 
fied, thus  explaining  the  uniqueness  of  the  solution. 

One  alternative,  therefore,  is  to  modify  the  original  problem  of 
(6.2)  and  solve:  (for  a given  fixed  T ) 


Min  J (Ac  (t)  + W (t)  o (t)]dt 
h2  0 


(6.40) 


It  is  shown  (Appendix  A 3.5)  that  (6.40)  has  a solution  given  by  an 
Euler-Lagrange  differential  equation.  Another  alternative  is  to  distort 
the  noise  characteristics  in  such  a way  that  the  resulting  h2(t,T) 
will  have  the  proper  excess  of  poles  over  zeros.  For  this  purpose,  we 
introduce  artificially  an  equivalent  autocorrelation  function  of  the 


Wt(T)  = (-1)S  °N  66(t-T) 


(6.41) 


where  e % 1 is  an  integer. 

By  analogy  with  the  L.T.I.  case  (recall  (6.36)),  we  then  try  to 


minimize : 

aout  - / 

in  place  of  (6.1).  Under  this  assumption,  (6.2)  becomes: 

Vt  , Min  { A c2  (t)  + W(t)2o2  } 

out 

n 


(6.42) 


(6.43) 


which  has  a unique  solution  as  shown  below.  It  should  be  recalled 

2 

that  in  L.T.V.  systems,  g and  p do  not  commute  and  therefore  o 

2 

^ Op  . The  weighting  function  W(t)  might  then  be  used  to  try  to 
partially  compensate  for  this.  The  remainder  of  this  chapter  is 
devoted  to  a detailed  discussion  of  the  above  two  approaches. 
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• 4 Solution  to  the  Filter  Problem  by  Distortion  of  the 
Noise  Characteristics. 

VI. 4. a Statement  of  the  Problem. 

The  minimization  of  (6.43)  was  solved  in  (2.5),  giving  for  the 
optimum  solution  h^tt.T) 

t 

tiT  YDD(t,T)  = I h2(t’^)  !YDD(C'T)  +YNN(^'T)Id5  (6.44) 

Using  (6.41)  and  recalling  (6.9)  that  = (Pj-P^l'p^  ^ is 


deterministic  because  the  input  signals  are  assumed  to  be  deterministic, 
we  get: 


t - t 

D2(t)D2(x)  = / h2(t,c)D2(5)D2(T)dC + (-l)eoNW(t)  / h2(t,?)6e(?-T)dC 


(6.45) 


t 3e 

D_  (t)D_  (T)  = / h (t  ,t)  D (O  D,  (T  ) dC  + (-1)0„W  ( t)  — — h.  (t  ,T  ) 
22  \ 2 22  N . e 2 

0 3t 


A necessary  condition  is  therefore: 
_e 


3t 


— h (t ,x ) = - y(t)D  (r)u(t-T) 
e 2 2 


(6.46) 


(6.47) 


Integrating  with  respect  to  r , and  requiring  that 


we  get: 


3t 


3e_1  1 1 

r h (t  ,x ) = y(t)  [ D,(t)  - D,(t)  ] u(t-T) 

e-1  2 22 


3e_1  I 

h,(t,T)  = 0 

. e-l  2 

3x  T = t 


(6.48) 


where  the  notation  D*(t)  = J D (£)dc  is  used. 


Let  D™(t)  = j D™  (5)d£  and  integrating  (6.48)  with  respect  to  T 

9 


2 2 

T h (t,T)  = Y (t)  t D,(t)  -D,(l) 

e-2  2 2 2 


(t-T)D2(t)  ] U(t-T) 


. 


gives : 


more  integrations  with  respect  to  t 


we  have 


D (t)  U(t-T)  = y (t)2(t,T)u(t-T) 


Inserting  (6.49)  into  (6.45)  gives 


(-1) 6 a^W(t)  + / Q(t,C)D  (c)dc 


with 


D,(t)  D,(C)dC 


jD,(Od? 


(t-C) D_  (C)  dc 


/ (t-C) 2D2 (C) dC 


/ (t-;)Dj(()d(  = D^(t) 


D (C)  dC' 


/ Q(t,C)D2  (OdC  = D2(t)D2(t)-  / D2(?)D2(C)de  - d' 


Noting  that 


I1-  ; 

« I 

i 


t r 

- / D®(C)D2<Ud;  = - D*(t)D,(t)  + / D,  <C)D,<C>dC 


2 2 


D®(t)D2(t)  + D®_1<t)D2(t)  - / D®"2(OD^(C)dC 


= - D2(t)D2(t)  + D2_1(t)D2(t)  - D2‘2(t)D2(t)  + ... 


+ (-i)e_2D2(t)D2 (t)  + (-D e_1  / 02(C)D2(C)dC 


(6.51) 


.(6.50)  becomes: 


w « - 

/ Q(t,?)D  (C)d?  = (-1)  / D (C)D®  (C)dC 

0 0 


(6.52) 


so  finally: 


Y(t>  = D2  (t)  / ( (-l)e"1o2W(t)  + (-1) e_1  / D2(OD2(C)d?  ) 


and 


h 2 ( t , "O 


(-l)e_1D  (t)  tD2(t)-D2(T)-(t-T)D®_1(t)+  D2'2(t)  + ...+  (-l)e_1  7^77,  D2(t)l 


02W(t)  + / D (OD®(C)dC 

N 0 2 2 


-U(t-T) 


(6.53) 


(6.53)  holds  for  e = 1,2,3...,  and  it  can  be  checked  that: 

,e-l 


h (t,t)  = (t,t)  = i4(t,t)  = 

3t 


at 


e-l 


h(t,t)  = 0 , 


implying  that  the  system  obtained  is  of  order  e + 1 ( e + 1 poles  and 
no  zeros) . 


i Ik- 


■Si .. 


1X8 

VI. 4. b.  Analysis  for  different  values  of  e 
VI . 4 .b  .a  e = 1 


When  e = l , (6.53)  becomes 


h2<tfT) 


D2(t)  [D2(t)  - D2(T)]u(t-T) 

O^W(t)  +D^(t)2/2 
N Z 


(6.54) 


If  W (t)0  = U is  a constant,  (6.54)  is  associated  (Appendix  A3. 6) 

N 


with  the  Differential  Equation  (D.E.): 


..  r 1 62  1 • r D2°2  °2  2 'D2  1 

[h2J  ! y+  [2V3d7  |y+  [D2«-2Y,t3^  " J y ' D2?X 


a=  y + 


y + 


where 


V + (D2)2/2 


y - D2gx 

(6.55) 


(6.56) 
'1 


Recalling  that  h = £•  (1+K)  1 with  I = p*g  , we  have:  l = p*p2‘*H2 


and: 


h = P,P21,J2(1+P’P21,!.2)"1 

= p-p21-i2- (l+i2)-1  ((l+i2)_1)"1(l+p-p21-ll2)-1 
= p-p21*h2* ( (l+p*P21* l2) (l+&2)  1)  1 
* p,P21*h2* ( (l+i2>  1 + p,P21,h2)  1 
As  (1+S.2)-;I  = 1 - h2 
h » p* p,1*h2* (l+h2 (p«p21-l) ) 1 


(6.57) 


As  we  know  from  L.T.I.  systems,  the  problem  of  sensor  noise  rejection 

at  the  plant  input  becomes  crucial  in  the  region  where  P(s)  has  its 

e 

asymptotic  behavior,  k/s  p , so  we  can  restrict  ourselves  to  gain 
factor  uncertainty  with  p*p^^  = kA.,  = X • The  differential  equation 


associated  with  h = Ah2- (l+h2  (A-l ) ) * is  obtained  as  follows: 


[ (A-l)  h2l 


y + a^y  + C^y 


= (A-l)D2gX 


[l+(A-l)h  ] 


9 + a^y  + a2y  = ( (A-l)  D2g+a2>  X+a^X+X 


[l+(A-l)h2I 


y + a^y  + [a2+ (A-l)  D2g]  y = x+ot  X+a  X 


th2(l+(A-l)h2)  ] : y + a^y  + [a2+ ( A-l ) D2gJ y = D2gX 


y + a1y  +[a2+(A-l)D2g]y  = D2gX 


(6.57.a) 


(6.58) 


In  a practical  problem:  . 

DJ(t) 

lim  D2 (t ) = K , lim  D2(t)  = 0 , lim  — - — = K and 
t=QO  t=°°  t=°° 

2 

lim  g ( t ) • Kt  = 2 
t=°° 

Therefore  the  asymptotic  behavior  of  the  D.E.  (6.58)  is: 


(6.59) 


This  is  a Cauchy-Euler  type  of  D.E.  Through  the  change  of 
variable  t = eU  , (6.59)  becomes: 


V + 3y  +2Ay  = 2AX 


(6.60) 


B.I.B.O.  stability  guaranteed  for  (6.60)  and  therefore  for  (6.59). 
We  can  then  apply  a theorem  by  Cesari  ( [C5 ] , p . 38 ) which  states  that  if 
(6.55)  tends  asymptotically  to  (6.59),  as  t -*  ® , and  if  (6.59)  has 

bounded  solutions,  so  has  (6.55).  Therefore,  we  conclude  that 
V A > 0 , (6.55)  secures  B.I.B.O.  stability. 


VI.4.b.B  e =2 


(6.53)  qives:  h2  (t,T)  = 


D2<t)  (D2(T)  -D2(t)  + (t-T)D2(t)] 


C2W(t)  + / D,(C)D,(OdC 
N o 2 2 


U (t-T) 


(6.61) 


F 
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Following  the  same  procedure  as  above  h(t,T)  is  associated 
with  the  D.E.: 

[ h]  : *y  + o^ftjy  + !*2(t)y  * a3(t)y  = lD2(t)g(t)X  (6.62) 


with 


A 2 

°1  ' 3D29  - 5 D 

A D9(t)  2 

g(t>  £ ^ 

”5*  / D2«'D2“'a«  A ,„1  , j2  „2  , 62  „,52,2 

0 a,  = 3D  g - 7 — D g - 4 — + 12  (— ) 

2 2 2 

a3  =>  XD2g-t  12(D2g)3 -6D2D2g2  - ^ (30  (D2g)  2 - 3D3g)  + 2lD2g  {—■)  2 

2 U2 

^23  ^2  2 °2  ^2  °2 

-!2(-)  -(^,*30,^3^  . 

2 2 2 2 2 


As  t -*« 


, (6.62)  is  equivalent  to: 


...  9 „ l8  • 61  61  w 

y + -V  + -^-y  + — y = x 

C t tJ  t 

Using  (t = eU)  this  Cauchy-Euler  D.E.  becomes: 


(6.63) 


•y  + 6 y + 11  y + 6AV 


6XX 


(6.64) 


(6.64)  has  bounded  solutions,  if  all  minors  A^  (Routh  criteria)  of: 


6 1 0 ! fin 

61  11  6 j > 0 , i.e.  j A > 0 

0 0 61  ! *■  1 < 11 


Therefore,  B.I.B.O.  stability  is  secured  for  (6.62)  if  0 < 1 < ll 

This  imples  that  k < 11  k2  , and  therefore  if  we  have  to  cope  with 

uncertainty  such  that  k > 11 k , B.I.B.O.  stability  is  secured 

max  min 

when  k . < k_  , which  is  a very  poor  choice,  as  far  as  the  system 

mm  2 1 

response  is  concerned,  (recall  (6.35)  and  the  discussion  in  section 
VI. 2. c) . 


..  /-A-— 


K; 

If  *i 

( 


L 
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VI . 4 .b  . y e = 3 


(6.35)  gives : 


D <t>  ID^  <t)  -D^T)  - (t-T  ) (t)  + --]  - (t)  1 

h2(t,T)  = : u(t-T) 


o^W(t)  + / D2(?)D2(C)d H 
0 


(6.65) 


h(t,T)  is  then  associated  with  a D.E.  whose  asymptotical  behavior  is 
given  by: 


fv1  (4)  16  ...  12  96  • 241  241  „ 

[h)  : y + — y +-=•  V + -j  y +—  Y = — T-  X 

t t t t 

which  is  equivalent,  after  the  change  of  variable  t = eU  , to: 
(41  [31 

y 1 1 +10  yl  J + 35  y +50  y +24Xy  = 24XX 

A sufficient  condition  for  stability  is  then  given  by: 


(6.66) 


(6.67) 


> 0 


0 < X < 


25 


10  1 0 0 

50  35  10  1 

0 24X  50  35 

0 0 0 24X 


Note  that  the  uncertainty  range  one  can  cope  with  efficiently  in  such 
a case  is  only  25/3  (recall  (6.35)). 

VI.4.b.6  Analysis  for  e > 3 and  Conclusions 

As  e increases,  we  can  develop  similar  relations.  To  secure 
B.I.B.O.  stability,  one  gets  sufficient  condition 

0 < X < * (6.68) 

with  lim  k =1.  It  becomes  once  more  clear  that  stability  is  conflict- 

e»>on 

ing  with  the  sensitivity  of  the  system  response  to  parameter  uncertainty. 
Note  that  our  filter  is  an  "all  pole"  filter  and  a big  improvement 
might  be  obtained  by  inserting  proper  "zeros".  We  have  not  pursued 
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this  direction,  since  the  filter  is  L.T.V. , and  it  is  rather  difficult 
to  introduce  those  zeros  in  a systematic  way. 

As  for  e > 2 , the  stability  problem  prevents  us  from  coping 

with  large  uncertainty,  the  corresponding  systems  are  not  so  useful. 
However,  for  e=l  , no  stability  problem  occurs,  and  therefore  the  idea 
arises  to  use  the  "second  order  system"  obtained  with  e = 1 , and  to 

doctor  the  "equivalent  transfer  function"  by  inserting  a L.T.I.  far-off 
zero-pole  package  in  order  to  come  out  with  the  desired  excess  of  poles 
over  zeros.  Without  loss  of  generality,  we  used  a package  similar  to: 


K (s+z) 


, 2 „ 2.  n 

(s  +2r  a ) s+w  ) 

P P P 


(6.69) 


where  n , z , ? , ai  are  to  be  suitably  chosen  by  the  designer. 

P P 

VI.  4. c.  Synthesis  for  Large  Gain  Uncertainty. 

VI.4.c.a  Philosophy  of  the  Synthesis  Procedure 

Let  the  uncertainty  range,  (k  . ,k  ] and  the  T.D.S.  be  given. 

min  max 

It  is  found  experimentally  that  if  k ,=  k (recall  (6.35))  is  used, 

1 max 

then  the  nominal  response  c^(t)  is  best  chosen  close  to  the  prescribed 
time  domain  upper  bound. 

Using  (6.35)  k = k . is  first  chosen.  This  leads,  in  general, 
z min 

to  the  slowest  possible  response  signal,  which  may  or  may  not  be 
satisfactory.  If  not  satisfactory,  k2  < ''min  is  adjusted  until 

the  response  signal  for  k = satisfactory.  This  can  be  under- 
stood physically  from  the  fact  that  specifying  a larger  uncertainty 

range  (k_  ,k,]  than  the  actual  [k  . ,k  ) means  that  i(t,T) 
z l min  max 

handles  a bigger  burden  due  to  uncertainty,  which  then  leads  to  an 
increase  in  its  "Bandwidth"  and  thus  a smaller  spread  Ac(t)  for  the 
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considered  ranqe  [k  . ,k  ] c:  [v  , k,  ] 
min  max  2 1 

As  a consequence,  the  sensor  noise  effect  at  the  plant  incut  is 
increased,  but  there  is  still  preserved  a certain  improvement  over 
L.T.X.  design,  as  seen  from  the  examples  below.  Note  that  the  factor 
W(t)  in  (6.1)  can  also  be  used;  however  this  is,  in  aeneral , only 
efficient  at  small  t 

In  a general  sense  k = k . does  correspond  to  the  lower  Time 

min 

Domain  Bound  but  k = k does  not  correspond  to  the  upper  one, 
because  as  k -»•  k^^  , corresponding  signals  have  the  tendency  to 

oscillate  around  the  nominal  signal  c^(t)  (see  for  example  Fig.  VI. 8 
and  VI. 10).  Therefore,  if  one  is  to  acheive  T.D.S.  this  can  certainly 
be  done  by  cut  and  try,  i.e.,  once  k 2 has  been  selected,  to  correct 
the  nominal  signal  c^(t)  in  order  to  satisfy  those  T.D.S.  and  then 
derive  a new  k^  to  satisfy  the  lower  bound  ar.d  so  on.  The  above 
'cut  and  try'  is  always  inevitable  whenever  a 'trade-off'  optimization 
criterion  is  used  as  the  basis  of  design,  whereas  the  primary  specifi- 
cations are  in  the  form  of  performance  bounds.  We  are  at  least  finding 
a systematic  cut  and  try  procedure. 

In  order  to  determine  the  L.T.I.  pole  zeros  package,  the  concept 
of  "frozen-time  open  loop  transmission"  is  used.  This  means  that  for 
each  fixed  time  t^  , the  L.T.V.  differential  equation  becomes  a L.T.I 
one  and  thus  one  can  define  the  associated  L.T.I.  open  loop  transfer 
function  A(s,t^)  (Fig.  VI. 7).  We  are  aware  of  the  fact  that  L.T.V. 
system  can  be  unstable  although  each  of  its  associated  A(s,t0)  leads 
to  a stable  design,  which  means  that  the  tool  used  here  is  certainly 
not  a very  accurate  one.  However,  it  turns  out  that  this  concepts  can 
be  applied  successfully  here  and  therefore  we  did  not  investigate  more 
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accurate  methods  like  the  application  of  the  EPLTI  concept  (see 
chapter  VII) . 

The  package  is  chosen  to  give  a very  small  contribution  to 

A(s,tg)  in  the  low  frequency  range.  Therefore,  it  suffices  to  consider 

the  worst  case,  i.e.,  the  A (s,t  ) with  the  biggest  bandwidth 

M 0 

(Fig.  VI. 7).  Thus  it  is  conceivable  that  a 'time-varying  far-off  pole 
i zero  package'  may  be  better. 

VI.4.C.6  Example  and  Results. 

Assume  a set  of  step  command  inputs  of  amplitude  <*  , with 

I 2 

<„.>  . = 1 , and  let  the  T.D.S.  be  as  shown  in 

ui  Average  over  l 

I Fig.  VI. 8. c.  The  plant  P (s)  = k/s  is  considered  with  k£  11,1001. 

| 

1 Furthermore,  it  is  assumed  that  the  strength  of  the  white  sensor  noise 

2 

is  a = 1 . Following  the  above  procedure,  the  nominal  system  response 

i N 

-t  2 

Cj(t)  = (1-e  ) is  selected  and  paired  with  = kmax = 100 

k = k . =1  and  W(t)  =1  are  then  first  chosen.  The  corresponding 

2 min 

system  response  to  a step  command  are  then  shown  in  Fig.  VI. 8. a for 
different  values  of  the  gain  factor,  and  it  is  seen  that  the  T..D.S.  are 
’ not  satisfied.  The  sensitivity  of  the  system  response  to  plant 

uncertainty  can  be  decreased  by  decreasing  W(t)  , giving  then  less 
j emphasis  to  the  noise  performances.  For  instance  with  W(t)  = .01  (Fig. 

► 

vi. 8. b)  over-design  is  acheived.  By  cut  and  try  W(t)  *.l  is  found 

i 

satisfactory  and  the  system  response  to  a step  command  are  then  shown  in 
Fig.  VT.fl.c  for  different  values  of  the  gain  factor  k 

Note  that  the  plant  is  first  order,  therefore  a far  off  zero  pole 
package  i s not  needed  in  the  present  case  since  the  closed  loop  trans- 
mission, described  by  (6.58)  with  all  terms  in  D/D  and  D/D  dropped, 
is  second  order,  i.e., 


VI. 8. a (W(t)  = 1. ) 


VI. 8. b (W(t)  = .01) 


Figure  VI.8.a  LTV  system  response  to  a Figure  VI. 8. b LTV  system  response 

command  input  for  different  plant  gain  to  a command  input  for  different 

factors  when  W(t)  =1.  [P(s)  = k/s] . plant  gain  factors  when  W(t)  = .01 

[P  (s)  = k/s]  . 


Figure  VT.8.C  LTV  system  response  to  a command  input  for  different 
plant  gain  factors  when  W(t)=.l  [P(s)=k/s] 


« y + 


2D2D2  • MD2)2  2D2°2 
■ — y + — - ■ ■ — - ■ - — - — y 

U+(D*)2/2  U+(D2r/2  u-MD2)2/2  J 


Vki 

where  here:  D„(t)  = c 

2 kl 


, t 

x(t)  , D2(t)  = / D2(C)d;  , 


U = W(t)cr2  / k 2 = .1  and  A = k/k_ 

N l <£ 


U+(D2)2/2 
(6. 69. a) 


From  (6. 69. a)  and  recalling  that  £ = h(l-h)  , we  get: 

...  •• , 2d!d2  . 2ok  "v2 

U+(D2)  /2  U+(D2)  /2  u+(D2)  /2 

The  plant  is  [pi :y  = kx  , so  [p  : ky  = x and  as  g = p 
we  have : 


, , . 2 2 . 2 2 1 2 
[g]  : y + J— - — y - — — y = — — - — 

U+(D*r/2  P+(D2)  n 2 U+(D2)  /2 


X (6.69.b) 


As  c (s)  = , «e  have:  ( c,  = t,  r ) 

1 s (s+1) (s+2) (s+3)  1 1 


T (s)  = 

V ' (s+1)  (s+2)  (s  + 3) 


and  T^ts)  is  therefore  associated  with  the 


[t^]  : y + 6 y + 11  y +6y 


(6.69.c) 


By  definition,  h^  (obtained  with  k = k^  ) is  associated  with  t^  , 
therefore,  the  prefilter  f is  such  that: 


Vf  = h 

f = 


(6.69.d) 


Using  (6. 69. a)  and  (6.69.c)  we  get: 


kl  <D2)2  - - • f 2D2D2  • 

[f  1 : — r— - — [ y + 6y  + lly  + 6y  ] = 6 X + - — - — X 

S ti+frr»^/s 


2 u+(D‘)72 

+ x (\  -(P2)2  . 2D252  \ 1 

^ k2  u+(D2)2/2  u+  (D^)  2/2  ' 


U+  (Dj)  /2 


(6.69.e) 
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The  mean  square  values  of  the  sensor  noise  effect  at  the  plant- 
input  are  shown  on  Fig.  VI. 9 for  different  plant  conditions.  The  results 


are  also  given  for  the  L.T.I.  design  [Appendix  A3.8.b]  which  achieves  the 


L.T.V.  and  L.T.I.  design  (P(s)  = k/s) . . 
same  T.D.S.  shown  in  Fig.  VI. 8. b. 

There  is  clearly  a big  improvement  in  the  noise  level  at  the  plant 
2 

input.  Note  that  ^ (t)  -*•  0 as  t -*■ 00  and  for  t -*•  0 . This  can 

bo  understood  from  the  fact  that  as  t ■* 00  , the  L.T.V.  network  behaves 

like  an  open  circuit  and  as  t»<0  the  signal  e (t)  (Fig.  VI. 1)  which 


Figure  VI. 10  LTV  system  response  to  step  command  inputs  for  different 
plant  gain  values  (P(s)  = k/s2)  . 


Figure  VI. 11  LTI  system  response  to  step  command  inputs  for  different 
plant  gain  values  (P(s)  = k/s2)  . 


(Fig.  VI. 7)  the  bandwidth  of  our  system  has  a maximum,  leading  to  an 


xtremum  in  the  noise  level,  both  occuring  roughly  at  the  same  time  as 


The  trade-off  between  noise  level  and  bandwidth 


the  maximum  of 


Figure  VI. 12  Mean  square  value  of  the  sensor  noise  effect  at 
the  plant  input  for  different  values  of  the  gain 
value  [P (s)  = k/s2l  . 


of  a system  needed  to  handle  fast  varying  signals  is  once  more  very 
clear. 


Let  us  now  consider  the  plant  P(s) 


k/s 


with 


k € [1  , 1001 


with  the  same  set  of  command  inputs  as  before.  The  T.D.S.  are  shown  on 

Fig.  VI. 10  for  this  case.  The  nominal  output  c^(t)  = (1-e  t) 3 is 

once  more  selected  and  associated  with  k = k,  =k  =100  . Usinq  the 

1 max 

above  procedure,  k = k . =1  is  taken,  and  by  cut  and  try  W(t)  = .1 

2 mm  J 

is  found  satisfactory  for  the  T.D.S.  to  be  matched.  However,  as  the 

plant  is  second  order,  we  need  here  a far-off  pole-zero  package  in  order 

to  be  able  to  realize  the  compensation  G . n = 1 is  selected  in 

2 

tl)  ('-  + Z) 

(6.69)  and  the  package  used  is:  — • ■ - — . The  addition  of 

z 2 2 

(s  +2 C oo  s+co  ) 

P P P 

such  a package  should  not  affect  too  much  the  signal  responses  obtained 
with  the  "second  order  system".  Therefore,  we  took  as  a criterion,  a 
phase  lag  of  at  most  5°  at  the  maximum  crossover  frequency 

= 20  rps.  (Fig.  VI. 7)  in  the  family  of  A(s,t  ) depicted  on  that 

max 

figure.  The  ’best'  package  is  then  found  to  be  the  one  for  which 
z = 70  rps.  = .4  , “p  = .75  * z . In  a similar  way  to  the  previous 

example  (equations  6.69.a-e)  the  D.E . for  f and  g can  be  obtained; 
this  step  is  skipped  here  to  avoid  lengthy  derivations.  System  output 
responses  are  plotted  in  Fig.  VI. 10  for  different  plant  conditions. 

Using  then  the  T.D.S.  shown  on  Fig.  VI. 10,  a L.T.I.  system  was  designed, 
for  which  the  system  responses  to  command  are  shown  in  Fig.  VI. 11. 

For  both  designs,  the  mean  square  value  of  the  sensor's  noise  at  the 
plant  input  is  plotted  in  Fig.  VI. 12  for  the  two  extreme  cases. 
Improvement  is  seen  for  all  k € [1  , 100] 
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VI. 5 Solution  to  the  filter  problem  by  means  of  the  Euler- 
Laqranqe  differential  equation. 

VI . 5 . a Statement  of  the  problem. 

It  was  explained  in  section  VI. 3. a,  that  the  minimization  problem 
of  (6.2)  has  no  solution  and  miqht  be  replaced  by  that  qiven  in  (6.40), 
which  has  a solution  (Appendix  A3. 5). 


Let 

P(s) 

= with  uncertainty  in  both  a and 

s+a 

in 

k 

Let 

P1  - 

k^/(s+a^)  be  paired  with  output  siqnal 

ci 

(t) 

For 

P2  = 

k2/(s+a2)  we  want  to: 

Min 

/[ 

^ 2 2 1 
(D ( t ) - J h (t,C)D(C)dO  W(t)  +o  (t)  | 

dt 

h2 

o *■ 

o 2 p-1‘  J 

2 

a 

P. 

,I.(t) 

is  defined  in  (6.1),  with  (6.1a) 

°(f"C)  (^h2(t'5)+a2h2(t'a)  and  YNH(;i'52)",!<(CrC2) 

2 

o 

A N 

(6.70)  becomes  ( u = — — ) 

k2 

t f t 2 t 3 

Min  / (D  ( t ) - / h (t,C)D(C)dC)  W(t)+p  / [ ^ h2  (t , ? ) + (t , C)  1 dC 

(6.71) 


h2  0 


As  shown  in  Appendix  A3. 5,  if  a^=0  , (6.71)  is  equivalent  to 

the  Euler-Laqranqe  D.E. 

2 

3 h2  t 


(t,i)  + D ( t ) [D (t)  - / h (t,C)d;] 
3t  0 


(6.72) 


No  closed  form  solution  exists,  in  qeneral,  and  (6.72)  must  be  solved 
by  numerical  methods.  This  is  certainly  a direction  to  pursue  in  some 
future  research,  but  as  we  want  to  avoid  numerical  methods,  we  assume 
h ^ ( t , t ) has  the  form: 


dt 
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h2(t,T)  = A(t)  (6  ( t ) - 9 (T ) )u(t-T)  (6.73) 

whore  Alt)  is  unknown  and  6(t)  is  given. 

VI. 5. b Solution  to  the  filter  problem  and  outline  of  the 
synthesis  procedure. 

As  shown  in  Appendix  A3. 7,  when  (6.73)  is  inserted  into  (6.70), 
the  Euler-Lagrange  differential  equation  is: 


a ^ te2(t)  - 26  (t)e1  (t)  + / e2u)d t ^ + 2A9  (t)  ^ te  (t)  - (t)^ 

a(  e (t)  ( te ( t ) - e,  (t) ) - a2 ( t@2 (t)  — 20 (t) 0 ( t)  + / 02(odc')  - J2(t)^ 

V 12V  1 o / V ) 


D(t)W(t)J(t) 

u 

t 

where  J(t)  = / (0  (t)  - 6 (O  )D(C)d?  and 
0 


0 <t> 


= / ©(Ode 

0 


(6.74) 


(6.74)  is  a singular  differential  equation  which  has  in  general  at  least 
one  unstable  solution,  because  the  coefficient  in  A is  always  negative, 
while  those  in  A and  A are  both  positive.  To  extract  the  stable 
solution  involves  here  too,  numerical  methods,  which  we  want  to  avoid. 
However,  it  seems  that  by  a proper  choice  of  A(0)  and  A (0 ) , we 

can  obtain  quasi-stable  solutions  of  (6.74),  that  is,  solutions  which 
can  be  considered  as  stable  over  some  finite  interval  [0,T]  , big 

enough  for  our  purposes. 

The  synthesis  procedure  is  then  very  similar  to  that  described 

in  section  VI.4.d.u.  c^(t)  is  chosen  to  lie  within  the  prescribed 

T.l>.S.,  close  to  the  upper  bound  and  is  paired  with  corresponding 

to  the  highest  gain  value  k = k .We  then  try  first  k ,=  k , 

max  7 1 max 

k k . An  before,  a cut  and  try  procedure  is  then  used  to  find 

2 mi  rt 

k^  and  W(t)  such  that  the  output  signals  corresponding  to 
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V £ [k  . ,k  ] lie  within  the  T.D.S. 
nun  max 

VI. 5. c Examples  of  application  and  results. 


VI. 5. c. a p(s)  = k/s 

The  plant  P(s)  = k/s  with  uncertainty  in  the  gain  factor 
k 6 [1,100]  is  considered  here.  This  example  was  already  considered 
in  section  VI.4.d.B.  The  T.D.S.  shown  on  7ig.  VI. 13  are  to  be  satis- 


I 

- 1 
i 

I I 

I ; 

i 


j 


fied  for  system  step  command  inputs.  The  nominal  output  response  to  a 

unit  step  command  is  taken  as  c^(t)  = (1-e  t)^  and  is  paired  with 

k - ’ - k = 100  . k_  = k . =1  is  then  selected,  0 (t)  in  (6.73) 

! rux  2 min 

• k2_kl 

is  taken  as  0(t)=D1(t)  , (0(t)=D(t))  , where : D (t)  = — - — c (t)  , 

1 ^ i 

2 2 2 

and  it  is  assumed  that  o = 1 , implying  that  y = o„/k.  = 1 In 

N N 2 

a manner  similar  to  the  procedure  described  in  VI.4.d.0.  W(t)  =1  is 
i r :_.t  chosen  and  by  cut  and  try,  W(t)  = .1  is  found  to  be  satisfactory. 
" erefore  the  ’optimum'  system  (under  the  present  assumption)  is 
' /nod  from  (6.73)  as  h^tt,!)  = A(t)(D^(t)  - (t  ))  u (t-r)  with 

) the  solution  of:  (recall  6.74) 

^t.D^(t)  -2D1(t)D2(t)  + / D^(c)d^  + 2AD(t)^tD1  (t)  -D^t)^ 

/.  D*(t)  \ D(t)D*(t) 

+ A ^D  (t)  (tD^  (t)  - D2  ( t ) ) - .1  j = - .1  — 
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so  q = ^2***2  assoc^ate<^  with  the  differential  equation: 


A /„i,2+AD_A 

AD  A 


(ql:  y-y^2  = - 5 ; + y ^(j?‘  + Xd  k7  x 


^ _ AD 
’ C2 


The  prefilter  f is  obtained  by  satisfyinq 
l • (l+ft  )_1-  f = t 


f = (1+21) •l‘1‘t1  = (l+t*1)*t1 

where  t is  the  specified  overall  transfer  function  associated  with 


Cl(s)  8 

P = P!  ’ 1’B"  T1<S>  “ R(s)  ( s+ 1 ) (st-2)  (s+4) 


in  the  present  case. 


The  system  response  to  step  command  inputs  are  shown  on  Fiq.  VI. 13 
for  different  plant  qain  factors.  It  should  be  noted,  as  stated  above, 
that  all  these  solutions  are  "unstable"  at  infinite  t , but  can  be 
reqarded  as  "stable"  over  some  finite  interval  of  time.  This  is  not 
inconvenient  if  one  has  to  satisfy  T.D.S.  over  a finite  interval  [0,T]  , 
as  is  often  the  case  with  L.T.V.  desiqns  (recall  that  L.T.V.  desiqns 
do  not  qive  identical  results  for  [0,T]  and  (a,  T+a]  ). 

2 

The  effect  of  white  sensor  noise  (of  power  spectrum  0^=1.  ) 

is  shown  on  Fiq.  VI. 14  for  different  plant  qain  factors,  and  is  com- 
pared with  the  results  qiven  by  a L.T.I.  network  which  achieves  the 
same  T.D.S. 

Vl.S.c.8  P(s)  = k/s+a 

Consider  the  plant  P(s)  = k/(s+a)  with  uncertainty  in 

a € (.1  , 1.]  and  with  no  uncertainty:  in  the  qain  factor.  Let  us 

select  k2  = 1.2  , k2  = l.  , and  assume  that  c^(t)  “ (1-e  is  paired 

2 

with  P, (s)  = k,/s+]  . Let  0=1  and  let  the  T.D.S.  to  step  command 

11  N 

input-,  bo  .is  shown  in  Fiq.  VI. 15.  Once  more,  by  cut  and  try,  W(t)  = .1  , 
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Figure  VI .14  Mean  square  value  of  the  noise  at  the  plant  input  due 
to  white  sensor  noise  for  the  LTV  and  LTI  design, 
is  found  satisfactory  giving,  for  different  values  of  a , the  responses 

to  a step  command  shown  in  Fig.  VI.15i  g and  f can  be  obtained  as 

in  the  above  section.  (The  details  are  therefore  skipped  here.) 

It  should  be  underlined  that  our  synthesis  technique  cannot  cope 

with  large  uncertainty  in  both  the  gain  factor  k and  the  pole  a 


(Note  that  this  was  also  true  with  the  first  approach.) 


Figure  VI. 15  LTV  system  response  to  a step  input  command  for 
different  values  of  the  pole  a 


Figure  VI. 16  Typical  LTV  system  response  to  step  command  inpute 
for  different  plant  conditions.  P(s)  = k/(s+p) 


This  is  due  to  the  fact  that  our  optimum  design  is  very  sensitive 


max 


because  it  is  "tuned"  to  the  two  extreme  cases 


does  not  help) 


(interchange  of  a 


max 


This  means  that  both  approaches  give  correct  system  responses  [Fig.  VI. 16) 


when  a « a 


13ft 


| 


I 

t 


! 
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at  k=k  . , it  is  seen  [Fig.  VI. 16]  that  an  intolerable  undershoot 

min 

occurs  for  some  a € (a  . ,a  ) which  increases  with  a . This  is 

min  max 

due  to  the  fact  that  the  "equivalent  open  loop  transfer  function"  needs 

more  "gain  factor"  when  a = a than  when  a = ct  . , in  the  time 

2 max  min  

interval  where  undershoot  occurs.  Suppose  that  by  means  of  the 
weighting  functions  W(t)  , the  "gain  factor"  of  the  "open  loop 
transfer  function"  is  increased  to  handle  the  undershoot  problem.  We 
then  have  a L.T.V.  system  which  no  longer  has  better  noise  performances 
than  a L.T.I.  design  which  achieves  the  same  T.D.S. 

VI. 6.  Conclusions 

This  chapter  has  presented  a first  attempt  at  using  optimum 
filter  techniques  to  handle  the  uncertainty  problem  by  means  of  LTV 
compensation.  While  the  problem  has  not  been  fully  solved,  in  the 
sense  of  providing  systematic  design  procedures  for  systems  with  large 
plant  uncertainty,  a respectable  beginning  has  been  made  for  achieving 
TDS  despite  large  gain  factor  uncertainty.  It  was  emphasized  that  the 
minimization  problem: 

Min(Ac(t)2  + W(t)o2  <t) ) , Vt 

h2 

had  no  analytical  solution  for  realistic  plants  (of  order  nil  ) . 

Two  alternatives  were  then  presented,  leading  to  realistic 
solutions  with  substantial  improvements  over  LTI  design  in  terms  of 
sensor  noise  rejection  at  the  plant  input.  We  can  conclude  that  this 
chapter  gave  us  a better  understanding  of  the  achievements  and  the 
limitations  that  can  be  expected  from  LTV  networks,  when  used  to  cope 
w;tn  the  uncertainty  problems.  Considerable  work  has  yet  to  he  done 
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to  arrive  at  a comprehensive  synthesis  theory  for  large  and  general 
plant  uncertainty.  A possible  direction  is  try  numerical  methods, 
inasmuch  it  was  noted  that  no  closed  form  solution  can  in  general  be 
expected. 

It  is  seen  that  the  improvements  were  much  more  spectacular  with 
N.L.  compensations  (recall  FORE),  as  might  have  been  expected.  One 
price  paid  is  obviously  in  terms  of  stability,  since  the  stability 
problem  is  usually  more  difficult  in  non-linear  systems  than  in  LTV 
systems.  Also,  one  must  design  for  specific  input  classes  in  non- 
linear systems  and  must  carefully  check  for  the  system  response  to  the 
occasional  signal  belonging  to  other  classes.  While,  here,  we  also 
tailor  the  design  to  a specific  class  of  inputs,  one  can  much  more 
easily  determine  its  response  to  other  input  classes.  However,  it 
should  be  mentioned  that  a significant  restriction  of  LTV  systems  is 
that  one  must  know  precisely  when  a specific  input  (of  a certain  class) 
begins  in  time. 
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CHAPTER  yij  . NOJlLimAJjJXEDBACjC  SYNTHES  I S_FgR  SYSTEMS  WITH 

LARGE  PLAEX  XGXQP^NCE^  FOR  PRESCRIBED  TIME 
DOMAIN  SPECIFICATION'S  OF  A I NEAR  T,VPE " . 


VII. 1.  Introduction ■ 

Vll.l.a  Generalities . 

It  was  shown  in  chapter  II  (section  4)  that  some  classes  of  N.L.T.V. 
plant  W can  be  characterized  by  a set  Pa[s]  , denoted  as  the  EPLTI , 
which  is  associated  with  the  system  input  ia(t)  and  defined  over  a 
set  of  acceptable  system  output  Ca[t]  . This  equivalent  character- 
ization of  the  plant  (that  was  proven  by  application  of  Schauder's  fixed 
point  theorem)  lends  itself  very  easily  to  quantitative  specifications 
and  is  therefore  a powerful  tool  in  the  synthesis  of  nonlinear  feedback 
systems.  This  technique  has  been  applied  successfully  to  the  quantitative 
synthesis  of  feedback  systems  that  included  L.T.V.  or  static  nonlinear 
uncertain  plants.  This  chapter  presents  the  first  application  of  the 
EPLTI  set  concept  to  dynamic  nonlinearities. 

The  synthesis  technique  for  L.T.I.  plants  that  was  reviewed  in 
chapter  II  (section  1)  implicitely  assumes  T.D.S.  of  a linear  "type", 
i.e.,  if  r^  denotes  a command  input  applied  to  the  system  of  Fig.  VII. 1 


-I 


Figure  V r 1 . 1 . Feedback  stiucture  vPn  nonlinear  uncertain 
planf  set-  W 
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and  ci  its  corresponding  output.,  then:  = Xr2  » c (t)  =Xc2(t)  (7.1). 

However,  it  is  conceivable  that  the  T.D.S.  are  of  a different  nature, 
i.e.,  linear  time  varying  (L.T.V.),  nonlinear  time  invariant  (N.L.)  or 
nonlinear  time  varying  (N.L.T.V.). 

a).  T.D.S.  of  a "L.T.V.  type"  means  that  the  T.D.S.  on  the  system 
impulse  response  are  different  for  different  instant  of  application  of 
the  impulse.  This  problem  has  been  discussed  elsewhere  [H8] , [H6] . 

S).  T.D.S.  of  a "nonlinear  type"  means  for  example  that  the  specifi- 
cations are  different  for  different  signal  amplitudes  or  the  T.D.S.  for 
a ramp  input  are  not  necessarily  the  integral  of  those  for  a step,  or 
more  generally  given  (T.D.S. K for  r^  and  r = g*^  , then 
(T.D.S. ) j ? g* (T.D.S. ) ^ . This  problem  is  considered  here. 

v)  . T.D.S.  of  a "nonlinear  time  varying  type " combines  a)  and  8)  . 

The  solution  to  this  problem  is  therefore  a combination  of  the  solutions 
to  a)  and  8)  . 


Vll.l.b  Nonlinear  time  domain  specifications  (N . L. T.D.S .) . 


Fig.  VII. 2 shows  the  acceptable  response  tolerances  for  a 


Figure  VI I. 2.  "Linear"  and  "nonlinear"  tolerances  on  steo  response. 
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unit  step  response  r(t)  =u(t)  . In  a linear  system,  those  for  r(t)  = 
,5u(t)  are  one  half  of  the  former,  as  shown.  In  a nonlinear  system 
they  may  be  as  shown  in  Fiq.  VII. 2,  of  the  kind  usually  associated  with 
time  optimal  response.  It  is  conceivable,  for  example,  that  one  might 
want  the  unit  ramp  response  to  differ  significantly  from  the  integral 
of  the  unit  step  response,  which  is  impossible  in  a linear  system.  One 
should  note  at  this  point  that  there  is  considerable  ambiguity  with 
N.L.T.D.S.  so  defined.  Indeed,  with  L.T.I.  system  the  response  specifi- 
cations are  the  same  when  the  command  input  jumps  from  0 to  1 , 

1 to  2 or  from  5.2  to  6.2 

This  is,  however,  no  longer  tree  for  N.L.T.D.S.  and  therefore  all 
the  specifications  that  will  be  considered  here  are  relative  to  the 
reference  r = 0 . Obviously,  one  can  also  incorporate  other  references 

than  zero  by  adding  some  other  contraints  which  then  increase  the 
complexity  of  the  problem. 

It  is  shown  that  it  is  possible  to  achieve  N.L.T.D.S.  of  the  type 
of  Fig.  VII. 3,  by  using  a L.T.I.  compensation  G(s)  (Fig.  VII. 1)  and  a 
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nonlinear  pr^filter  F 

VII. l.c  Nonlinear  prefilter. 

The  problem  raised  by  the  synthesis  of  such  a nonlinear  prefilter 
F introduces  us  to  the  research  area  of  "open  loop  nonlinear  synthesis", 
for  which  only  primitive  solutions  exist  up  to  now.  The  problem  is: 
given  a finite  set  of  input  signals  R = {r},  design  F such  that  to 
each  input  r ^ there  corresponds  a specified  output  u^tt)  (Fig.  VII. 1). 
One  obvious  necessary  condition  for  the  existence  of  such  an  F is  that: 
if  r^tt)  = r..(t)  for  t € [ 0 , t ^ ] then  u^(t)  = u^(t)  for  the  same 
interval . 

VII. 2.  Synthesis  procedure. 

VII. 2. a Philosophy  of  the  design  procedure. 

Following  section  II. 4,  the  a - EPLTI  transfer  Function  set 
associated  with  input  ia(t)  6 I over  the  admissible  output  set  ca[t] 
and  the  nonlinear  plant  set  W , is 

p0l[s]  = {PU  (s)  , z0  € Za  [t]  , w . € W}  , which  is  well  defined 
1 v v 1 

for  all  frequencies  u>  € fl  = [0,“[  . It  was  shown  [H2]  that  frequency 

domain  specifications  suffice  to  guarantee  T.D.S.,  therefore  the 
synthesis  technique  uses  frequency  domain  concepts. 

In  the  frequency  domain,  there  is  associated  with  each  system 
input  i^tt)  € I a set  Q(u)  of  permissible  equivalent  overall  transfer 
function  Ta(jll0)  [see  [S2|].  The  synthesis  technique  for  L.T.I.  systems 
that  was  discussed  in  chapter  II  (section  1)  can  now  be  applied  and 
Equations  (2.1)  and  (2.2)  are  rewritten  here  as 

I G ( ju>)  p’tjw)  I 

v e - sup  Ain  I ! * A In!  Ta(jw)|  (7.2) 

1 l+G(j(J,)F^(ju.) 


mm m 
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rhe  inputs  considered  are  steps,  so  FT(s)  = X FL(Xs)  and 


C.  (s) 

T.(s)  4 ^ 

1 


R.  (s) 
1 


X C. (Xs) 
1 

X^R. (As) 
1 


= T.(Xs) 
i 


(7.5) 


Thus  it  is  only  necessary  here  to  concentrate  on  the  T.D.S.  of  a 
particular  step  command  (say  r^  = 1)  and  translate  them  into  frequency 
domain  specifications  (Fig.  VII. 5). 

The  frequency  domain  specifications  associated  with  r^(t)  *=  Xr^(t) 
can  then  be  deduced  from  (7.5)  , (Fig.  VII. b, 

PLANT  W - is  given  by: 


W : u -*  z k^z  + k^z  = u 


(7.6) 


with  uncertainty  k^  , k^  € [1,10]  independently. 

VII.2.b.B  Characterization  of  the  EPLTI  set. 

The  acceptable  system  output  C^ft)  is  represented  here  by  a 


second  order  system  step  response,  in  order  to  obtain  an  analytic 
expression  for  Pa[s],  For  step  inputs  of  amplitude  r , 
a-Cn“>nt 


C (t) 
r 


_ A e sn”n"  • - r,  ~ 2 

r 1 1 


v 


n sin  (u>  / 1-C  t + cos 

2 n n 


Our  T.D.S.  translate  into: 


.2<C  i»  <4.  and  . 5 < w /1-C  < 2 . 

n n n n 


(7.7) 


If  ui  = u)  /TT"2 *  , then: 

On  n 


C (t)  = 
r 


r2  / 1+  — 

V 2 (1-C 


- e"2^n<unt  + 

n2)  2C1-V> 


e 2^nt0nt  Cos2u)0t 


- 2e  n n (cosw^t  + 


n 


C _ \ 

, sinu  t)  + ---n  „■  e n<1)n  sin2w  t)  (7.8) 

ATT2  0 /TT2  °J 

n n 


and  •'  (C  2(t ) ] ^ C 2 (s ) = 
r r 


r22w  4(3s+4C  u > ) 
n n n 


2 2 2 2 

s(s+2C  u)  ) (s  +2C  co  s+co  ) (s  + 4C  u)  s+4oo  ) 

nn  nnn  nn  n 
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Using 


u(s)  = 


(7.6) 


2 

roi 

n 


2 2 
s(s  + 2r,  m s+tu  ) 
n n n 


r 

i 


k2S 


2 

r2 a)  k.  ( 3s+4£  w ) 

n 1 n n 

2 2 
(s+2£  to  ) (s  +4C  u)  s+4oj  ) 
n n n n n 


1 

J 


(7.9) 


a Cr(s) 

and  P,s,=— 


2 2 
(S  + 2C  hi  ) (s  +4£  hi  S+4 hi  ) 
n n n n n 


2 2 2 

k„s(s+2£  hi  ) (s  +4t  hi  s+4w  )+2oj  k,r(3s+4c  co  ) 

2 nn  nnn  nl  nn 


(7.10) 

with  £ ,ai  given  by  (7.7)  and  for  a = 1,2,..  5, 6.  EPLTI  sets  are  shown  in 
n n 

Fig.  VII  6. a for  various  hi  and  a . Note  that  for  small  hi  the 

EPI.TI  set  is  composed  of  two  disjoint  regions  (for  +r  and  -r) 

2 

because  of  the  term  in  z in  W (see  7.6).  Such  phenomena  does 
obviously  not  occur  with  L.T.I.  plants.  At  very  high  frequencies, 
howi  vor,  the  equivalent  templates  (Fig.  VII. 6. a)  tend  toward  a straight 
line  similar  to  L.T.I.  plants,  due  to  the  high  frequency  gain  factor 
uncertainty. 

In  order  to  characterize  the  plant  W for  d = ±1  , ±2  , +4 
wc  used  a second  order  model  system  step  disturbance  response, 
leading,  for  P^ts)  , to  an  identical  formula  as  (7.10),  with  sub- 
script d in  place  of  r , and  with  r = -d 

How  are  r and  u.  then  chosen? 
d d 

The  specification  on  the  maximum  tolerable  overshoot  is  trans- 
lated, ,ir  usual  (see  Appendix  Al)  into  a minimum  tolerable  damping 
value.  For  instance,  here  i .36  for  30«  maximum  overshoot,  or 
equivalently  to: 


sup 

i 


G(  ju)P  a(  jm) 

1 +G  ( jid ) P / ( jaj) 


< 3.4  db 


(7.11) 


whore  i ranges  over  all  possible  plants  and  ft  = 7, 8..., 12. 


PLANT  TEMPLATES 


Figure  VII. 6.a.  Templates  at  ui  = .5,2,10  rps . 

due  to  step  command  inputs 
r = i 1 , i2  , *4  . 


w *10  i 


Unfortunately,  id  cannot  be  known  beforehand,  leading  to  an  unavoid- 

d 

able  cut  and  try  procedure  that  will  be  discussed  later. 

This  specific  problem  is  easier,  in  the  above  sense,  if  there  ai 
specific  bounds  on  the  disturbance  response,  for  a class  of  input 
d isturbance;;, when  disturbance  attenuation  imposes  a bigger  feedback 
burden  than  satisfying  command  response  specifications. 
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Figure  VII. 7.  Derivation  of  the  bounds  on  G(jw) 


Deorc  or; 


some  frequencies. 


- Bounds  on  G due  to  disturbances. 

If  one  is  to  derive  G(s)  that  satisfies  the  bounds  on  G due  to 
command  inputs  only,  one  ends  up  with  a final  design  for  which  the  cross- 
over frequency  of  the  equivalent  open  loop  transfer  function  6 [4,80], 

Therefore  for  such  a design,  one  guesses  that  sa  4 rps  and 

min 

■j>d  « 80  rps  for  the  disturbance  response  approximation.  It  is  obvious 
max 

that  those  values  represent  approximate  minimum  bounds  for  both 


“d  . and  "d 
min  n 


because  any  G(jw>)  which  will  also  take  care  of  the 


idditional  disturbance  specifications  can  only  be  more  conservative,  and 

therefore  , can  only  be  biaoer.  So,  one  guesses  some  range  for 

d 

. , | say  (8,400)1  and  derives  corresponding  templates  P°[sl  for 

. 7,4,.  ..,11'  which  are  plotted  on  Fig.  VII. 6. b for  some  w 


52 


'"•orrrrpondinq  bounds  B . (co  ) due  to  disturbances  are  then  derived  and 

d E 

i 'dotted  lines)  in  Fig.  VII. 7.  A new  G can  be  derived,  as  will 

be  explained  later,  and  therefore  one  ends  up  with  a set  of  possible 
valuer  for  . Conceptually  this  is  repeated  until  the  real  range  of 

corresponds  to  the  one  predicted. 

VII. 2. b. 6 Derivation  of  G(jw)  . 

Using  (7.3),  the  final  bounds  B(w)  on  G due  to  the  set  I of 
inputs,  (L.e.  for  a = 1 , 2 , . . . , 12)  are  then  obtained  and  plotted  on 
Fin,  VII. 8 for  some  ui  . The  optimum  G(ju>)  lies  on  its  bound  at  each 
>■>  and  an  approximation  (Fig.  VII. 8)  obtained  is: 

. 14.  1010  (s+1.25) (s+2) (s+120)2  ,, 

G(S)  ~ (7.1 2) 

s (s+5) ( s+60) (s+300) ( s+400) (s  +800s+10  ) 


l-'iuiirc  VI  1.8.  Final  bounds  on  G and  derivation  of  G 
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for  an  excess  of  3 poles  over  zeros. 

The  equivalent  open  loop  transfer  functions  corresponding  to  the 
nominal  plant  Kl=l  , K2=l  are  then  plotted  on  Fig.  VII. 9 for  the 
values  r=  1,2,4  of  the  command  step  inputs.  We  have  (Fig.  VII. 9) 

[20  , 400]  which  is  within  the  predicted  range  [5,400],  implying 
some  overdesign,  and  therefore  one  should  begin  again  the  design  with 
a new  prediction,  say  [10,400].  However,  in  the  present  example, 
removing  the  points  corresponding  to  6 [5,20]  leads  to  new  templates 
(Fig.  VII. 6. b in  dashed  lines)  which,  in  turn,  leave  the  final  bounds 
B (<u ) unchanged.  Therefore  G(jio)  derived  previously  remains  unchanged 
too,  and  the  final  design  is  given  by  (7.12). 

Note  that  at  high  frequency  all  three  loop  transfer  functions 

(Fig.  VII. 9)  are  identical.  This  is  due  to  the  fact  that  pa[s]  is 
independent  of  r at  high  frequency,  as  seen  from  (7.10). 

VII.2.b.e  Derivation  of  the  nonlinear  prefilter. 

Conceptually,  one  can  derive  3 different  prefilters  corresponding 
to  the  3 different  linear  time  invariant  designs  that  were  obtained. 

The  method  presented  in  [HI]  and  in  [A2]  can  then  be  used  to  actually 
build  the  nonlinear  prefilter. 

However,  the  relation  (7.5)  in  T suggests  that  one  can  use  a 
similar  relation  for  F , namely: 

F.(s)  = F. (As)  with  A such  that:  r.(t)  = Ar.(t)  (7.13) 

3 1 3 1 

Indeed  in  a typical  L.T.I.  design,  L(joj)  copes  with  plant  uncertainty: 
therefore,  at  least  for  hiqh  gain  factors,  L/l+L = 1 from  the  command 
input  i>oint  of  view,  and  therefore  T(s)  ± F(s)  . If  the  reference  is 
taken  as  the  unit  step,  then  (7.13)  is  rewritten  as: 

Fr(s)  = F^rs) 


(7.14) 
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G(s)P. (s) 
3 


inf  \ & ..  ( rs ) 

j rr*N  i+G(s)Pa(s) 


Min  T^rs) 


(7.16) 


Taking  into  account  that  the  inequality  sign  > in  place  of 


(7.16)  implies  overdesign,  F^(s)  is  thus  chosen  as  to  satisfy: 


F (s)  = Min  *£  (s)  (7.17)  for  r = 1,2,4 

1 w r,N 

r 

It  is  found: 


F^s) 


6.25  x 60 


(7.18) 


(s+60) (s  +4.75s+6.25) 
which  is  plotted  in  Fig.  VII. 10. 


One  needs  now  to  check  that  for  all  to  , and  for  r- 1,2,4 


sup 

j 


G(s)P . (s) 


F (s) 
r 


l+G(s)P.  (s) 


< Max  T (s) 
" r 


(7.19) 


when  j ranges  over  all  possible  plants,  which,  in  the  present  case, 
is  easily  satisfied.  If  this  were  untrue,  i.e.,  if  3 some  frequency 


range  ft  on  which  (7.19)  is  violated  for  some  r , then  one  must 
a 


use  a smaller  tolerable  variation  A TMju))  over  that  range  and  modify 


the  design  accordingly.  . This  procedure  is  certainly  convergent,  but 
may  lead  to  some  overdesign,  the  price  that  has  to  be  paid  in  return  for 
such  a simple  nonlinear  prefilter. 

VII. 3.  Results. 


VII. 3. a.  Command  inputs. 

As  shown  in  Fig.  VII . 11 . a. & b. , the  nonlinear  feedback  system 
satisfies  the  N.L.T.D.S.  for  |r|  = 1,2,4  . Because  of  the  smoothness 

of  the  nonlinearity  it  can  reasonably  be  expected  that  the  N.L.T.D.S. 
will  be  satisfied  for  all  |r|  6 [1,4]  . For  example,  the  system 

response  to  a step  r = 3 is  plotted  in  Fig.  VII. 12. 

However,  nothing  can  be  said  for  Jr!  >4  and  |r|  < 1 . The 


v- 

. . 


TIME  (SECONDS) 


Figure  VII .15.  Nonlinear  system  step  disturbance  response  for 
different  plant  parameter  conditions.  !d|  = 4 . 
.(a)  K2  = 1 
(b)  K2  = 10  . 
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system  response  to  truncated  ramp  and  sinusoidal  command  inputs  are 
shown  in  Fig.  VII. 13  and  VII. 14.  a & b respectively.  They  are  seen  to 
be  quite  reasonable,  although  the  amplitude  of  the  sine  wave  for 
example,  varies  between  -5  and  +5. 

VII. 3. b.  Disturbance  inputs. 

From  Fig.  VII. 15.  a and  b it  can  be  seen  that  the  specifications 
of  30%  maximum  overshoot  for  the  step  disturbance  responses,  is  satis- 
fied over  the  range  of  plant  uncertainty. 

For  step  disturbance  greater  than  6 , however,  the  system  becomes 
unstable  for  some  plant  parameter  values.  This  result  is  expected 
because  the  resulting  equivalent  linear  plant  representation  predicted 


instability  for  the  design  used. 

Fig.  VII. 16  and  VII. 17. a to  d present  the  system  responses  to 


16  3 


truncated  ramp  and  sinusoidal  disturbances  respectively.  They  are  seen 
to  be  reasonable,  although  the  amplitude  of  the  input  signal  varies 


between  - 5 and  + 5 . (Recall  that  the  design  was  guaranteed  for 


Figure  VII. 17.  Nonlinear  system  response  to  a sinusoidal  disturbance 
d = 5sinugt  for  kl  = 10,  k2  = l.  (a)  ug  = 1 rps, 

(b)  tog  = 10  rps,  (c)  nig  = 100  rps,  (d)  ug  = 1000  rps. 


TIME  (xIO) 

Figure  VTT.'8.  Nonlinear  system  response  (3)  to  a combination  of  a 

1 ! ’ step  command  (1)  and  step  disturbance  (2)  for  the  plant 

condition  Kl  = 10  , K2  = 1 . 


'>■“***■*%  XT 


VII. 3. c. 


Other  system  inputs. 

A combination  of  a step  command  and  a step  disturbance  leads  to 
satisfactory  results,  although  the  superposition  theorem  does  not  hold 
here. 
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Figure  VII. 19.  (a)  white  sensor  noise  (system  input). 

(b)  nonlinear  system  response  (2)  for  K1 = 10  , K2 
to  a step  command  (1)  r * -4  in  presence  of  (a) . 
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Figure  VII. 20.  (a)  Disturbance  input 

(b)  Nonlinear  system  response  (2)  for  K1  = 1,  K2  = 10  to 
a step  command  r = -4  in  presence  of  (a) 

(c)  Nonlinear  system  response  (3)  for  K1  = 10,  K2  = 1 to 

a step  command  r = 1 (1)  and  the  disturbance  (2). 
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VII. 4.  Conclusions . 

The  EPLTI  concept  was  used  successfully  in  the  design  of  a non- 
linear feedback  system  to  achieve  N.L.T.D.S.  It  was  shown  that  a L.T.I. 
compensation  G and  a nonlinear  prefilter  F handle  this  nonlinear  problem 
and  achieve  quantitative  specifications  both  on  the  command  and  on  the 
disturbance  signals. 

Nonlinear  specification  on  the  disturbance  was  not  emphasized 
here  but  this  could  easily  be  handled  by  our  design  technique. 

Applications  of  the  synthesis  procedure  to  L.T.I.  plants  is  straight- 

! 

forward,  leading  to  the  synthesis  of  L.T.I.  loops  with  nonlinear 

! prefilters,  to  ach;eve  N.L.T.D.S. 

t 

The  problem  of  a nonlinear  prefilter  has  widely  opened  and  supolied 

i 

considerable  motivation  in  the  research  area  of  nonlinear  open  loop 
synthesis.  Some  primitive  solutions  are  presented  in  Appendix  A2  , 
where  it  is  shown,  when  possible,  (i.e.,  for  a large  class  of  problems), 
how  to  use  the  same  relationship  on  the  set  {F°}  that  characterizes 
the  set  of  {Ta}  , where  Ta  is  a representative  nominal  T from 
the  set  of  acceptable  transfer  functions  for  r“  . This  approcah 

■ 

t has  the  advantage  of  simplicity  in  both  the  derivation  and  implement- 

M 

ation  of  F , although  it  might  sometimes  lead  to  some  overdesign,  a 
price  that  may  have  to  be  paid. 

4 

In  our  design,  a L.T.I.  G(s)  was  used  which  also  accounts  for 
the  simplicity  of  the  design  procedure.  This  obviously  leads  to  over- 
design  in  the  bandwidth  of  the  effective  L.T.I.  open  loop  transfer 
function  for  some  inputs.  Certainly  a nonlinear  G would  be  much 
better  in  this  sense  and  this  should  definitely  be  one  research  direction 
to  pursue.  Such  a nonlinear  G should  present  a different  "transfer 

k* 


function"  for  each  command  input  and  disturbance  considered. 

In  such  a case,  let  (e)  be  the  set  of  inputs  to  G when 

°1 

ial  = r^  . There  is  a set  because  of  the  set  W , i.e.  each  w in 

W gives  a different  z . It  is  necessary  that  the  nonlinear  G acts 

like  Ga-'-(s)  a fixed  transfer  function  to  the  set  {e}  , like  G°2(s) 

aj_ 

to  the  set  {e}  , etc. 

— “2 

If  disturbances  are  also  to  be  considered,  then  the  problem 
becomes  even  more  difficult. 

Therefore,  one  should  be  aware  that  if  the  problem  of  a nonlinear 
prefilter  is  a tough  one,  the  problem  of  a nonlinear  G is  definitely 
a vm-y  difficult  one  for  which  even  primitive  solutions  do  no  exist. 


Al-1 


A justification  is  presented  here  for  the  use  in  chapters  III  & 

go  2e_tds 

VII  of  a delayed  second  order  model  T (s)  = as  a means 

A 2 __  2 

s +2CUJ  S+GO 
n n 

of  predicting  the  overshoot  in  the  L.T.I.  system  step  response.  Since 
mathematical  equivalence  can  obviously  not  be  made,  some  examples  will 
be  considered.  We  will  show  that  if  Tft(s)  is  used,  the  predicted 
overshoot  is  the  same  as  that  obtained  experimentally.  Using  Equation 
(3.1a)  and  (3.1b),  Fig.  Al.l  is  first  derived  which  gives  the  well 
known  relations  that  exist  between  the  damping  factor  t , the  over  - 

shoot  and  the  peaking  in  |Tft  (jto)|  of  a delayed  second  order  system. 

Peaking  in 
20.  db  of 

18.  1+La' 

16. 

14. 

12. 

10. 

8. 
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4. 

2. 

damping 

.1  .2  .3  .4  .5  .6  .7  .8 

Figure  Al.l  Overshoot  in  the  step  response  and  peaking  of  | L/l+L| 
versus  damping  for  a second  order  system. 
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Consider  the  high  order  open  loop  transfer  function  L^(s) 


(7  zeros  - 12  poles)  used  in  chapter  III,  Fig.  III. 13,  corresponding 


to  minimum  gain  factor  k = l . As  k increases  from  1 to  k = 1000, 

max 


the  corresponding  L (s)  is  obtained  by  shifting  the  0 db  line  in 

1 I * 


Fig.  III. 13  by  the  amount  - 20  log  k . Table  A1.2  gives  the  peaking 


of  |li  / 1+l.  . versus  k . Fig.  Al.l  is  then  used  to  give  the 
1 ' k X'k  db 


TABLE  A1.2 


1 1 

1 Ll.k  1 

■ 1 1 

predicted  overshoot 

! 

actual  overshoot! 

_! 

peaking  in 

1 1+Li,k  ; 

i 

1 

0. 

- 

1 

\ ^ 

2 

0. 

- 

- 

t i 

3 

0. 

- 

! 

5 

.3 

.09 

.09 

: i 

7 

1.35 

.17 

.17 

10 

2.93 

.27 

.26 

20 

6.63 

.46 

.44 

30 

9.1 

.56 

.53 

50 

12.4 

.68 

.64 

70 

14.6 

.73 

.70 

a 

i l 

i 100 

17. 

.78 

.75 

V 

200 

20. 

.86 

.82 

* ' 

i 300  i 

21. 

.86 

.84 

*. 

t* 

5 500 

21. 

.86 

.86 

700  | 

21. 

.86 

.88 

• 

^ 1000  j 

, 

1 

20. 

.86 

.89 

predicted  overshoot  in  Table  A1.2,  which  is  then  compared  with  the 
experimental  values  (Fig.  A1.3.a-c).  Note  the  good  agreement. 


despite  the  shape  of  L^(jw)  fFig.  III. 13).  Indeed,  L^tju)  has 


many  corner  frequencies  in  the  low  frequency  range  and  its  shape  in 
the  high  frequence  range  (w > 600  rps)  is  far  from  resembling  a second 


L-5 


order  system.  Despite  those  two  facts,  it  is  seen  that  the  prediction, 
as  far  as  the  step  response  overshoot  is  concerned,  is  good. 

This  can  be  understood  by  considering  the  root  locus  of  1+L,  , 

as  k varies,  which  is  plotted  on  Fig.  Al.4.  At  small  k , the 
dominant  pair  of  closed  loop  poles  are  located  on  arc  2 (the  poles  of 
arc  1 are  very  close  to  the  zeros,  even  at  k = 1,  so  their  effect  can 
be  neglected,  while  those  of  arc  3 are  far, relative  to  those  on  arc  21 . 
The  pair  located  on  arc  3 becomes  dominant  at  intermediate  and  high 

gain  values  of  k , the  poles  of  arc  2 being  then  very  close  to  the 
zeros.  The  effect  of  the  two  far-off  pole  pairs  on  arcs  4 and  5 can 
be  neglected,  pee  [H7]).  Therefore  a second  order  model  for  T(s)  can 
be  used  here  and  this  will  remain  true  in  general  for  systems  with  large 
uncertainty,  because  of  the  universal  character  of  the  resulting 
optimal  loop  functions  ( (H2 ) , [S2])  for  plants  with  large  high  frequency 
gain  factor  uncertainty,  to  which  this  work  is  devoted.  Note  that  it 
is  implicitely  assumed  that  the  two  following  cases  do  not  appear,  as 


examples  of  open  loop  Corresponding  system  responses 

transfer  functions.  c.  to  a step  disturbance. 
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discussed  below. 

1st  case:  ! L/L+ 1 i qualitatively  as  shown  (see  L ) in  Fig.  A1.5.a 
implying  the  system  step  disturbance  response  c^(t)  of  Fig.  A1.5.b. 

In  such  a case  it  is  obvious  that  the  overshoot  cannot  be  characterized 
by  the  above  formula.  Note  that  for  such  a system  we  would  have  a root 
locus  for  1+L  qualitatively  as  shown  on  Fig.  A1.5.C.  Such  a phenomena 


f.  ‘or  in  his  far-off  pole-zero  package  [F2]  in  order  to  descrease  'l| 
v y fast.  Suffice  it  to  say  that  (1)  the  L.T.I.  system  response 
depicted  in  Fig.  Al.S.b.  is  certainly  not  desirable. in  general,  (2)  such 
minor  peakings  in  |l/1+l|  should  be  avoided,  especially  if  there  are 
bending  modes  present,  (3)  minimization  of  k does  not  necessarily 
mean  minimization  of  the  sensor's  noise  effect  at  the  plant  input,  sc 
a reasonable  compromise  is  to  use  characteristics  like  L^fju) 

(Fig.  A1.5.a)  rather  than  ( jou) 

2nd  case:  L(jw)  is  qualitatively  as  shown  in  Fig.  Al.6.a  . Such 

L - shaping  can  occur  in  designs  where  the  plant  template  at 
becomes  "fatter"  than  at  , causing  the  bound  on  to  be 

significantly  higher  than  on  I.(jruj)  . The  effect  of  such  a loop  in 


open  loop  transfer 
function  L(juj) 


Fig.  Al.6.a  gives  rise  to  two  peaks  in  |l/1+l|  as  shown  in  Fig.  A1.6.b 


and  is  usually  achieved  by  a pair  of  complex  zeros  and  complex  poles  in 


L(jw)  as  shown  in  Fig.  A1.6.C.  It  should  be  recalled  that  such  pheno 


mena  is  less  likely  at  high  frequency  because  P(s)  -*■  k/  e 


at  such  frequencies,  at  least,  the  poles  of  L/l+L  on  arc  1 are  practic 


ally  cancelled  out  by  the  complex  zeros  of  L 


and  therefore  the  complex 


pair  lying  on  the  arc  2 can  be  considered  as  dominant  and  a second  order 


approximation  is  then  legitimated.  In  those  cases  where,  at  low 


approximation,  a more  complex  model  must  be  derived  unless  overdesign 


Figure  A1.6.C  Root  locus  of  1+L 
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The  problem  considered  here  is  to  design  a system  which  achieves 
N distinct  output  signals  u,  (t)  . . .u„(t)  when  it  is  excited  resDect- 

1 N 

ively,  by  N different  signals  R (t)...R^(t) 

A2 . 1 The  Input  signals  are  "qualitatively  identical". 

By  "qualitatively  similar"  signal  set,  it  is  meant  that  : 

V i,j  3 K € ft  such  that  R^(t)  = KR^  (t)  (A2.1),V  t . This  case  was  consid- 
ered in  chapter  VII  where  the  input  signals  were  steps  of  amplitude  r . 

Because  some  empirical  formula  u (t)  = ru, (t/r)  could  then  be  found 

r 1 

on  the  output  signals,  there  existed  a very  simple  prefilter  whose 
transfer  function  had  the  simple  expression: 

Vs)  i 

F (s)  = F (rs)  with  F (s)  A --  - ■ — and  R,  (s)  = — • This 

r 1 1 R^  (s ) 1 S 

concept  will  now  be  extended. 


As  another  example,  assume  that  the  set  of  inputs  consists  of  impulse 
R.  (t)  = r.6(t-f)  , r.  € R , which  are  to  produce  the  set  of  outcuts 

ill 

shown  on  Fig.  A2.1.a  where  V r^  the  maximum  value  is  y occuring  at 


a time  t inversely  proportional  to  the  amplitude  of  the  input  signal 

— £ 

Suppose  that:  Uf(t)  = Brt  e (A2.2.a),  for  which 

t = ~ and  ^ ~ • Therefore  U (s)  = — — — — and  as  R(s)  =r 

r (s+ra) 

6 

we  have:  F (s)  — (A2.2.b). 

r (s+ra) 

More  generally,  if  u (t)  = g(r)t  e then  F (s)  = /r  _ 

r r (s+f  (r) ) 2 


This  transfer  function  is  easily  implemented  (Fig.  A2.1.b)  on  an  anolog 

computer,  the  amplitude  r being  obtained  in  this  case,  by  integrating 

t 

the  input,  i.e.,  r = / R(OdC 
0 

As  a third  illustration,  let  the  set  {u}  of  desired  outputs  in 
response  to  step-inputs  be  as  shown  on  Fig.  A2.2.a. 


Figure  A2.2.a.  Time  domain  sooeif ications  as  a function  of  the 
amplitude  r. 


SSliAL'a’ 


which  is  also  easily  implemented  as  shown  in  Fig.  A2.2.b 


Figure  A2.2.b  Block  diagram 


These  simple  examples  suggest  that  when  there  exists  a relation 


u (t)  = f(u, (t),r)  then 


[F  (s),r]  which  can  be  written  in  the  general  form 


F (s)  is  easily  implemented,  as  shown  in  Fig.A2.3  provided  the 


amplitude  r is  known  , i.e.  F is  "tuned"  to  this  input  set,  so  the 


inputs  must  be  restricted  to  this  set 


Up  to  now  we  only  considered  steps  or  impulses  for  which  the 


ermi nation  of  r was  straightforward.  In  the  general  case,  let 


the  signal  R(t)  of  amplitude  r be  described  by  the  differential 


The  left-hand  side  of  (A2.6)  is  easily  obtained,  as  it  only  involves 

integrations  of  the  input  signals,  and  so  is  yMt)  too.  Therefore 

V t > 0 , r can  be  estimated  and  its  value  is  then  used  to  compute 

all  7 ( r ) , p (r)  of  (A2.4).  It  is  important  to  recall  that  the 
i 1 

above  is  restricted  to  'qualitatively  similar'  input  sets. 

A 2 . 2 The  input  signals  are  linearly  independent. 

Here,  consider  the  more  general  case  of  n independent  input 
signals  R^  , . . . , R^  , i.e., 

X,R,  (t)  + X R,(t)  + ...  +X  R (t)  =0  **  X = X = ...»  X =0  . (A2.7) 

11  22  nn  12  r 

Consider  the  system  of  Fig.  A2.4  where  the  sampling  period  T is 
supposed  to  be  very  small.  We  have:  [of  the  modified  Z transform  in 
chapter  IV] 

U*  (z)  = A*  (z)  R*  (zN)  + A*  (2)R*(zN,^jp)  + ...  +A‘  (z)R*(zN,^) 

1 2 N K+l  N 


(A2.8) 


Such  an  equation  can  be  written  for  each  of  the  N pairs  (R^,IT)  giving 
N equations  with  N unknowns  (of  course,  if  N is  small  there  isn't 
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Figure  A7.4.  System  conf iouration 


enough  flexibility  for  good  matching  - but  one  can  use  (N  + M)  terms 


and  choose  M of  these  arbitrarily  to  suit  the  approximation,  solving 


Using  (A2.7)  the  rows  are  independent  and  (A2.9)  has  a unique  solution 


which  then  solves  the  problem  for  the  case  considered 
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APPENDIX  .A3:  DETAILS  OF  THp  L.T.V. 


A3.1  Derivation  of  from  h2<t,x) 

In  any  realistic  L.T.V.  system,  the  closed  loop  impulse  response 
h2(t,f)  is  the  solution  of  a linear  differential  equation  of  order  N 
with  variable  coefficients,  where  the  input  x to  be  considered  is  a 
unit  impulse  occuring  at  t — T This  implies  [S5]  that 


h2(t,x)  = A ( t ) • B (t)  Tu  (t-x) 


(A3 . 1) 


with  A ( * ) and  B(*)  two  N-dimensional  vectors.  By  anology  with 


L.T.T.  systems,  if  h2<t,x)  is  of  order  N , so  is  the  open  loop 


impulse  response  i (t,x)  , and  therefore 

^ a (t)  • 6 (t)Tu  (t-r) 


(A3. 2) 


with  a(')  and  6 ( * ) two  N-dimensional  vectors. 

and  h2  are  related  through  (6.11),  i.e.. 


/ '»  ,(t  .O  + «(t-0)h  (C,T)dC  = H-(t,x) 


(A3. 3) 


Noting  that  x < £ < t due  to  the  causality  of  h2  and 
using  (A3.1),  (A3. 2)  in  (A3. 3)  gives: 


and 


T<t  A(t)-B(x)  +a(t)  / 8(()  A(()d(B(x)  =a(t)*8(T)  (A3. 4) 

x 


Let 


Aft)  = f R(0T A(c)dc  an  N»N  array 

0 


(A3.f) 


Then  (A  t. 4)  1 1<  •comes:  ( A(t)  + a(t)  A (t)  ] •B(t)T  = a(t)*[B(x)T+A(x)B(x)Tl 


and  identifying  terms  in  t and  in  x gives: 
f .*  (t ) * A (t ) [ 1 — A (t)  ) _1 


(A3. 6a) 


B(t)  = B(t)  [ 1 — A (t)  ] ' 


(A  i.f  b) 


Using  (A3. 6b)  in  (A3. 5)  and  differentiating  with  respect  to  time  gives 

A (t)  = [ 1-A  (t) ] B (t) TA (t)  (A3. 7) 

A particular  solution  is  A (t)  = 1 

In  section  VI. 2. b it  was  found  that 

h (t,T)  =D  (t)ir‘1(t)*D  (T)Tu(t-x)  with  * (t)  = U u + / D _ (t ) TD.  (C ) dC 

0 2 2 

(recall  6.19)  so,  from  (A3.1),  A (t)  = D2  ( t ) tt_1  (t)  and  B(t)=D2(t) 

The  homogenous  part  of  (A3. 7)  becomes: 

A (t)  + A (t)  [D2  (t)  rD2  (t)-TT-1  (t)  ] = A (t)  + A(t)ff  (t)ir_1(t>  = 0 

which  implies  that  the  homogenous  solution  is: 

Ap(t)  = ^C*tt  '*'(t)  where  K is  a non-singular  matrix  , 

independent  of  time.  Therefore,  (A3. 7)  gives  A(t)  = 1 + IK  it'V) 

and  replacing  in  (A3. 6) : 

r a (t)  = A (t)  (-JCi i'1(t))'1 
6 (t)  = B (t)  (-  tt_1  (t) ) T 
Therefore  from  A3. 2,  A3.1, 

*2(t,T)  = D2  (t)ir-1  (t)  ( — tt  (t)^C_1)  (-  tt-1  (t)  )D2  (t  )Tu  (t-r) 

= D2  ( t ) tt  1 (t)  *D2  (t)Tu  ( t— T ) (A3. 8) 

A3. 2 Derivation  of  h(t,f)  from  ft(t,T) 

Let 

h(t,t)  = (t)  *B  (t)  Tu  (t-T)  (A3. 9) 

and 

e.  ( t , T ) = (t)  •(-t1(T)Tu(t-T)  (A3. 10) 

where  in  the  present  case  A^  and  B^  are  the  unknowns,  while 
<i^  , lJ.j  are  given.  (A3. 3)  to  (A3. 6a  R b)  can  then  be  rederived  with 


- ■ 1 


V j 

* 'I 

{ 


A3. 3 


subscript  1 

Therefore,  (A3. 6a)  and  (A3.6o)  are  rewritten  as: 


A^  (t)  = ai(t)  !U  -A  (t)  1 


B,  (t)T  = tD  -A  (t)  ] _16  <t)T 


(A3. 11a) 
(A3. lib) 


and  (A3. 5)  becomes  A (t)  £ / B^C)  A^OdC  = / 6]. ( ^ a1  (?)  1 ^ -A <C)1  dC 


or  by  differentiating: 

, T 


(A3. 12a) 


A (t)  = B (t)  c^ft)  ('ll- A (t)  ] 

A particular  solution  is  A(t)  = U 

The  open  loop  impulse  response  Mt,t)  associated  with  P € 'P 


is  (Eq.  6.9.b)  i = p * p^1  • %2  = jp  • l2  at  p = k . So 


f.(t,T)  = Ai2(t,T)  = AD2(t)’TT(T)  1*D2  (x)Tu(t-T)  by  using  (A3. 8)  and 


defining  X = kA0  • By  identification  with  (A3. 10)  , a^(t)  = XD_,(t) 


T -1  T 

and  (t)  = 11  (t)  D2  (t)  and  therefore  the  homogenous  part  of 


(A3, 12a)  becomes: 


A (t)  + 61(t)Ta1(t)A(t)  = A (t)  + Xw(t)_1D2(t)TD2(t)A(t) 


-1  ■ 


= A(t)  +Xir(t)  ir(t)A(t)  = 0 


(A3. 12b) 


from  (6. 19) 


(A3. 12b)  has  no  closed-form  solution,  in  general,  so  let  A^  (t) 


denote  the  solution  to  (A3. 12b).  Then:  A (t)  = +A^  (t)  §C.  where 

flC  is  a nonsingular  matrix.  From  (A3. 11a  & b)  and  using  the  above 
definitions  of  a,  and  R 


'1  pl  ' 

At(t)  = XD2(t)  (-Ax(t)  IK  ) 

^ BjJt)7  = - flC'1A1(t)"1  Tl(t)'1D2(t)T  . 

Ilonre,  from  (A3. 9),  the  closed  loop  impulse  response  associated  with 


* -L.  . . 


where  A,  is  such  that 


In  the  special  case  where  the  input  R is  deterministic,  or  when 


D (t)  is  a one  dimensional  vector  m(t)  becomes  a number  and  so  is 


A. (t)  (instead  of  a matrix).  Then  the  solution  to  (A3. 13a)  is 


and  (A3. 13)  becomes,  in  this  special  case 


A3. 3 Derivation  of  the  system  response  c(t)  associated  with  p-k 


Let  c(t)  be  the  system  response  to  a command  r 


when 


Prom  Fig.  VI. 3 


For  the  sake  of  simplicity,  let  N=1  ; this  imples  that  the  closed 


loop  impulse  response  of  the  bottom  part  of  Fig.  VI. j is  given  by 


(A3. 13b) . Under  the  above  assumption,  (Fig.  VI. 3) 


Ac (t)  = D (t)  - / h(t,C)D(i;)dt 


(recalling  the  definition  of  D from  (6.9)) 


Replacing  (A3. 13b)  in  (A3. 15)  gives: 

D (t)  t 

Ac  (t)  •*  D(t) / XD 2 (C) (C) 


D(e)d  C 


(A3. 16) 


6 00 


(A3. 16)  can  bo  f .rittei.  as: 


Ac(t)  = D (t) yr-  / Xtt(5)X'1d  (C)2dC 

ir(t)  0 

recalling  (6.19)  we  have  D2<C)2  = ir  (O  so 


Ac (t)  = D(t) 


-"“’[(Tiff)1]”-  D'‘>fe)‘ 


(A3. 17) 


(A3. 18) 


Using  (A3. 15a),  (A3. 14)  becomes 

c,l>  - V"  * ^ (?rSf) 


So: 


(A3. 19) 


A3. 4 Effect  of  white  noise  at  plant  input. 

Under  the  above  assumption  that  N = 1 , we  can  use  (6.1.b)  and 

(6. 12.d)  to  obtain  the  mean  square  value  of  the  noise  at  the  plant 
input  when  P = P2  , namely: 


P.I. 


2 0 0 

By  anology  with  (6.12.d), 


0 = - i • h 
k 


1_  t t 

2 //h2(  t,Vh2(t,(2)^(V9d(ldC2 


(A3. 20) 


(A3. 21) 


when  the  plant  is  any  P = kC  f , associated  with  the  closed  loop 
’.'-pulse  response  h of  (A3. 13b).  Therefore,  the  mean  square  value  of  the 
ttni-.e  at  the  plant  input,  becomes  in  this  case:  (from  (6.1.b)  and  (A3. 21)) 
t t 
k2  0 0 


2 i t t 

Vt.  = ~2  / ^ h(t'Cl,h(t'l;2,YNN(!:i'C2,dCl<l!:2  (A3. 22) 


i.>'l  the  sensor  noise  be  a white  stationary  process  of  strength  c2  , 


l .e.  , 


Ywv*r,?'r,p^  = ,0  (A3. 22)  becomes: 


a2  t 


j = / h(t,C)2dC  • Using  (A3. 13b), 


x21d 

k 0 ir(t) 


2 2 
0„  D (t)  * t 
N 2 r 


Op  _ = — j / tt(0  TT(t)d;  , because  (6.19)  ir(t)  =D  (t)‘ 

’ k^  w(t)  0 

and  (6.21a)  X = k/k. 


Finally: 


vt)2  r i .(o )2>-1 


«p.i.  (tl  - 


2 2X-1  [ w(t) 


if  X / i, 


(A3. 23) 


f»  h!l!l 

„2  *(t) 


\ if  (0)  / 


if  X = ij 


As  X = k/k2  » for  any  t = tQ  , op  (tQ)  is  a function  of  k 

Let  u = 2X  , then  for  X / <j  , i.e.  u^l  , 

2 ,«.»  Zll  r,  r,-:2  -L-  f_L_  ^(O)11"1  \ 

°P.I.(t)  ' .2  °2(t)  u-1  \ir(t)  J 


Therefore, 


3°P.I.  = _2_  30P.I.  _ _2_  2 dJ  (u) 

3k  'k2  3U  'k2k2  2 du 


j(u)  = — — lS21 'j  aH  derivatives  being  taken  at 

u-1  \w(t)  ...u  J 


fixed  t^ 


dJ  _ 1 

du  2 

III.  I 1 


(u-1)2  I w(tQ)U  "V  W(0) 


(A3. 23a) 


which  can  be  written  as: 


dJ 

1 

. u-1 

T (0) 

f 1 -( 

"Vf1 

1 1 On 

r(v 

u-1 
'l  1 

du 

(u-1)2 

w(t0)U 

L 1 V 

*(o)  ; 

\ »(0) 

) J 

Since  ir  (t 

A tQ 

) § ir(0)  + / D 

0 

(02dC 

(recall  6.19),  we  have 

*(t  ) 

X=  .(0)  51 

Noting  that 

1 - X 

U 1 + log 

xU  1 < 0 for  all  x > 0 

, we 

conclude  that 

dJ 

< o f n y 

all  ii 

3 PI  pi  V IS  o * 

refore 

3<I.  . . 

all  V 

, ^ u tor 

du 

ai.1  U 

ana  tne. 

3k  ° for 

all  K 

f wnicn 

imples  that 

2 

°P.I. 

(tQ,k)  is  a decreasing  function 

of  k 

for  given  tQ 

Note  ■ 

that  as 

u-1  , 

(i.e. , 

X - S ) then  u 

= 1+e 

(with  t 

positive?  or  negative)  , and 


1 tt  (0) e 

I 1 

- XE  + e log  X 1 

2 , 1+e 

c tt (tQ) 

1 t(0)£ 

I 1 

2 

/ 1 J.P  l!  p\l#  /On  w \ 

2 , ,1+E 

e if  (t0) 

^ Ite  Jinx— c \ x>nx j 

it  (0)E  (£nx)  2 

•V1*' 

after  using  a serie 

) . i 

■ 

tnf (V  ^ I" 

) TMt0) 

\ tt(0)  / J 

3a;  T k„ 

— ^ — arc  both  continuous  at  X = S (k  = — ) . 

This  explains  the  shape  of  curves  A in  Figs.  VI. 5. a & b. 

A 3 . 5 Solution  to  a certain  filter  problem. 

The  problem  in  (6.40)  is  to  find  hjtt.T)  such  that: 

Min  / (Ac2(t)  + W(t)o2  (t))dt  , i.e.,  if  we  take  P»  — 
h2  0 P*X*  s 

T / t t 3h  v 

Min  / f (D_,  (t)  - J h (t,C)D  (OdC)  + U / ~ <t,C>  d;  dt  (A3. 24) 

h n ' - n i - n ““  / 


h2  (, 


with  u 


o w(t) 

N 


( W(t)  is  assumed  constant), 


Let  us  use  a variational  argument,  i.e.,  h2  = h2  0 + ^ ' w*th  ^ an 


arbitrary  function  satisfying  f(t,t)  =0  . 

T 


T r t 

I = / [(D2(t)  - / (h2  0 (t , C ) + f (t,C)D2(i;)d5)2  + 


t r 3h_ 

M / I 
0 

T r 


3f 


V 


i i r~  (t*c)  + ?t (t>C)  j dc 


dt  = IQ  + 6l  with 


<5l  = - 2 / (D  (t)  - / h (t.C)D.(OdC)  / f (t,C)D,(C)dC  + 

n >■  * * r,  *■ 


t 9h  t 

+ 2 U | — 2-,-°  (t.O  ||  (t.C)dCj  dt  + 0(f2)  where  0 (f2) 


indicates  a positive  function  of  f(t,C)  , which  is  zero  when  f=0  . 

A necessary  condition  for  an  optimum  is  therefore  that  6i  = 0 


We  can  write: 


t 3h 


T 3h„ 


j = j dt  / -Jhl  (t,o  i|  (t,C)dC=/  dc  / (t,C)  |f  (t,C) 


3t 


3 L 


0 0"  ” 0 C 

by  permuting  the  order  of  integration.  Integrating  by  parts. 


3t 


dt 


J * / dr. 


3h 


jf--  (t,o  f ( t , c ) 


t=T 


t=r 


- / dC  / 


T 3 h 


2,0 


(t,C)f (t,;)dt 


C 3t 


T 3h 


t 3 h„ 


/ — “ (T.Of  (T,r.)dC  - J dt  / 


0 


(t,C)  f (t,c)dc 


A3. 9 


So, 


61  = 


r T t r t ^ 

2 [ / dt  / [ (D2(t)  * / h2,0(t'c)D2(OdOD2(i;)  +u  ItS.  (t,o] 


1 T 3h 

,OdC  ♦ 2 u / (T,C)f  <T»C>dC 


f (t 


As  61-0  , it  implies  that  the  optimal  solution  is  given  by: 

. 2 


( 8 h2<t.T)  r t , 

2 + [°2(t)  ‘ / h2  (t.ODjdi;  J D2(t)  =0  Vt>T 


at 


h2(t,t)  = 0 


V t 


(A3. 25) 


V 


3h_ 

i 

3t 


(T ,c)  = 0 , ¥ t 


This  proves  that  when  P(s)  =|  , a solution  to  (6.40)  exists.  This 

solution  is  not  in  closed  form,  even  when  in  this  highly  idealized  case 

of  r>2(0  constant.  However,  numerical  solutions  exist. 

Note  that  when  P(s)  =—  , -i  in  (A3. 24)  is  then  replaced  by 

s 

j)  n ^ 2 2j\ 

— h7(t,c)  and  — in  (A3. 25)  is  replaced  by  ~ , while 


3t 


at 

3h2 

(T,c)  = 0 is  replaced  by 


3mh, 


3t 


(T,a  =0  , V C , V m = 1,2 n 


A3. 6 Obtention  of  the  differential  equation,  knowing  h2(t,t) 
Let  (6.54) 

D2(t)(D^(t)  -dJ(t)] 

h2(t'T)  " 2 i 5 u(t-T)  2 g(t)  (D,(t)  -D,(T))u(t-T) 

uNM(t)*n'(t)  /2  2 2 

h2(f,,)  ’ g h2(t'T)  f q(t)02(t)u(t-T)  , from  (6.48) 

h2(t,x)  = (~)  h?  (t,r)  + ~ h2(t,i)  + (qD,)’u(t-T)  + gD^f(t-t) 


A3. 10 

Let  y (t)  = h (t,T  > , x (t)  = o(t-T  ) , and  we  get: 


• . * (go,) ' . 

y = la)  y t 1 y * — (v  - 2 


g ' g 

so  the  D.E.  is: 


(y  - | y)  + gD2x 


••  • f D 

y + y ^ - 


D„  • \ / • • D \ 

. JL. 2 S.)*y(  2(2)2  ♦ 2 _2  . 2 ) . 

D,  g / y \ g g o2  g / 


gD2x 


As 


g = d2  - °2g  and  | = <5T)  ' D2q  + (f>  ( 57  ‘ 2gD2  ) 


D.  , D 

- 2 2 = 2D^g  - 3 — 


and 


2 ( 


^ - ! = 2 (&2  ♦ (dV  - 2 ^ D^gV  &2-?1  D 


g 02  g 


2g  ; -D2 


D,  2 1 


- ‘57’ ' ‘ °25  ‘(s?  ‘ 2’D0(  57  - DS’) 


" V ‘ ; 5J  d2’  ' 5J  * 3 '5J>2 


and  therefore  h2(t,x)  is  associated  with  the  D.E.: 


+ y ( 2D^g  - 3 -2-  ) + y ( D. 
2 


2 1 


,°2  2\ 


y + yv2D2g_  3 57 ; ^ y v °2g  - 2 °2g  - 57 + 3<^  } = gD-x 

(A3. 26) 


A3. 7 Derivation  of  the  Euler-Lagrange  Equation  for  a special  case. 


We  want  to 
T r/  t 


2jr  1 


Min  / j ( 0 (t ) - / h (t,c)D  (?)dr,  j W(t)+u  / [ ^ h,  (t,C)  +a_h_  (t  ,C)  1 dC 
h2  0 LV  2 0 2 2 l 0 ’fc  2 2 2 J 


dt 


1 


u / [ A(t>2(0(t)-e(cx)2+Act)2e(t)2+2A(t)A(t)e  (t)  o (t) -e  (O)+a2A(t)2(9(t)-0(O)2 
o L i 


+ 2a2A(t)A(t)  (9(t)-0(C))2+2a2A(t)29(t)  (9  (t)  -9(0)  ] dC 

. t 2 t 

Let  i)i(t)  = / (9  (t)  -0(C))  dC  , n (t)  = / (9(t)  -9(C))dc 
0 0 

/ . 2 2*2  ••  22  • 2* 
So  = u y A ^ + A t9  + 2AA9n  + o^A  V + 2a2AA\(i  + 2a2A  0n 


(A3. 27) 


Lot  J(t>  = J (0  (t)  - 0 (C))D  (C)dC  , then  (A3. 27)  is: 

0 

t r 2 1 T 

Min  / (D  (t)  - A (t)  J ( t ))  W(t)  + i (t)  dt  = Min  / F(t)dt 

J an 


The  Euler  equation  is  then  derived  as: 


3F  _ _d_  ^ 3F  \ 

m = dt  v 3-  y 


(A3. 28) 


which  is: 


• 2 2 • 

J(D2-AJlW  + u (At0  + A9n  +a2<.i(i  to^AY  + 2a2A6nl 


= u — lAi|»  + A9n  + a A\Ji] 

clt  d 

u ^ Ki()  + aJ>  + A(9n  + a2^)  + A(0n  + 0n 


(A3. 29)  is  rewritten  as: 


0n  + “2^) 


(A3. 29) 


•••••  • •••••  • 2 • 2 2W 

Mi  + A (U>  + 0n  + - a2ij<  - 9n)  + A(0n  + 9n  + - 2a29n  - c^ili  - t9  - J — ) 


After  simplifications:  (n = 9t  , i = 28n  ) 


2 

AiJ:  + A'Jt  + A ( 0n  - a2ij 1 - ) 

2 u 


D JW  t 

— or:  0 (t)=  J etc)  dC 

0 


A(t02  - 2O01 +u>)  +2A9(te-91)  + A(9(t9  - 9X)  - a2(t92  - 29©1  +(P) 


W ( ) ^ 2 °2W 

- -■  y J(n(t)  -0(C))D,(C)dCy  =~  . 


L 
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A 3 . 8 Design  details  for  both  the  LTV  and  LTI  systems. 
A3. 8. a LTV  system 


The  L.T.V.  closed  loop  impulse  response  h2<t,T)  of  (6.27)  is 


associated  with  the  differential  equation: 


< it  2 

[h  i : y + t 


V* 


\ 2 f 2 

Jy  = D AX  where  t(t)  = u + J D (?)  dC 


-1 


Therefore  the  open-loop  impulse  response  l = h (1-h  ) is 


characterized  by: 

in 


: y + fe-57-^r)y  = D2/1TX 


-1 


so  g = P2  ’i2  = n2/k2 


and : 

lg]: 


y + 


. ^2  °l\  _n]r 

\"  D * J y ' 


A C1(S)  -1 

The  overall  transfer  function  t^  is  such  that  T^  (s)  = 


for  c^(t)  = 1-e 


-t 


? 1 

y + 

y = 

X 

N 

t-  \ 

since 

the  prefilter 

S- j 

LI 

J 

f = h 

Hence 

i1*^ 

f i 
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# 

°2 

4 

r • 
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v* 
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/ 7T 

°2 
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I? 
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k2 

Ifl  : 
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y “ 

*1 
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A3 . 8.b  LTI  Systems 

A two  degree  of  freedom  structure  (Fig.  A3. 2)  is  considered. 

The  LTI  system  roughly  achieves  (as  seen  in  Fig.  VI. 6. b) , for 
P = k € [1,100]  , the  T.D.S.  shown  in  Fig.  VI. 6. a for  the  LTV  design. 

The  LTI  synthesis  procedure  given  in  [S2,  H2)  is  used  with 

. 5.69  1Q6(s+8.)  (s+1.57)  (s2+2  ■ 34s+l . 6) 

_ - " g "2  5 

(s  +1.65+1.1) (s  +4s+8) (s  +60s+5500. ) (s  +63s+6000.) 

, . 3396.  (s+1.5)  (s+1.7)  (s+4.67) 

1 2 5 

(S+.36) (s+226. ) (s  +1.32s+.67) (s  +4s+8) 

rf  the  plant  P = k/s  with  k € [1,100]  is  considered  and  the 
T.D.S.  are  those  depicted  in  Fig.  VI. 8. c,  following  [S2,  H2] , it  is 
found  that 

G(s)  11.43  106(s+.6)  (s+15.4)  (s2+l ■ 5s+ . 72) 

(s+.8) (s+34 . ) (s2+2. 3s+1.8) (s2+63. s+6600) (s2+67s+7400. ) 

p(s)  = 290.  (s+.55)  (s+1.43) 

(s2+s+1.5)  (s+1) 2 (s+31. ) 


O-^O  N 


!•' igure  A 1.2  Two  degree  freedom  structure 


[Bl]  BODE,  H.W.,  Network  Analysis  and  Feedback  Amplifier  Design. 

Van  Nostrand,  N.Y.,  (1945). 

[B2]  BOOTON,  R.C.,  An  Optimization  theory  for  time-varying  linear 
systems  with  nonstationary  statistical  inputs.  Proc.  IRE, 
p.  977,  981,  1952. 

[Cl]  CLEGG,  J.C.,  A Nonlinear  integrator  for  servomechanisms,  Trans. 
A.I.E.E.  Part  II,  Appl.  Ind.,  1958 

[C2]  CHESNUT,  H.,  Obstacles  to  progress  in  Non-linear  control. 

Trans.  IRE,  AC-5,  p.  59,  1958. 

[C3]  COURANT , R.  and  HILBERT,  D. , Methods  of  Mathematical  Physics. 
Intarscience  Publishers,  N.Y.,  1966. 

[C4]  COCHRAN,  J.,  Analysis  of  Linear  Integral  Equation.  MacGraw- 
Hill,  1972. 

[C5]  CES  A ill , L.  , Asymptotic  Behaviour  and  Stability  Problems  in 

Ordinary  Differential  Equations.  Springer-Verlag,  Berlin, 
1959. 

[Dl]  DAVENPORT,  W.  and  ROOT,  W. , Random  signal  and  Noise.  MacGraw- 
Hill , 1958. 

[El]  EMELYANOV,  S.,  ULANOV,  G.,  VIKTOROVA,  V.  and  AGAFONOV,  V., 

Synthesis  of  a variable  structure  system  for  automatic  control 
of  thickness  of  hot-rolled  metal.  Remote  & Control,  p.  262, 
1970. 

[FI]  FLEISMEK,  P.,  Design  of  passive  adaptive  linear  feedback 
systems  with  varying  plants.  Trans.  IRE,  AC-7,  1962. 

[Gl|  GIBSON,  J . , Nonlinear  Automatic  Control.  McGraw-Hill,  1963. 

[G2]  C.ELB.  A.  and  VANDER  VELDE,  W.  , Multiple-input  describing 

functions  and  nonlinear  system  design.  McGraw-Hill,  1968. 

|C3]  GANTMACHER,  F. , Application  of  the  theory  of  matrics.  Inter- 


im 


R2 


[fill 

[H2] 

[H31 

[H4] 

[H51 

[H6] 

[117] 

[H8] 

[119] 

IM10I 

MU 

[Kll 


science  Publishers,  1959. 

HOROWITZ,  I.,  Synthesis  of  feedback  systems  with  non-linear 

time  varying  uncertain  plants  to  satify  quantitative  perfor- 
mance specifications.  Proc.  I.E.E.E.,  Vol.  64,  No.  1, 
p.  123-130,  1976. 

HOROWITZ,  I.  and  SIDI,  M. , Synthesis  of  feedback  systems  with 
large  plant  ignorance  for  prescribed  time  domain  tolerances. 
Int.  J.  Control,  Vol.  16,  p.  287,  1976. 

HOROWITZ,  I.,  Optimum  loop  transfer  function  in  single  loop 
minimum  phase  feedback  systems.  Int.  J.  Control,  Vol.  18, 
p.  97,  1973. 

HOROWITZ,  I.  and  KRISHNAN,  K. , Synthesis  of  nonlinear  feedback 
system  with  significant  plant  ignorance  for  prescribed 
system  tolerances.  Int.  J.  Control,  Vol.  19,  p.  689,  1974. 

HOROWITZ,  I.  and  ROSENBAUM,  P.,  Non-linear  design  for  cost  of 
feedback  reduction  in  systems  with  large  parameter 
uncertainty.  Int.  J.  Control,  Vol.  21,  p.  977,  1975. 

HOROWITZ,  I.,  ROSENBAUM,  P.,  SIDI,  M.  and  WU,  M. , A synthesis 
theory  for  non-linear  systems  with  plant  uncertainty.  6th 
Triennial  World  Congress  IFAC.  Part  IC,  37.2,  1975. 

HOROWITZ,  I.,  Synthesis  of  feedback  systems.  Academic  Press, 
1963. 

HOROWITZ,  I.,  A synthesis  theory  for  linear  time  varying  feed- 
back systems  with  plant  uncertainty.  Trans.  I.E.E.E.,  AC-20, 
p.  454,  August  1975. 

HOROWITZ,  I.  and  SIDI,  M. , Optimum  Synthesis  of  non-minimum 
phase  feedback  systems  with  plant  uncertainty.  Int.  J. 
Control.  To  appear. 

IIOI.ZMANN,  E.O.,  Non-linearity  in  process  systems.  Trans.  IRE, 
AC-5,  p.  61-64,  1958. 

JURY,  U. , Theory  and  application  of  the  z-transofrm  method. 

Wiley  R Sons,  N.Y.,  1964. 

KHARAROV,  V.S.,  Optimal  control  of  variable  structure  systems. 


Remote  & Control,  p.  1603-1609,  1970. 

[K2]  KAILATH,  T. , A view  of  three  decades  of  linear  filtering  theory.  1 

Trans.  I.E.E.E.,  XT-20,  p.  146,  March  1974. 

[K3]  KAPLAN,  W. , Operational  methods  for  linear  systems.  Wesley,  1962. 

[LI]  LINDORFF,  D. , Theory  of  sampled  data  control  systems.  Wiley  S 
Sons,  1965. 

[L2]  LANIN,  J.  and  BATTIN,  R. , Random  processes  in  automatic  control. 

McGraw-Hill,  1956. 

[Ml]  MIRA,  C. , Etude  des  Systems  Asservis  non-lineaires . Toulouse 
University,  1968. 

[PI]  PAPOULIS , A.,  Probability,  random  variables  and  stochastic 
processes.  McGraw-Hill,  N.Y.,  1965. 

[51]  SHINBROT,  M.,  A generalization  of  a method  for  the  solution  of 

A 

the  integral  equation  arising  in  optimization  of  T.V.L.S. 
with  non-stationary  inputs.  Trans.  IRE,  IT-3,  1957. 

[52]  SIDI,  M.,  Synthesis  of  linear  feedback  systems  for  prescribed 

time  domain  tolerances.  Weizmann  Institute  Thesis,  1973. 

[53]  SEVELY,  Y.,  Asservissements  echantillonnees  lineaires. 

Toulouse  University,  1968. 

[54]  SCHWARTZ,  L. , Mathematics  for  the  Physical  Science.  Wesley, 

1966. 

[55]  S0L0D0V,  A.,  Linear  Automatic  control  systems  with  varying 

parameters . 

[S6J  STUBBERUD,  A.,  Analysis  and  synthesis  of  linear  time  varying 
systems.  University  of  California  Press,  1964. 

[Ul]  UTKIN,  V.,  Equations  of  the  slipping  regime  in  discontinuous 
synthesis  I.  Remote  & Contol,  p.  1157,  1972. 

[Ul]  UTKIN,  V.,  Equations  of  the  slipping  regime  in  discontinuous 
synthesis  II.  Remote  & Control,  p.  211,  1972. 

[Wl]  WIENER,  N.,  The  Extrapolation,  Interpolation  & Smoothing  of 
Stationary  Time  Series",  Wiley,  N.Y. , 1948. 


■ i"i 


REPORT  DOCUMENTATION  PAGE 

NIC  AD  INSTRUCTIONS 
m-'TOKK  COMPI.KTINC,  NORM 

1 REPORT  NUMBER  * f GOVT  ACCESSION  NO. 

AFCSR-7H-  7 7 - 12  2 4 ^ a — ^ 

3 RECIPIENT’S  CAT  ALOC  NUMBER 

4.  TITLE  (and  Subtitle) 

REDUCTION  OF  THE  COST  OF  FEEDBACK  IN 
SYSTEMS  WITH  LARGE  PARAMETER  UNCERTAINTIES 

s.  TYPE  OF  REPORT  6 PERIOD  COVERED 

I nterim 

r 

6.  PERFORMING  090.  REPORT  NUMBER 

7.  AuTHORfsj 

Patrick  Rosenbaum 

I saac  Horowitz 

8.  CONTRACT  OR  GRANT  NUMBERS; 

AFOSR  76-2946  T 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

University  of  Colorado  \ « 

L'iepartment  of  Electrical  Engineering 

Boulder,  Colorado  80309 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  d WORK  UNIT  NUMBERS 

61102F 

2304/A1 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Air  Force  Office  of  Scientific  Research/NM 
Bolling  AFB  DC  20332 


14.  MONITORING  AGENCY  NAME  & ADDRESSri/  different  from  Controlling  Office) 


13.  NUMBER  OF  PAGES 
200 


15.  SECURITY  CLASS,  (of  this  report j 

UNCLASSIFIED 


a.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


6.  DISTRIBUTION  STATEMENT  (of  this  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  from  Report) 


19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 


f0  ABSTRACT  7 C nntlntie  on  reverse  side  If  rwTTTTTr * ' >r) 

This  work  deals  with  ttie  synthesis  of  feedback  systems  to  achieve 
specified  performance  tolerances, despite  large  uncertainty  in  a con- 
strained  part  of  the  system,  denoted  as  the  plant’.  Part  of  this  work 
deals  with  linear  time  invariant  (lti)  plants  whore  the  ^*cost  of  feod- 


DD  1473  EOlTION  OF  I NOV  65  IS  OBSOLETE 

1INCI.ASS1  FI  ED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  f»>i tr  Car*  F.n farad) 


back'"^  if  lti  compensation  is  used,  is  primarily  in  the  bandwidth  of  ' 


1*03 


SECURITY  CLASSIFICATION  OF  THIS  P AGEfHTian  Cola  Entered) 


, 


20.  ABSTRACT  (Continued) 


the  feedback  loop  being  much  larger  than  that  of  the  system  as  a whole 
- making  the  system  very  sensitive  to  sensor  noise.  Here,  the  object ive 
is  to  reduce  the  loop  bandwidth  by  means  of  non-lti  compensation. 

The  result  is  a very 

significant  reduction  in  loop  bandwidth  and  with  it,  system  sensitivity 
to  sensor  noise. . - 


Another  method  of  non-lti  synthesis  is  by  means  of  linear  time- 
varying  compensation,  applicable  to  a certain  problem  class.  The 
solution  is  not  in  general  available  analytically,  but  is  found  for 
certain  cases  and  exhibits  in  these  reduced  'cost  of  feedback t"T~ 

”^he  plant  can  be  nonlinear  with  large  uncer- 
tainty giving  a set  of  nonlinear  plants  W . The  concept  is  use  of 
a set  P which  is  precisely  equivalent  to  W providing  the  output: 
is,  in  both  cases,  a member  of  an  acceptable  output  set  A 


UNCI.ASS  1 1;!  HI) 


